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THREE-DIMENSIONAL  EFFECTIVE  PROPERTY  AND  STRENGTH  PREDICTION 
OF  THICK  LAMINATED  COMPOSITE  MEDIA 

1  INTRODUCTION 

1.1  Objective  and  Motivation 

The  objective  of  this  work  was  to  develop  an  engineering  design  tool  capable  of  determining 
the  effective  three-dimensional  thermo-mechanical  properties  and  ultimate  strength  characteristics 
of  laminated  composite  media.  Specifically,  a  suitable  homogeneous  representation  is  sought  to 
replace  the  actual  heterogeneity  of  the  laminated  composite  material.  Homogeneous 
representations  in  the  form  of  three-dimensional  anisotropic  stiffiiess  matrices  and  the  typical 
thermo-mechanical  engineering  constants  of  laminated  media  are  developed.  In  addition,  a 
laminate  failure  analysis  that  identifies  specific  ply  failure  modes  and  provides  ultimate  laminate 
strength  predictions  (defined  in  terms  of  average  or  effective  laminate  stress)  has  been 
incorporated  in  the  predictive  capabilities. 

Modeling  ply-level  details  is  often  an  unrealistic  approach  in  the  design  and  analysis  of 
multi-layered  laimnated  composite  structures  because  of  computational  limitations  and  time 
constraints.  Consequently,  homogeneous  (or  “smeared”)  material  representations  for  “thick” 
laimnated  materials  are  often  highly  useful.  Such  representations,  when  used  in  conjunction  with 
traditional  structural  analysis  techmques  (e.g.,  the  finite  element  method),  enable  accurate 
prediction  of  the  effective  or  “average”  stress  and  strain  distributions  that  develop  within  a 

laminated  composite  structure  under  mechanical  and/or  thermal  loading  (Chou,  Carleone,  and  Hsu 
1972). 

The  failure  analysis  of  thick  laminated  composite  structures  should  not,  however,  be  based 
solely  on  average  stress  and/or  strain  distributions.  In  practice,  it  is  a  generally  accepted  fact  that 
any  realistic  laminate  failure  assessment  must  be  based  on  ply-level  stress  or  strain  results.  A 
reasonable  approximation  of  local  ply-level  stress-strain  results  can  be  obtained,  however,  by 
applying  the  average  stress-strain  state  as  a  local  set  of  boundary  conditions  onto  a  laminated 
media  boundary  value  problem.  In  this  manner,  ply-level  stress-strain  details  throughout  an 
entire  structure  can  be  computed  and  subsequently  employed  in  a  realistic  failure  assessment 
methodology. 


One  such  approach,  which  is  sometimes  referred  to  as  “smearing-unsmearing,”  is  depicted 
in  Figure  1 .  A  representative  sublaminate  configuration  for  the  composite  structure  is  first 
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identified  (see  Step  1  in  Figure  1).  In  the  most  general  sense,  this  material  could  be  considered  as 
heterogeneous  and  anisotropic.  A  set  of  equivalent  or  effective  homogeneous  properties  for  this 
representative  sublaminate  configuration  is  then  computed  (see  Step  2  in  Figure  1).  This  step  is 
referred  to  as  the  “smearing”  of  the  properties.  Although  this  is  typically  done  with  an  analytical 
solution  technique  (e.g.,  classical  laminated  plate  theory  [Whitney  1987]),  a  numerical  approach 
may  also  be  employed. 


Laminated  Composite 
Structure 


Appiy  Lamina 
Failure  Criteria 
and 

Postprocess 

Results 


Representative 
Sublaminate 
(Heterogeneous  & 
Anisotropic) 


Equivaient 
Representation 
(Homogeneous  & 
Anisotropic) 


Step  5 . 


I  Calculate  Ply-by-Ply 
Stress  and  Strain 


Structural  Analysis  to  Compute 
Stress  and  Strain  Distributions 

_ ^  ■ 

Figure  1.  Smearing-unsmearing  analysis  methodology. 

Next,  the  assumption  is  made  that  the  effective  or  homogeneous  thermo-mechanical 
properties  of  the  representative  sublaminate  configuration  can  be  used  to  replace  the  actual 
heterogeneity  of  the  laminated  material  in  the  structural  analysis.  A  typical  structmal  analysis  is 
then  conducted,  employing  the  effective  thermo-mechanical  properties  as  input  to  obtain  the 
average  stress  and  strain  distributions  within  the  structure  imder  a  prescribe  loading  (see  Step  3  in 
Figure  1).  It  is  important  to  note  that  the  stress  and  strain  distributions  here  represent  the 
average  values  (i.e.,  smooth  and  continuous),  since  the  laminated  material  has  been  replaced  with 
an  equivalent  set  of  homogeneous  properties.  They  do  not  reflect  the  actual  stress  and  strain 
distributions  (i.e.,  ply  by  ply)  associated  with  the  heterogeneity  of  the  material.  The  ply-level 
stress  and  strains  states  are  required  to  accurately  assess  failure  within  the  laminated  composite 
structure. 
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An  assessment  of  ply-level  stress  and  strain  results  can  be  obtained  straightforwardly  by 
applying  the  average  stress-strain  state  as  a  local  set  of  boundary  conditions  onto  a  laminated 
media  boundary  value  problem.  Again,  this  can  be  conducted  analytically  or  numerically  (see 
Step  4  in  Figure  1).  This  step  is  referred  to  as  the  “unsmearing”  of  the  average  stress-strain  state. 
Ultimately,  the  local  ply-level  stress-strain  states  throughout  the  structure  can  be  calculated  and 
employed  in  an  appropriate  lamina  failure  criterion  to  assess  overall  mechanical  performance  (see 
Step  5  in  Figure  1)  of  the  structure. 

A  major  result  of  this  work  is  the  codification  of  a  three-dimensional  laminated  media  model 
into  a  user-friendly  computer  program  environment  (LAM3D)  for  implementation  as  an 
engineering  design  and  analysis  tool  for  thick  laminated  composite  structures.  In  the  subsequent 
sections,  the  theoretical  backgroimd  of  the  laminated  media  theory  and  failure  analysis  employed 
in  LAM3D  are  presented.  Effective  thermo-mechanical  property  predictions  and  ultimate 
laminate  strength  estimates  from  LAM3D  are  validated  with  results  given  in  the  published 
literature  and  available  experimental  results.  Several  illustrative  example  problems  are  presented 
to  summarize  typical  input  and  output  parameters  required  by  LAM3D  and  to  demonstrate  the 
utility  of  the  code  for  general  laminate  analysis  purposes. 

1.2  Problem  Statement 

In  this  work,  the  analytical  model  developed  by  Chou,  Carleone,  and  Hsu  (1972)  is  used  to 
(1)  compute  the  effective  (homogeneous)  three-dimensional  stress-strain  constitutive  relations 
and  thermo-mechanical  properties  of  a  laminated  media  composed  of  an  arbitrary  number  of 
orthotropic  layers  or  plies,  (2)  compute  the  three-dimensional  ply-level  stress  and  strain  states 
within  the  laminated  media  subjected  to  three-dimensional  mechanical  loading  and/or  uniform 
thermal  loading,  and  (3)  generate  ultimate  laminate  strength  predictions  based  on  various  lamina 
failure  criteria  and  a  progressive  ply  failure  analysis  methodology. 

Three  different  lamina  constitutive  relations  are  considered:  (1)  a  general  orthotropic 
material  assumption  that  requires  nine  independent  elastic  constants,  (2)  a  transversely  isotropic 
material  assumption  that  requires  five  independent  elastic  constants,  and  (3)  a  simplified  form  of 
the  transversely  isotropic  case  requiring  only  three  independent  elastic  constants.  The  first  two 
material  constitutive  relations  require  in-plane  and  out-of-plane  material  properties  to  describe 
the  three-dimensional  stiffriess  matrix  for  a  single  lamina.  The  third  material  model  only  requires 
in-plane  properties  and  is  therefore  useful  in  practical  applications  when  three-dimensional 
material  properties  may  not  be  available  because  of  limited  experimental  characterization  data. 
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A  total  of  eight  lamina  failure  criteria  have  been  incorporated  into  the  failure  analysis  to 
predict  ultimate  laminate  strength  at  first  ply  failure  or  in  the  presence  of  progressive  ply  failure. 
Thermal  residual  ply  stresses  attributable  to  uniform  or  non-uniform  thermal  laminate  loading 
have  also  been  incorporated  in  the  strength  predictions.  Ultimate  laminate  strengths  as  w^ell  as 
the  particular  modes  of  failure  in  the  individual  plies  are  predicted. 

1.3  Literature  Survey 

Much  attention  has  been  given  to  determining  the  three-dimensional  effective  or 
homogeneous  properties  of  laminated  media  for  structural  analysis  purposes.  This  literature 
survey  presents  an  overview  of  relevant  research  conducted  in  this  area. 

Early  efforts  focused  on  obtaining  equivalent  material  stiffness  constants  for  laminated 
media  composed  of  individual  layers  that  are  isotropic  (White  and  Angona  1955;  Postma  1955; 
Rytov  1956;  Behrens  1967).  Salamon  (1968)  conducted  a  similar  analysis  but  computed 
effective  compliance  constants  for  the  laminated  media.  For  the  most  part,  the  approaches  taken 
in  these  early  works  employed  dispersion  techniques  and  elastic  wave  propagation  theories. 
These  researchers  foimd  the  effective  stiffiiess  constants  of  the  laminated  media  to  possess 
transverse  isotropy. 

Sun,  Achenbach,  and  Hermann  (1968)  have  also  investigated  the  problem  of  effective 
property  determination  of  laminated  materials.  They  have  taken  a  continuum  theory  approach  in 
which  the  microstructure  of  an  alternately  layered  medium  is  modeled.  Achenbach  (1970)  later 
showed  how  the  results  of  the  continuum  theory  approach  reduce  to  earlier  generated  results 
under  appropriate  assumptions  related  to  the  details  of  the  microstructure. 

Christensen  (1988)  developed  a  three-dimensional  lamina  constitutive  relationship  based  on 
only  three  independent  properties  rather  than  the  five  typically  associated  with  transversely 
isotropic  material  models.  He  used  this  simplified  lamina  constitutive  theory  to  develop  a  three- 
dimensional  laminate  theory  that  closely  resembles  the  classical  two-dimensional  laminated  plate 
theory  for  the  in-plane  directions.  Out-of-plane  directional  properties  for  the  laminate  are 
identical  to  those  of  any  individual  layer  since  they  are  assumed  to  be  independent  of  in-plane 
ply  orientations.  These  results  led  directly  to  a  three-dimensional  lamina  failure  theory  that 
separates  fiber  and  matrix  dominated  failure  modes.  The  composite  lamina  failure  theory  closely 
resembles  the  Von  Mises  criterion  often  used  for  isotropic  materials. 
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Enie  and  Rizzo  (1970)  developed  a  convenient  set  of  expressions  for  determining  the 
effective  three-dimensional  mechanical  properties  of  balanced  laminates.  They  employed 
classical  laminated  plate  theory  for  calculating  in-plane  properties.  Out-of-plane  properties  \vere 
predicted,  based  on  uniform  interlaminar  shear  and  interlaminar  normal  stress  distributions  and 
uniform  in-plane  strains  assumptions.  Model  predictions  were  foimd  to  be  in  reasonable 
agreement  with  finite  element  results  for  several  different  composite  laminate  configurations. 

Pagano  (1974a,  1974b)  developed  a  three-dimensional  laminate  model  within  the  framework 
of  anisotropic  heterogeneous  elasticity  theory  to  predict  the  effective  response  of  composite 
laminates  to  mechanical  and  thermal  loading.  His  work,  which  was  not  restricted  to  isotropic 
layers  or  the  assumption  of  a  repeating  sublaminate,  considered  the  representative  volume  to 
contain  the  entire  thickness  of  the  laminate.  Laminate  symmetry  was  also  not  required,  and 
consequently,  his  constitutive  “stress-strain”  relations  included  moment  and  curvature  effects. 

Chou,  Carleone,  and  Hsu  (1972)  conducted  a  control  volume  approach  to  yield  a  closed 
form  solution  to  the  problem  of  effective  homogeneous  property  determination  for  a  laminated 
media  composed  of  individual  layers.  Unlike  the  works  of  White  and  Angona  (1955),  Postma 
(1955),  Rytov  (1956),  Behrens  (1967)  and  Salamon  (1968),  which  required  the  individual  layers 
to  be  isotropic,  Chou,  Carleone,  and  Hsu  (1972)  permitted  general  anisotropy  of  the  layers.  The 
analysis  is  based  on  the  assumptions  that  all  interlaminar  stresses  are  continuous  across  ply 
interfaces  and  that  all  in-plane  strains  are  continuous  through  the  thickness  dimension  of  the 
representative  volume  element  (i.e.,  the  repeating  sub-laminate  configuration). 

Sun  and  Li  (1988)  developed  closed  form  expressions  to  predict  the  equivalent 
homogeneous  anisotropic  representation  of  a  laminated  composite  medium.  Their  analysis  is 
remarkably  similar  to  the  earlier  work  presented  by  Chou,  Carleone,  and  Hsu  (1972).  Reduced 
expressions  for  the  elastic  constants  are  derived,  which  greatly  simplify  the  analysis.  Effective 
property  predictions  for  several  laminates  are  generated  and  found  to  be  in  reasonable  agreement 
ivith  the  results  of  Enie  and  Rizzo  (1970). 

Roy  and  Tsai  (1992)  formulated  a  boundary  value  problem  approach  and  used  the  Airy 
stress  function  solution  technique  to  estimate  three-dimensional  effective  properties  for 
laminated  composite  media.  They  note  that  since  the  works  of  Pagano  (1974a,  1974b);  Chou, 
Carleone,  and  Hsu  (1972);  and  Sun  and  Li  (1988)  assume  the  distributions  of  interlaminar  shear 
stress  to  be  constant  through  the  thickness  of  the  laminate,  their  effective  interlaminar  shear 
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stiffiiesses  are  independent  of  stacking  sequence,  which  is  not  in  complete  agreement  with  the 
experimental  observations  of  Roy  and  Kim  ( 1 989). 

The  model  developed  by  Chou,  Carleone,  and  Hsu  (1972)  was  selected  for  implementation 
in  the  present  work,  based  on  the  combined  merits  of  its  accuracy  and  computational  efficiency. 
Their  model  predictions  have  been  validated  with  the  results  of  other  investigators.  In  addition, 
their  model  provides  a  convenient  closed  form  solution  for  enhanced  computational  efficiency, 
which  is  a  strong  consideration  in  future  work  for  the  development  of  an  engineering  tool  for  the 
design  and  analysis  of  thick  laminated  composite  structures  (e.g.,  finite  element  pre-  and  post¬ 
processors). 

2.  ANALYSIS 

2.1.  Nomenclature 

The  laminate  nomenclature  adopted  in  this  work,  including  ply  stress  and  ply  strain 
notation,  is  consistent  with  that  in  most  composite  material  primers  for  two-dimensional 
laminate  analyses  (Whimey,  Daniel,  and  Pipes  1982;  Whitney  1987).  The  three-dimensional 
version  of  this  notation  is  reviewed  in  this  section. 

An  orthogonal  principal  material  coordinate  system  (1, 2, 3)  of  a  single  ply,  or  lamina,  is 
illustrated  in  Figure  2.  The  principal  material  directions  include  the  fiber  or  longitudinal  (1- 
direction),  the  transverse  (2-direction)  and  the  normal  (3-direction).  A  laminate  or  representative 
sublaminate  configuration  is  composed  of  N  individual  plies  and  is  defined  with  respect  to  an 
orthogonal  (x,  y,  z)  global  laminate  coordinate  system  (see  Figure  3).  The  (1,2)  plane  of  each  ply 
is  assumed  to  be  co-planar  with  the  x-y  laminate  plane.  The  orientation  angle  of  each  ply,  0, 
defines  the  right-hand  rotation  of  the  (1, 2, 3)  ply  coordinate  system  about  the  3  axis  with 
respect  to  the  fixed  (x,  y,  z)  system  (see  Figure  3).  The  “stacking  sequence”  or  “lay-up”  of  the 
laminate  describes  the  orientation  pattern  of  all  plies  in  the  laminate.  The  geometric  mid-plane  of 
the  laminate  is  positioned  at  the  z=0  location.  The  bottom  ply  is  designated  as  Ply  1  and  the  top 
ply  is  designated  as  Ply  N.  A  notation  is  used  to  define  the  lay-up  in  which  ply  angles  are 
separated  by  slashes,  and  brackets  are  used  to  enclose  all  the  plies  in  the  laminate.  For  example, 
the  nomenclature  [45/90/0/-45/-45/0/90/45]  defines  stacking  sequence  for  the  eight-ply  laminate 
depicted  in  Figure  3. 
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3-  Normal 


Figure  2.  Lamina  (principal  material)  coordinate  system. 


Figures.  Laminate  nomenclature. 


When  stress  and  strain  are  referenced  in  the  principal  (1, 2,  3)  material  coordinate  system, 
they  are  denoted  with  an  accent  and  are  represented  according  to  the  following  conventions.  For 
stress,  the  notation  is  given  by 

cr’i=  <722  <733  <723  <713  <712]  for  (1) 

Similarly,  the  notation  for  strain  is  given  by 

e’i=  [e’ll  e’22  e’33  e’23  e’l3  £12]^  for  (i=l,6)  (2) 

The  subscripts  1 1, 22,  and  12  refer  to  the  longitudinal,  transverse  and  shear  components, 
respectively,  of  both  in-plane  stress  and  strain.  The  subscripts  33, 23,  and  13  refer  to  the 
interlaminar  normal  and  two  interlaminar  shear  components,  respectively,  of  the  out-of-plane 
stress  and  strain.  This  contracted  notation  is  schematically  illustrated  in  Figure  4.  for  the  stresses 
acting  on  the  faces  of  an  idealized  material  element  of  a  unidirectional  ply.  The  notation  for  ply 
strains  would  be  represented  in  an  analogous  manner. 

When  stress  and  strain  are  referenced  within  the  global  (x,  y,  z)  laminate  coordinate  system, 
they  are  denoted  with  a  bar.  Stresses  are  expressed  according  to 


Figure  4.  Contracted  notation  for  stress. 
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for  (i=l,6) 


(3) 


Cy  Cz  O'yz  O’xy] 

and  strains  by 

£  i  “  [  £  X  c  y  c  z  £  yz  £  xz  £  xy  ]  for  0  ~  1 »  6)  (4) 

The  subscripts  x,  y,  and  xy  represent  the  in-plane  components  of  stress  and  strain  defined  in  the 
global  (x,  y,  z)  laminate  coordinate  system.  The  subscripts  z,  yz,  and  xz  represent  the  out-of¬ 
plane  components  of  stress  and  strain. 

2.2  Lamina  Constitutive  Relations 

The  three-dimensional  Hooke’s  Law  linear-elastic  stress-strain  constitutive  relation  for  an 
individual  lamina  is  written  in  the  following  contracted  form 

a’i  =  C’ije’j  for  (i,j  =  1, 2, 3, 4, 5, 6)  (5) 

in  which  C’ij  represents  the  stiffiiess  matrix  of  the  material  defined  in  the  principal  material 
coordinate  system.  For  a  material  that  possesses  general  anisotropy,  21  independent  elastic 
constants  are  required  to  characterize  the  stiffiiess  matrix.  Various  assiunptions  about  the 
material  symmetry  of  continuous  fiber-reinforced  composites  can  be  made  to  reduce  the  number 
of  independent  constants  required  to  characterize  the  lamina  stiffiiess  matrix.  The  three  types  of 
lamina  constitutive  models  considered  in  this  work,  which  are  simplifications  of  the  general 
anisotropic  case,  are  summarized  in  the  following  sections. 

2.2.1  Orthotropic  Material  Model 

An  individual  lamina  of  unidirectional  composite  material  possess  three  mutually 
orthogonal  planes  of  elastic  symmetiy  (see  Figure  2)  and  is  said  to  be  “orthotropic.”  The 
orthotropic  lamina  constitutive  model  is  the  most  general  type  of  material  model  considered  in 
this  work.  The  elastic  symmetry  of  orthotropic  materials  reduces  the  number  of  independent 
elastic  constants  required  to  characterize  the  lamina  stiffiiess  matrix  from  21  to  9.  These  elastic 
constants  are  defined  in  terms  of  the  following  typical  engineering  properties:  the  extensional 
moduli  El,  E2,  and  E3;  the  shear  moduli  G23,  G13,  and  G12;  and  the  Poisson’s  ratios  V23,  V13,  and 
V12.  The  orthotropic  lamina  stiffiiess  matrix  is  symmetrical  (i.e.,  C’lj  =  C  ji  for  i,  j  =  1, 2, 3, 4, 5, 
6)  and  takes  the  following  form  (Whitney,  1987): 
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C'li  C’i2  Ci3  0  0  0 

C'l2  C'22  C23  0  0  0 

C'i3  C'23  C’33  0  0  0 

0  0  0  C'44  0  0 

0  0  0  0  C'55  0 

0  0  0  0  0  C'66 

The  non-zero  stiffness  coefficients  of  the  orthotropic  lamina  stiffness  matrix  are  defined  as 
follow; 

C'li  =  (1-v||E3/E2)Ei/V 
C'i2  =  (V12  +  V13  V23E3/E2)E2/V 
C'l3  =  (Vi3 +Vi2V23)E3/V 

C'22  =  (1-Vi3E3/Ei)E2/V 

C'23  =  (V23  -  V12V13  E2/Ei)E3/V  (7) 

C'33  =  (1-Vi2E2/Ei)E3/V 
C'44  =  G23 

C'55  =  Gi3 
C'66  =  Gi2 

in  which 

V  =  1  -  Vi2(Vi2E2/Ei +2  V23  V13E3/E1)  -  Vi3  E3/E1  -  V23  E3/E2 

2.2.2  Transversely  Isotropic  Material  Model 

It  is  generally  acceptable  to  assume  that  the  composite  lamina  also  possess  material 
symmetry  in  the  (2, 3)  plane  of  the  principal  material  coordinate  system  (see  Figure  2).  The 
material  is  therefore  said  to  be  “transversely  isotropic”  and  the  following  conditions  are  assumed 
to  hold  true: 

E3  =  E2 

Gi3  =  Gi2  (8) 

Vl3=Vi2 


10 


In  this  case,  the  number  of  independent  elastic  constants  is  reduced  from  nine  (orthotropic)  to 
five  (transversely  isotropic).  The  lamina  stiffiiess  matrix  for  the  transversely  isotropic  material 
model  takes  the  same  form  as  that  given  in  equation  (6).  Using  a  subset  of  the  coefficients 
defined  in  equation  (7),  the  matrix  is  written  as  (Whimey,  1987) 


C’li  C’i2  Ci2  0  0  0 

C'i2  C'22  C23  0  0  0 

C'i2  C'23  C'22  0  0  0 

0  0  0  0  0 
0  0  0  0  C'66  0 

-  0  0  0  0  0  C'66- 


(9) 


Now  the  reduced  set  of  engineering  constants  required  to  characterize  the  lamina  stiffiiess  matrix 
includes  only  Ei,  E2,  G12,  V12,  and  V23. 


2.2.3  Simplified  Transversely  Isotropic  Material  Model 

A  further  reduction  in  the  number  of  elastic  constants  required  to  characterize  the  lamina 
stiffiiess  matrix  is  also  possible  (Christensen,  1 988).  Christensen  has  suggested  that  the  number 
of  elastic  constants  may  be  reduced  to  only  three,  using  the  following  two  restrictions  on  the 
lamina  stiffiiess  matrix  coefficients 


C'22-C'23_ 

2 


C'66  and  C'23  =  C'i2 


(10) 


In  terms  of  the  engineering  properties,  these  restrictions  or  conditions  are  met  if  the  following 
relationship  is  found  true: 


Gi2= 


{l-ni2}  E2 


2{1- 


ni22  E2 
El 


} 


(11) 
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Christenson  (1988)  showed,  with  experimental  evidence,  that  these  assumptions  are  valid  for 
several  typical  composite  material  systems.  If  Equation  (11)  applies,  the  lamina  stiffness  matrix 
assumes  the  following  form: 

C'li  Ci2  C’i2  0  0  0 

C'l2  C'22  C'i2  0  0  0 

C'i2  C'i2  C22  0  0  0 

0  0  0  C'66  0  0 

0  0  0  0  C'66  0 

0  0  0  0  0  C'66 

in  which  only  the  in-plane  engineering  properties  Ei,  E2,  and  V12  are  required  to  characterize  the 
lamina  stiffiiess  matrix.  Motivation  for  this  material  model  is  because  empirical  characterization 
of  the  in-plane  material  properties  is  significantly  more  straightforward  than  that  of  the  out-of- 
plane  properties. 


rransformations  Between  the  Principal  and  Laminate  Coordinate  Svstems 


Transformation  of  ply  stress  and  strain  between  the  principal  (1, 2, 3)  and  global  (x,y,z) 
coordinate  systems  is  required.  The  global  ply  stresses,  o  i,  can  be  expressed  explicitly  in  terms 
of  the  principal  ply  stresses,  a’i,  and  the  ply  orientation  angle,  6,  (see  Figure  3). 

Mathematically,  this  transformation  is  accomplished  with  the  following  second  order  tensor 
transformation: 


ffi=  [T(e)]ja  (13) 

in  which  the  stress  transformation  matrix  is  given  by 

m2  n^  0  0  0  2mn 

rP-  m2  0  0  0  -2mn 

0  0  1  0  0  0 

0  0  0  m  -n  0 

0  0  0  n  m  0 

-mn  mn  0  0  0  (m2  -  n- 


and  m  =  cos  6,  n  =  sin  6.  Similarly,  global  ply  strains  are  obtained  according  to 


ei=  [T(9)]fje’ 


(15) 


in  which  the  strain  transformation  matrix  is  given  by 


[T(e)]  y  = 


”  m2  n2  0  0  0  -mn  ~ 
n2  m2  0  0  0  mn 

0  0  1  0  0  0 

0  0  0  m  -n  0 

0  0  0  n  m  0 

— 2mn  2mn  0  0  0  (m2  -  n2)- 


(16) 


The  lamina  stress-strain  constitutive  relationship,  defined  in  the  global  (x,y,z)  laminate 
coordinate  system,  is  written  explicitly  as 


CTi  =  Cijej  for  (i,j  =  l,2,3,4,5,6)  (17) 

in  which  the  bars  denote  that  the  quantities  are  defined  in  the  global  coordinate  system. 

Combining  equations  (5),  (13),  (15),  and  (17),  it  can  be  shown  that  the  lamina  stiffness  matrix, 

Cy,  can  be  expressed  explicitly  in  terms  of  the  principal  lamina  stiffness  matrix  and  the  ply 
orientation  angle,  0,  through  the  following  expression: 

Cij=  [1(0)]  [j  C'ij  {[T(0)]^}-1  (18) 

2.4  Effective  Laminated  Media  Constitutive  Relations 

The  theoretical  background  of  a  stress-strain  constitutive  model  for  predicting  the  effective 
material  constants  of  laminated  media  is  summarized  in  this  section.  The  detailed  theoretical 
development  of  this  model  is  presented  elsewhere  (Chou,  Carleone,  and  Hsu  1972),  and 
therefore,  only  the  most  significant  features  are  reviewed  here. 

First,  an  N  layered  laminate  (i.e.,  representative  sublaminate  configuration),  such  as  the  one 
depicted  in  Figure  3,  is  identified.  The  “laminate”  is  considered  to  be  representative  of  a  “small 
element”  (i.e.,  material  point)  in  a  larger  structural  body.  The  objective  is  to  obtain  a  single 
constitutive  relationship  that  can  be  used  to  define  the  effective  (i.e.,  homogeneous)  stress-strain 
response  of  the  laminate.  The  following  expression  is  used  to  represent  the  effective  stress-strain 
constitutive  relationship  for  the  laminate 
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ai*  =  Cij*ej* 


for  (i,j  =  1,2,3,4,5,6) 


(19) 


The  superscript  is  used  here  to  denote  the  “average”  or  effective  laminate  stress  and  strain 
quantities. 

In-plane  strains  are  assumed  to  be  uniform  (i.e.,  constant  within  each  ply)  and  equal  to  the 
effective  strains  of  the  laminate.  Mathematically,  this  is  expressed  as 

£  ik  =  e*  for  (i  =  1, 2, 6;  k  =  1, 2, ...,  N)  (20) 

in  which  £  represents  the  strain  in  the  k*  ply  of  the  laminate  (see  ply-numbering  convention  in 
Figure  3). 

To  ensure  stress  continuity  across  ply  interfaces,  all  ply  stress  components  associated  with 
the  out-of-plane  direction  (i.e.,  z-direction)  are  assumed  to  be  uniform  and  equal  to  the 
corresponding  effective  stresses  in  the  laminate.  Mathematically,  this  is  expressed  as 

o-^=oC  for  (i  =  3,4,5;k=l,2,...,N)  (21) 

in  which  ct  jk  represents  the  stress  in  the  k^k  piy  of  the  laminate. 

All  remaining  effective  laminate  stresses  and  strains  are  assumed  to  be  the  volume  average 
of  all  their  corresponding  ply  stress  and  strain  components,  respectively.  Mathematically,  these 
assumptions  are  expressed  as 


N 


I  vk  Cjk  for 
k=l 

(i  =  3,4,5) 

(22) 

N 

I  vk  CTjk  for 
k=l 

(i=l,2,6) 

(23) 

in  which  Vk  is  the  ratio  of  the  original  (i.e.,  undeformed)  volume  of  the  k*^k  piy  over  the  original 
volume  of  the  entire  laminate. 

Equation  (17)  is  re-written  below  using  the  superscript  notation  to  emphasize  that  it 
applies  to  each  individual  ply  in  the  laminate: 

aik  =  Cijk  ejk  for  ( i,  j  =  1, 2, 3, 4,  5, 6;  k=l,2, ...,  N)  (24) 
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Equations  (19)  through  (24)  represent  12N+6  linear  algebraic  equations  in  12N+12  unknowns. 
Chou,  Carleone,  and  Hsu  (1972)  found  that  the  solution  to  Equations  (19)  through  (24)  yields  the 
following  effective  three-dimensional  stress-strain  constitutive  relation  which  can  be  used  as  an 
equivalent  (i.e.,  homogeneous)  representation  for  the  laminated  media  where  the  coefficients  in 
the  laminate  stiffness  matrix,  Cij*,  are  given  by 


N 

Cij*=  IVk 
k=l 


Cijk- 


C33L 


CnkCsjk  ^ _ 

Csjk  Cssk  X51i 


C33^  -I 


for  (i,j  =  1,2,3,  6) 


(25) 


Cij*  =Cji*  =0  for  (i  =  l,2,3,6;j  =  4,5) 
and 


C-* 


vN  yk  o  k 

^k=l  Cij 


Sk=l  (^44*^  Cssk  -  C45k  C54k) 


for  (i,j  =  4,5) 


(26) 


(27) 


in  which 


C44k  C45k 

C54k  Cssk 


C44k  Cssk  -  C45k  C54k 


(28) 


The  effective  stress-strain  constitutive  relation  for  the  laminated  media  is  therefore  given  by 
Equations  (19)  and  (25)  through  (28).  The  effective  laminate  stiffiiess  matrix  described  is 
included  in  the  LAM3D  code  predictions. 

2.5  Effective  Laminated  Media  Engineering  Constants 

In  the  following  sections,  the  constitutive  relation  just  presented  will  be  used  to  define  the 
effective  thermo-mechanical  properties  (i.e.,  typical  engineering  constants)  for  the  laminated 
media.  These  quantities  are  often  viewed  as  a  more  convenient  form  of  input  for  structural 
analyses  purposes  (e.g.,  finite  element  methods). 
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2.5.1  Mechanical  Properties 

A  set  of  effective  three-dimensional  engineering  constants,  which  can  be  used  to  represent 
the  laminated  media  as  a  homogeneous  continuum,  is  now  presented.  The  engineering  constants, 
defined  in  the  global  laminate  coordinate  system,  include  the  extensional  young’s  moduli.  Ex,  Ey, 
and  Ez;  the  shear  moduli,  Gyz,  Gxz,  and  Gxy;  the  major  Poisson’s  ratios,  Vyz,  Vxz,  Vxy;  and  the 
minor  Poisson’s  ratios,  Vzy,  Vzx>  Vyx  (see  Figure  3  for  coordinate  reference).  Note  that  the 
Poisson’s  ratios  presented  here  are  defined  in  accordance  with  typical  conventions  (i.e.,  Vij  =  -£j  / 

Ei). 


The  constitutive  relation  equation  (19)  (which  is  in  the  “stiffness”  form)  is  first  re¬ 
expressed  in  its  “compliance”  form.  Inverting  both  sides  of  equation  (19),  the  effective  laminate 
compliance  relationship  becomes 


Ei*=Hij*Oj*  for  (i,j  =  1,2,3,4,5, 6) 

(29) 

in  which  Hjj*  is  termed  the  laminate  compliance  matrix,  simply  the  inverse  of  the  laminate 
stiffness  matrix  (i.e.,  H  jj*  =  Cij*"*).  From  this  point,  the  effective  mechanical  engineering 

constants  can  be  expressed  explicitly  in  terms  of  the  compliance  coefficients,  Hy* 
extensional  moduli  are  given  by 

.  The  effective 

Ex  =  1/Hii* 

(30) 

Ey  =  I/H22* 

(31) 

and 

Ez=  I/H33* 

(32) 

The  effective  shear  moduli  are  given  by 

Gyz=  I/H44* 

(33) 

Gxz=  I/H55* 

(34) 

and 

Gxy=  1/H66* 

(35) 

The  effective  major  Poisson’s  ratios  are  given  by 

Vyz  =  -  H23*  /H  22* 

(36) 

Vxz  =  -  H13VH11* 

(37) 
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(38) 


Vxy=  -Hi2*/Hii* 

and  the  effective  minor  Poisson  ratios  are  given  by 


^zy  ■“  ■ 

H23VH33* 

(39) 

Vzx  “  • 

H13VH33* 

(40) 

Vyx  =  - 

H12VH22* 

(41) 

The  effective  mechanical  properties  just  described  are  included  in  the  LAM3D  code  predictions. 
2.5.2  Thermal  Expansion  Coefficients 

Effective  thermal  expansion  coefficients  can  also  be  derived  for  the  laminated  media.  The 
principal  direction  thermal  expansion  coefficients  for  the  individual  lamina  are  first  rotated  in  the 
global  laminate  coordinate  system.  The  strain  rotation  matrix  is  used  (see  Equation  (16)): 

ai=  [T(0)]ya’  (42) 

in  which,  again,  the  bar  notation  denotes  the  laminate  coordinate  system  and  the  prime  denotes 
the  principal  material  coordinate  system.  A  quantity  that  essentially  represents  a  non- 
dimensional  (average)  uniform  thermal  loading  condition  on  the  laminate  is  first  computed 
according  to 


Njth*  =  I  CijkajkVk 
k=l 


for  (i,j  =  l,2,3,4,5,6) 


(43) 


The  effective  laminated  media  thermal  expansion  coefficients,  defined  in  the  global  laminate 
coordinate  system,  are  then  given  by 


ax*  = 

=  H,j*  Njft* 

for 

0  = 

1,. 

..,6) 

(44) 

ay*  = 

■  H2j’  Njth  • 

for 

0  = 

1,. 

..,6) 

(45) 

PI 

N 

II 

’H3j*  Nj*’ 

for 

0  = 

1,. 

..,6) 

(46) 

OCyz* 

=  H4j’  Njtl>  • 

for 

0  = 

1,. 

..,6) 

(47) 

Otxz 

=  H5j*  Nj®* 

for 

0  = 

1,. 

..,6) 

(48) 

ttxy 

=  Hsj*  Nj*  • 

for 

0  = 

1,. 

..,6) 

(49) 
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The  effective  laminate  thennal  expansion  coefficients  just  described  are  included  in  the  LAM3D 
code  predictions. 

2.6  Plv-Level  Stress  and  Strain  Calculations 

To  assess  failure  within  the  laminated  media  because  of  mechanical  and/or  thermal  loading, 
all  ply-level  stress  and  strains  values  must  be  computed.  The  following  sections  describe  this 
procedure. 

2.6.1  Applied  Mechanical  Loading 

The  assumption  here  is  that  the  applied  mechanical  loading  on  the  laminated  media  ( a  i*)  is 
known,  uniform  and  represents  the  “average”  or  “effective”  stress  acting  on  the  sub-laminate 
configuration.  The  associated  “effective”  or  “smeared”  laminate  strains  (ej*)  can  be  obtained 
directly  from  Equation  (29).  From  the  assumption  made  in  Equation  (20),  all  in-plane  strain 
values  (defined  in  the  global  (x,  y,  z)  coordinate  system)  for  Plies  1  through  N  are  therefore 
known.  Similarly,  from  the  assumption  made  in  Equation  (21),  all  out-of-plane  stresses  for  Plies 
1  through  N  are  known  (also  defined  in  the  global  (x,  y,  z)  coordinate  system).  The  out-of-plane 
ply  strains  and  in-plane  ply  stresses  remain  to  be  determined. 

Sim  and  Liao  (1990)  derived  the  following  expression  for  determining  the  remaining  out-of- 
plane  ply  strains: 


Cssk  C34k  Cssk 

-1 

03k 

C3ikC32*^  C36k 

£lk 

= 

C43k  C44k  C45k 

04k 

- 

041^042^046^ 

£2*^ 

-C53k  C54k  Cssk  - 

-Gfk- 

-C5ikC52’^  056*^- 

-£6^_ 

Once  all  the  ply  strains  are  known,  the  remaining  in-plane  ply  stresses  can  be  calculated 
straightforwardly  through  the  following  relation: 
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Ciik  Ci2k  Cisk  Ci4k  Cis^  Cie^ 

02k 

= 

C2lk  €22^  €23^  024^  C25k  €26^ 

-aek- 

”C6lk  €62^  €53^  064^  Cesk  - 

£3*^ 

£4*^ 

£5^ 


(51) 


The  mechanical  ply  stress  and  strains  just  described  are  included  in  the  LAM3D  code 
predictions. 


2.6.2  Thermal  Loading 

Ply  stresses  and  ply  strains  resulting  from  thermal  loading  are  independent  of  the 
mechanical  ply  stresses  and  strains  but  are  computed  in  an  analogous  manner.  Thermal  ply 
stresses  result  from  the  fact  that  the  individual  plies  within  the  laminate  have  different  thermal 
expansion  coefficients  in  different  directions  but  are  constrained  to  deform  or  strain  uniformly. 
To  determine  thermal  ply  stresses  and  strains,  the  effective  or  average  thermal  stress  for  the 
laminated  media  is  first  computed  for  a  given  prescribed  thermal  loading.  In  this  formulation,  the 
thermal  load  for  each  ply  in  the  laminated  media  (i.e.,  AT*^)  can  be  different.  The  effective 
thermal  stress  acting  on  the  laminate  is  given  by 

U  *  N  _ 

(ap)  =  S  CijkajkykATk  for  (i,j  =  1,2,3,4,5,6)  (52) 

k=l 


The  actual  thermal  strain  or  deformation  that  the  laminate  undergoes,  (ei*  )*,  is  computed  in  a 
similar  manner  to  the  mechanical  strains,  through  the  following  relationship: 

(eith)*=  Hij*(Ojth)*  for  (i,j  =  l,2,3,4,5,6)  (53) 

As  with  the  calculation  of  the  mechanical  ply  strains,  the  in-plane  thermal  ply  strains  are 
assumed  to  be  constant  and  equal  to  the  actual  in-plane  thermal  laminate  strains: 

(eith  )k  =  (e -th )*  for  (i  =  l,  2, 6;  k  =  1, 2, ..., N)  (54) 
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in  which  (£  represents  the  thermal  strain  in  the  k*  ply.  The  out-of-plane  thermal  ply 
stresses  are  assumed  to  be  constant  and  equal  to  the  effective  out-of-plane  thermal  laminate 
stress 

(Oith)k  =  (ojth)*  for  (i  =  3,4,5;k=l,2,...,N)  (55) 

The  remaining  out-of-plane  thermal  ply  strains  are  then  straightforwardly  computed  through  a 
relationship  similar  to  that  given  in  Equation  (50): 

(is'h)*^  C33kC34kC35‘'  (^3*)*  €31^032*^03  5!^  (il*)* 

(i4^h)k  =  C43*'C44*'C45‘^  (^4*)*  ‘  041^042*^046^  (£2**^)* 

_(e5th)k_  Lc53kC54kC55k-l - (a5th)*J  ^C5i*^C52*C56*‘-L(  £5*  )* 

At  this  point,  the  stress-free  expansional  strains  (i.e.,  the  strains  that  would  be  induced  in 
each  ply  if  they  were  not  constrained  to  deform  together)  are  subtracted  from  the  actual  laminate 
thermal  strains.  The  resulting  thermal  ply  strains  are  actually  those  that  will  be  used  to  compute 
the  remaining,  unknown,  in-plane  thermal  ply  stresses  and  are  expressed  according  to 

(ejth  )k  =  (ejth  )k .  a jk  at  k  for  ( i  =  1, 2, 3, 4,  5, 6;  k=l,2, ...,  n)  (57) 


The  out-of-plane  thermal  ply  stresses  are  given  in  Equation  (55),  and  the  remaining  in-plane 
thermal  ply  stresses  are  given  by 


OiikOi2*^Oi3*^Oi4*^Oi5*^Oi6*^ 
021*^  022*^  023*^  024*^  €25*^  026*^ 
06 1*'  062*^  063*^  064*^  065*^  066*^ 


(eith)k 

(e2th)k 

(i3th)k 

(  £4*^  )k 
(e5th)k 

■(e6*)*'' 


(58) 


The  thermal  ply  stress  and  strains  just  described  are  included  in  the  LAM3D  code  predictions. 
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2.7  Laminated  Media  Strength  Predictions 


The  overall  failure  assessment  of  the  laminated  media  or  laminate  is  based  on  the  stress 
and/or  strain  states  that  exist  in  the  individual  plies  as  well  as  ply-level-based  (lamina)  composite 
failure  criterion  and  some  type  of  failure  solution  strategy.  For  a  given  set  of  mechanical  and/or 
thermal  loads  imposed  on  the  laminate,  the  individual  ply  stresses  and  strains  are  first  calculated 
according  to  the  procedure  outlined  in  the  previous  sections.  Any  one  of  a  number  of  lamina 
failure  criteria  can  then  be  used  to  make  a  failure  assessment  in  all  the  individual  plies  within  the 
laminate.  An  overall  laminate  strength  prediction  is  then  made  by  considering  the  ply-level 
failure  assessment  of  all  plies  in  conjunction  with  some  type  of  failure  solution  strategy.  A 
failure  solution  strategy  is  simply  an  assumption  or  set  of  assumptions  that  define  how  load  is 
redistributed  fi’om  the  failed  plies  to  the  other  plies  within  the  laminate. 

In  the  following  discussion,  the  specific  lamina  failure  criteria  used  in  LAM3D  are 
described.  Various  failure  solution  strategies  also  employed  in  LAM3D  are  then  presented  (i.e., 
first  ply  and  progressive  ply  failure  solution  strategies).  The  section  closes  by  describing  the 
manner  in  which  overall  laminate  strength  predictions  (quantitative  and  qualitative  laminatp: 
failure  assessments)  are  presented  by  LAM3D. 

2.7.1  Lamina  Failure  Criteria 

Several  articles  review  the  failure  theories  for  composite  materials.  Nahas  (1986)  gives  an 
overview  of  failure  theories  for  composites  and  isotropic  materials  proposed  before  1986.  Hahn 
and  Kallas  (1992)  have  presented  a  detailed  analysis  of  several  more  recent  composite  failure 
theories.  This  report  will  briefly  review  the  eight  lamina  failure  criteria  that  have  been  included  in 
the  LAM3D  code. 

2.7. 1.1  Von  Mises-Hencky  (isotropic  materials  only) 

The  Von  Mises-Hencky  or  maximum  distortion  energy  theory  predicts  that  a  ductile 
material  under  a  general  stress  state  will  yield  when  its  shear  distortional  energy  (the  total  strain 
energy  minus  the  strain  energy  attributable  to  change  in  volume)  is  equal  to  the  shear  distortional 
energy  imder  simple  tension  (Hertzberg,  1989;  Nahas,  1986).  This  yield  theory  is  only  valid  for 
isotropic  materials  and,  consequently,  is  not  generally  appropriate  for  composite  materials.  It 
has  shown  good  correlation  with  test  data  for  metals  under  multi-axial  loading  since  it  includes  the 
interactive  effects  of  all  the  stress  components.  The  Von  Mises-Hencky  theory  is  applied  by 
first  calculating  the  equivalent  or  effective  stress  ( crj  acting  on  a  material  element  (Hertzberg, 

1989)  according  to 
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(59) 


This  failure  theory  predicts  that  the  material  will  fail  when  the  effective  stress  becomes  greater 
than  the  material  yield  strength,  <T,.,  or  when 

o>a^..  (60) 


2.7. 1.2  Maximum  Stress 

Perhaps  the  most  widely  used  failure  criterion  for  unidirectional  composites  is  the 
maximum  stress  failure  criterion  (Tsai,  1987;  Nahas,  1986),  which  predicts  that  a  material  will  fail 
when  the  magnitude  of  the  stress  in  any  direction  exceeds  its  corresponding  allowable  level  in  that 
direction.  This  criterion  is  valid  for  both  isotropic  and  anisotropic  materials.  However,  it  does 
not  consider  interactions  between  the  various  stress  components  and,  therefore,  has  the  potential 
to  be  inaccurate  for  multi-axial  stress  states.  The  most  significant  advantage  of  this  failure 
criterion  is  that  it  identifies  the  specific  mode  of  failure  within  a  ply.  Failure  in  any  principal 
direction  of  the  material  is  predicted  when  any  of  the  following  conditions  exist: 


if  CTi>0  and  if  <t,>X1T,  then  the  failure  mode  is  fiber  tension  (61a) 

if  a,<0  and  if  |  <7,|>X1C,  then  the  failure  mode  is  fiber  compression  (61b) 

if  cT2>0  and  if  <T2>X2T,  then  the  failure  mode  is  matrix  tension  (61c) 

if  CTj  <0  and  if  |  (T2  |>X2C,  then  the  failure  mode  is  matrix  compression  (6 1  d) 

if  (T3>0  and  if  <T3>X3T,  then  the  failure  mode  is  matrix  tension  (61e) 

if  CT3  <0  and  if  |  <Tj  |>X3C,  then  the  failure  mode  is  matrix  compression  (6 1  f) 

if  |a4|  >  X23,  then  the  failure  mode  is  interlaminar  shear  (61e) 

if  IG5I  >  XI 3,  then  the  failure  mode  is  interlaminar  shear  (61f) 

if  loel  >  X12,  then  the  failure  mode  is  in-plane  shear  (61g) 


In  Equations  (61a)  through  (61g),  O’,  through  o^  are  the  six  principal  ply  stresses,  XI T  is 
the  tensile  strength  in  the  1 -direction  (longitudinal),  XI C  is  the  compressive  strength  in  the  1- 
direction,  X2T  is  the  tensile  strength  in  the  2-direction  (transverse),  X2C  is  the  compressive 
strength  in  the  2-direction,  X3T  is  the  tensile  strength  in  the  3-direction,  X3C  is  the  compressive 
strength  in  the  3-direction,  X23  is  the  shear  strength  in  the  23-plane,  XI 3  is  the  shear  strength  in 
the  13-plane,  and  X12  is  the  shear  strength  in  the  12-plane. 
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2.7. 1.3  Maximum  Strain 


The  maximum  strain  failure  criterion  (Tsai,  1987;  Nahas,  1986)  predicts  that  a  material  will 
fail  when  the  strain  in  any  direction  exceeds  its  allowable  level.  This  failure  criterion  is  similar  to 
the  maximum  stress  criterion,  except  that  it  accounts  for  some  of  the  interactions  between  the 
stresses  that  are  attributable  to  the  Poisson’s  effects  in  the  material  (i.e.,  stresses  in  the  1-  and  3- 
directions  will  affect  the  strain  in  the  2-direction).  The  failure  criterion  is  applied  in  the  exact 
same  manner  as  the  maximum  stress  failure  criterion.  The  mechanical  strains  in  the  six  directions 
(ei,  £2,  £3,  £4,  £5,  and  £5)  are  compared  to  their  corresponding  maximum  strain  allowables  (Y 1 T. 
YIC,  Y2T,  Y2C,  Y3T,  Y3C,  Y23,  Y13,  and  Y12)  in  the  same  manner  as  described  for  the 
maximum  stress  criterion,  for  example. 


if  e,>0  and  if  e,>Y  1 T,  then  the  failure  mode  is  fiber  tension  (62a) 

if  £]  <0  and  if  |  e,  |> Y 1 C,  then  the  failure  mode  is  fiber  compression  (62b) 

if  Cj  >0  and  if  >Y2T,  then  the  failure  mode  is  matrix  tension  (62c) 

if  £2  <0  and  if  |  |>  Y2C,  then  the  failure  mode  is  matrix  compression  (62d) 

if  £3>0  and  if  £3>Y3T,  then  the  failure  mode  is  matrix  tension  (62e) 

if  £3<0  and  if  |  £3|>Y3C,  then  the  failure  mode  is  matrix  compression  (62f) 
if  |£4|  >  ‘¥23,  then  the  failure  mode  is  interlaminar  shear  (62e) 

if  IesI  >  Y13,  then  the  failure  mode  is  interlaminar  shear  (62f) 

if  |£6|  >  Y12,  then  the  failure  mode  is  in-plane  shear  (62g) 


2.7. 1.4  Hydrostatic  Pressure  Adjusted 

The  hydrostatic  pressure-adjusted  failiue  criterion  is  actually  a  modified  version  of  the 
maximum  stress  failure  criterion  that  accounts  for  stress  interactions  under  compressive  loading. 
The  criterion  is  based  on  empirical  observations  reported  in  the  literature  (Hahn  and  Kallas, 
1992).  The  failure  criterion  essentially  states  that  all  the  compression  strength  allowables  in  the 
principal  directions  (XIC,  X2C,  and  X3C)  and  the  shear  strength  allowables  (X23,  X13,  and 
XI 2)  increase  with  the  hydrostatic  pressure  state  (HP)  that  exists  vvithin  the  23  planes  of  a  ply. 
A  review  of  this  phenomenon  is  discussed  in  a  separate  report  (Hoppel,  Bogetti,  Gillespie  1995). 

For  the  particular  version  of  this  criterion  presented  in  this  report,  several  assumptions  are 
made.  First,  the  “hydrostatic  pressure  state”  in  the  23  planes  referred  to  here  (HP)  is  not  that 
which  is  defined  in  the  traditional  sense.  It  is  assumed  to  be  equal  to  either  (1)  the  minimum  of 
IC2I  and  loal  (when  both  are  compressive),  (2)  the  average  of  02  and  03  (when  this  average  is 
compressive),  or  (3)  zero  when  either  of  these  two  conditions  is  not  satisfied.  In  addition,  tensile 
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strength  allowables  are  assumed  to  be  independent  of  this  hydrostatic  pressure  influence,  and  the 
strength  allowables  in  the  2-  and  3-directions  are  assumed  to  be  the  equal  (i.e.,  X2T=X3T, 
X2C=X3C),  It  is  further  assumed  that  all  shear  strength  allowables  are  equal  (i.e., 
S23=S13=S12). 

Material  strength  here  is  assumed  to  be  a  bi-linear  function  with  the  respect  to  the 
hydrostatic  pressure  state  in  the  23  planes  of  a  ply.  This  relationship  is  expressed  explicitly  as  a 
function  of  the  assumed  hydrostatic  pressure  state,  HP,  according  to 

X1C(HP)=X1C(0)+ML1«HP  for  (HP<LTP)  (63a) 

and 

X1C(HP)=(X1C(0)+ML1*LTP)+ML2-HP  for  (HP>LTP)  (63b) 

in  which  XIC(O)  is  the  usual  longitudinal  compression  strength  allowable  (i.e.,  X1C(0)=X1C); 
MLl  and  ML2  are  the  two  slopes  describing  the  bi-linear  relationship,  and  LTP  is  the 
hydrostatic  pressure  state  in  which  a  change  in  slope  of  the  XIC  versus  HP  relationship  occurs. 
Similarly,  for  the  transverse  directions, 

X2C(HP)=X3C(HP)=X2C(0)+MT1«HP  for  (HP<TTP)  (63c) 

and 

X2C(HP)=X3C(HP)=(X2C(0)+MT1*TTP)+MT2*HP  for  (HP>TTP)  (63d) 

in  which  X2C(0)  is  the  usual  transverse  compression  strength  allowable  (i.e.,  X2C(0)=X2C= 
X3C);  MTl  and  MT2  are  the  two  slopes  describing  the  bi-linear  relationship,  and  TTP  is  the 
hydrostatic  pressure  state  in  which  a  change  in  slope  of  the  X2C  (or  X3C)  versus  HP 
relationship  occurs.  For  the  shear  directions,  the  following  relationships  are  used: 

X23(HP)=X13(HP)=X12(HP)=S+MS1*HP  for  (HP<STP)  (63e) 

and 

X23(HP)=X13(HP)=X12(HP)=(S+MS1«STP)+MS2»HP  for  (HP>STP)  (63f) 

in  which  S  is  the  usual  shear  strength  allowable  (i.e.,  S=X23=X13=X12);  MSI  and  MS2  are  the 
two  slopes  describing  the  bi-linear  relationship,  and  STP  is  the  hydrostatic  pressure  state  in 
which  a  change  in  slope  of  the  S  (or  S23  or  SI  3  or  S12)  versus  HP  relationship  occurs. 
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2.7. 1.5  Tsai-Wu  Quadratic  Interaction 


The  Tsai-Wu  quadratic  interaction  or  tensor  polynomial  failure  criterion  (Tsai  and  Wu 
1971)  accounts  for  the  interactive  effects  of  a  multi -axial  stress  state.  Failure  is  predicted  when 
the  following  condition  occurs: 

F,cr,  +  +  a,)  +  +  F^^(al  +  cr^,)  +  2F,,ct^(<7,  +  a,)  + 

2F^,a^G,  -H  2F +  F6,((T^  +  >  1  (64) 


in  which  g\  through  05  are  the  principal  stresses  in  the  lamina.  The  constants  Fi,  F2,  Fi  1,  F22, 
Fi25  F23,  F44,  and  F66  are  defined  by  the  following  expressions: 


1  1 

XIT  '  XIC 
1  1 

X2T  ’  X2C 
1 

(X1T)(X1C) 

1 


(X2T)(X2C) 

2(^22  -  F33) 

1  1 


(65a) 

(65b) 

(65c) 

(65d) 

(65e) 

(65f) 


The  constants  Fi2  and  F23  are  determined  experimentally.  Methods  to  determine  these  constants 
are  described  in  Tsai  and  Wu  (1971)  and  Jiang  and  Teimyson  (1989).  This  theory  assumes  that 
the  material  is  transversely  isotropic  in  the  principal  1-2  plane  of  the  composite.  The  major 
drawback  of  this  failure  criterion  is  that  it  does  not  distinguish  among  the  various  potential  modes 
of  failure. 


2.7. 1.6  Christensen’s  Criterion 

Christensen  (1988)  proposed  a  strain-based  failure  criterion,  which  identifies  failure  as 
being  either  fiber  dominated  or  matrix  dominated  while  considering  the  multi-axial  stress  state  for 
matrix-dominated  failure.  This  criterion  identifies  three  distinct  failure  modes  for  a  composite 
lamina:  fiber  tension,  fiber  compression,  and  matrix  failure.  Christensen’s  failure  criterion  has 
been  translated  into  a  stress-based  failure  criterion  by  Hahn  and  Kallas  (1992),  and  the  stress- 
based  failure  criterion  is  employed  in  this  work. 
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Fiber  tension  and  fiber  compression,  respectively,  are  predicted  to  occur  when  either 
Equation  66a  or  66b  is  satisfied. 


(Y1T)(E,) 

O]  -  V|2(y2  ■  ^13^3 

(YICXEQ 

Oi  -  V12G2  -  V13G3 


<  1  (if  (T,  -  VijCTj  -  Vi3CT3>0) 

<  1  (if  <T,  -  VijtTj  -  V,3(J3<0) 


(66a) 

(66b) 


in  which  YIT  is  the  tensile  failure  strain  for  the  lamina  in  the  1  direction,  YIC  is  the  compressive 
failure  strain  for  the  lamina  in  the  1  direction,  and  Ei  is  the  elastic  Young’s  modulus  in  the  fiber 
direction.  Matrix-dominated  failure  is  predicted  to  occur  when  the  following  condition  exists: 

Aa,  +  B(a2  +  Cj)  +  Ccrf  +  D(cr2  +  cTj)  +  E(7,((T2  +  <73)  + 

+  Ga,^  +  H(<T5  +  (Tg)  <  1  (67) 

in  which  the  coefficients  A,  B,  C,  D,  E,  F,  G,  and  H  are  given  by 


A  = 

a(l-2v,2) 

K^Ei 

(68a) 

B  = 

a(l-V2i-V23) 

(68b) 

C  = 

2(1 +  v„)^ 

SK^Ef 

(68c) 

D  = 

3K2e2  ^21  ^21 )  (!■  ^21)^23  ■^^23] 

(68d) 

E  = 

3K^E,E,  K-1-2V21  +  V23)(1+V,2)] 

(68e) 

F  = 

3K^E^  ^^■^■*■^^21  +  2v2]  )  —  2(2  +  V21  )V23  “  V23] 

(68f) 

G  = 

1 

2K^G2^3 

(68g) 

H  = 

1 

2K^Gi^2 

(68h) 

The  constants  K  and  a  are  experimentally  determined  material  parameters,  Ei  and  E2  are  the 
elastic  moduli  in  the  principal  and  transverse  directions,  G12  and  G23  are  the  in-plane  and  out- 
of-plane  shear  moduli,  respectively,  and  the  vy’s  are  the  usual  Poisson’s  ratios  for  the  lamina. 
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2.7. 1.7  Feng’s  Failure  Criterion 

Feng  (1991)  also  proposed  a  strain-based  failure  criterion  that  differentiates  between  fiber- 
dominated  and  matrix-dominated  failure  under  multi-axial  loading.  This  criterion  determines 
failure,  based  on  the  strain  invariants  in  the  lamina.  Matrix-dominated  failure  is  predicted  to 
occur  when  the  following  relation  exists: 

AiJi  +  AiiJi^  +  A2J2  -  1  ^  0  (69) 

in  which  Ai»  Aj  1,  and  A2  are  empirically  determined  parameters  and  Ji  and  J2  are  the  strain 
invariants  given  by 

Ji=ei+e2  +  e3  (70a) 

h  =  {[(ei  -  £2)2  +  (£3  -  £2)2  +  (£1  -  £3)2]  /  6}  £42  -h  £52  +  £^2  (70b) 

in  which  £1,  £2,  £3^  £4^  £5^  and  £5  are  the  ply  strains.  Fiber-dominated  failure  is  predicted  to  occur 
when  the  following  is  satisfied: 

A5J5  -I-  A55J52  +  A4J4  -  1  >  0  (71) 

in  which  As*  A55,  and  A4  are  experimentally  determined  parameters,  and  J4  and  J5  are  the  strain 
invariants  given  by 

J4  =  £4^  +  £5^ 

J5  =  £i 

2.7. 1.8  Modified  Hashin’s  Failure  Criterion 

Hashin  (1980)  proposed  a  stress-based  failure  criterion  for  composite  materials,  which 
considers  the  tri-axial  stress  state  for  matrix  failure  modes  but  only  considers  the  uni-axial  stress 
state  in  the  fiber  direction  for  fiber-dominated  failure  modes.  In  the  fiber  direction,  this  criterion 
is  the  same  as  the  maximum  stress  failure  criterion.  Gipple,  Nuismer,  and  Camponeschi  (1995) 
suggested  modifications  of  the  failure  criterion  proposed  by  Hashin  in  which  they  assume  that 
compressive  matrix  failure  occurs  because  of  a  shearing  mechanism  rather  than  through 
compression.  The  modified  Hashin’s  failure  criterion  is  currently  being  used  by  several 
government  agencies  and  contractors  and  is  presented  in  this  report. 

For  the  fiber-dominated  failure  modes,  the  modified  Hashin  failure  criterion  distinguishes 
between  tension  and  compression  according  to  the  following  conditions: 


(72a) 

(72b) 
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if  <7,>0  and  if  cr,  >X1T,  then  the  failure  mode  is  fiber  tension 
if  (T,<0  and  if  |  <t,|>X1C,  then  the  failure  mode  is  fiber  compression 


(73a) 

(73b) 


Matrix-dominated  failure  is  predicted  when  either  of  the  following  two  conditions  exists: 


[XITJ  1X2 


X2C 


>  1  (for  >0  (matrix  tension))  (74a) 

uy  VAizy 

^  >  1  (for  <0  (matrix  compression))  (74b) 


XI2; 


in  which  the  normal  (Onn),  normal-transverse  (Ont),  and  normal-longitudinal  (Oni)  stresses  are 
evaluated  according  to  the  following: 
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Here,  the  angle  P  defines  an  orientation  in  the  2-3  plane  in  which  the  maximum  matrix  stress  state 
in  the  lamina  exists. 


2.7.2  First  Ply  and  Progressive  Ply  Failure  Theories 

While  the  failure  assessment  of  a  single  ply  or  lamina  is  clearly  defined  by  the  particular 
failure  criterion  used,  the  strength  or  failure  assessment  of  the  overall  laminate  is  somewhat  less 
straightforward.  For  example,  for  a  given  set  of  mechanical  and/or  thermal  loading  imposed  on 
the  laminate,  any  number  of  various  failure  modes  could  occur  in  one  or  more  plies.  A  single 
failure  or  multiple  failures  on  the  ply  level  do  not  necessarily  result  in  catastrophic  failure  of  the 
laminate.  In  fact,  depending  on  the  laminate  stacking  sequence  and  loading,  a  laminate  will  often 
undergo  several  minor  ply-level  failures  before  catastrophic  failure  in  the  laminate  occurs. 


The  term  first  ply  failure  refers  to  an  approach  wherein  the  overall  laminate  strength  is 
determined  by  the  lowest  point  of  laminate  load  that  causes  any  failure  mode  to  occur  in  any  ply 
vsdthin  the  laminate.  Although  this  approach  is  not  uncommon,  it  generally  provides  laminate 
strength  predictions  that  are  far  too  conservative. 
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A  more  practical,  and  often  more  accurate  approach  is  referred  to  as  a  progressive  ply 
failure  theory.  In  this  situation,  the  laminate  load  is  increased  to  a  point  when  failure  is  first 
predicted  in  any  mode  within  a  ply  within  the  laminate  (first  ply  failure).  This  load  level  is 
noted.  The  corresponding  lamina  stiffness  or  stiffnesses  for  that  ply  are  reduced  to  an 
insignificant  value  and  the  laminate  is  re-loaded  until  the  next  failure  is  detected.  This  allows  ply- 
level  stresses  to  redistribute  within  the  laminate  and  simultaneously  prevents  stress  levels  from 
increasing  in  the  previously  detected  failure  modes  of  failure  plies.  The  procedure  of  loading  to 
failure  and  reducing  corresponding  stif&iesses  is  continued  in  an  iterative  manner  until  the 
laminate  can  no  longer  support  the  intended  applied  mechanical  loading.  Iterations  are 
discontinued  when  laminate  strains  reach  some  predefined,  arbitrarily  large  value  (representing 
laminate  failure).  The  ultimate  laminate  strength  is  defined  as  the  largest  load  level  reached  during 
the  loading  strategy.  This  load  level  does  not  necessarily  have  to  be  the  last  load  level  reached 
during  the  solution  strategy.  In  fact,  load  levels  rise  and  fall  during  the  iterative  solution  strategy 
as  they  are  dictated  by  the  complex  interactions  of  load  transfer  between  the  different  modes  of 
all  the  plies  within  the  laminate.  The  particular  mode  of  failure  and  the  actual  ply  that 
corresponds  to  the  largest  load  level  reached  during  the  solution  strategy  are  identified  as  the 
critical  mode  and  critical  ply  responsible  for  ultimate  laminate  failure. 

2.7.3  Failure  Assessment  Summary 

Laminate  strength  predictions  can  be  based  on  first  ply  or  progressive  ply  failure  analyses. 
For  either  approach,  an  ultimate  load  level  can  be  predicted  as  a  measure  of  the  quantitative 
performance  of  the  laminate.  It  is  often  desirable  to  describe  the  ultimate  laminate  load  or 
laminate  strength  in  terms  of  a  safety  factor  (SF).  Such  a  safety  factor  is  defined  as  the  ultimate 
load  that  a  laminate  could  sustain  (as  defined  by  first  ply  or  progressive  ply  failure)  divided  by 
the  applied  loading  on  the  laminate.  This  ultimate  load  would,  of  course,  be  some  scalar  multiple 
of  the  applied  three-dimensional  loading.  According  to  this  definition,  safety  factors  less  than 
one  would  represent  laminate  failures. 

A  qualitative  description  of  failure  is  also  very  useful  in  understanding  laminate  strength. 
Corresponding  to  the  ultimate  load  predicted  for  the  laminate,  a  criticeil  mode  of  failure  and  a 
critical  ply  can  be  defined.  This  information  is  important  since  it  facilitates  the  redesign  or 
optimization  of  laminate  architectures  for  particular  loading  requirements.  The  LAM3D  code 
generates  laminate  strength  predictions  in  terms  of  safety  factor  and  provides  critical  failure  mode 
and  critical  ply  information. 
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2.8  T.AM3D:  A  Three-Dimensional  Laminated  Composite  Media  Analysis  Software  Program 

2.8.1  Program  Options 

The  theory  that  has  been  described  in  this  report  has  been  codified  into  a  user-friendly 
FORTRAN-based  computer  program  called  LAM3D.  This  section  describes,  in  detail,  the  input 
and  output  format  specifically  for  version  3.0  of  LAM3D.  The  code  is  intended  to  be  used  for 
three-dimensional  laminated  media  analysis  and  calculates  (1)  effective  properties,  (2)  ultimate 
strengths  under  mechanical  and/or  thermal  loading,  and  (3)  the  effective  stress  versus  strain 
response  of  the  laminate  because  of  mechanical  loading  under  progressive  ply  failure.  The 
laminate  may  be  composed  of  different  types  of  ply  materials  (e.g.,  hybrid),  and  a  non-uniform 
thermal  loading  may  be  superimposed  on  the  laminate  before  mechanical  loading.  The 
progressive  ply  failure  analysis  is  incremental  in  loading.  At  every  point  where  ply-level  failure 
is  detected  during  the  loading  history,  any  or  all  of  the  effective  laminate  or  ply-level  details 
(properties,  stresses,  and  strains)  are  calculated  and  are  available  for  examination  (output). 

Failure  analysis  information  available  for  output  includes  the  safety  factor,  the  critical  ply,  and 
the  associated  failure  mode. 

Any  self-consistent  set  of  units  may  be  used  in  the  LAM3D  (version  3.0)  analysis. 

Typically  for  SI  imits,  pascals  (Pa)  are  used  for  pressure,  meters  (m)  for  distance,  kilograms  per 
meter  cubed  (kg/m^)  for  density,  and  degrees  Celsius  (°C)  for  temperature.  For  English  units, 
pounds  per  square  inch  (Ib/in^)  are  used  for  pressure,  inches  (in.)  for  distance,  slugs  per  inch 
cubed  (slug/in^)  for  density,  and  degrees  Fahrenheit  (°F)  for  temperature. 

The  program  LAM3D  (version  3.0)  is  executed  in  an  interactive  mode  and  uses  a 
convenient  data  base  format  for  conducting  parametric  studies.  Details  of  the  LAM3D  (version 
3.0)  program  execution  are  described  in  a  subsequent  section. 

2.8.2  Input  Summary 

2.8.2. 1  Material  DataBase 

The  LAM3D  (version  3.0)  data  base  input  file  is  used  to  define  the  laminate  to  be  analyzed. 
The  sample  data  base  presented  in  Appendix  A  will  be  used  to  illustrate  the  data  base  file  format. 

In  the  data  base,  the  data  are  separated  by  commas  and  spaces  on  each  line.  The  exclamation 
point  (!)  represents  a  comment  statement;  all  the  information  to  the  left  of  the  exclamation  point 
is  data  which  are  read  directly  as  input  into  the  program;  the  information  to  the  right  is  ignored 
by  the  program.  The  data  base  is  subdivided  into  three  sections:  the  first  consists  of  three 
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header  cards,  the  second  describes  laminate  details,  and  the  third  is  essentially  a  “material 
library.” 


2.8.2. 1.1  Section  1:  Header  Cards 

/ 

The  first  and  second  header  cards  are  set  to  the  number  1.  The  third  header  card 
defines  the  number  of  different  types  of  materials  that  are  listed  in  the  data  base  file.  (In  this 
example,  this  card  is  set  to  3,  since  three  materials  are  listed:  aluminum,  IM7-graphite/8551- 
epoxy,  and  S2-glass/3501-epoxy.) 

2.8.2. 1.2  Section  2:  Laminate  Architecture 

The  second  section  of  the  data  base  defines  several  important  parameters  for  the 
laminate  analysis.  The  first  card  of  this  section  contains  a  region  identification  number 
(REG_ID),  the  failure  criterion  identification  number  (FAIL_CRT),  the  maximum  number  of 
iterations  for  progressive  failure  (ITERS),  and  a  text  parameter  flag  (ITEXT).  The  region 
identification  number  (REG_ID)  is  always  set  to  the  number  1.  The  failure  criterion 
identification  number  (FAIL_CRT)  defines  the  failure  criterion  that  will  be  used  for  all  the  plies 
in  the  laminate.  A  value  between  1.0  and  8.0  is  required  with  (1.0)  corresponding  to  the  Von 
Mises-Henecky  failure  criterion,  (2.0)  to  the  maximum  stress  failure  criterion,  (3.0)  to  the 
maximum  strain  failure  criterion,  (4.0)  to  the  hydrostatic  pressure-adjusted  failure  criterion,  (5.0) 
to  the  Tsai-Wu  quadratic  interaction  failure  criterion,  (6.0)  to  Christensen’s  failure  criterion,  (7.0) 
to  Feng’s  failure  criterion,  and  (8.0)  to  the  modified  Hashin  failure  criterion.  The  number  of 
iterations  (ITERS)  is  used  to  define  the  maximum  number  of  times  the  program  will  allow  the 
laminate  to  be  loaded  during  the  progressive  ply  failure  analysis  (a  value  of  five  is  usually 
sufficient).  The  text  parameter  flag  (ITEXT)  is  reserved  to  accommodate  future  enhancements  of 
the  LAM3D  (version  3.0)  code.  It  should,  however,  be  set  equal  to  zero. 

The  second  card  in  this  section  contains  three  orientation  parameters:  (BETA), 

(PHI),  and  (SI).  These  parameters  are  reserved  to  accommodate  future  enhancements  of  the 
LAM3D  (version  3.0)  code  and  should  all  be  set  equal  to  zero.  All  spaces  designated  as  ( — ) 
throughout  the  data  base  file  are  also  reserved  to  accommodate  future  enhancements  of  the 
LAM3D  (version  3.0)  code  and  should  all  be  set  equal  to  zero.  Consequently,  the  fourth 
parameter  on  this  line  should  be  set  equal  to  zero. 

The  third  card  defines  the  number  of  plies  in  the  laminate  (NPLY)  and  the  region  type 
of  the  laminate  (REG_TYP).  The  parameter  (NPLY)  defines  the  number  of  input  cards  (lines) 
that  LAM3D  (version  3.0)  expects  to  read  before  it  reads  input  from  the  material  library  section. 
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The  region  type  parameter  (REG_TYP)  is  reserved  to  accommodate  future  enhancements  and 
should  always  be  set  to  1 .0.  Note  that  the  total  number  of  cards  in  this  section  will  be  equal  to 
(NPLY+3). 


The  following  (NPLY)  number  of  cards  are  used  to  define  the  individual  ply  materials 
(MATJD),  thicknesses  (THICKNESS),  orientation  angles  in  degrees  (ANGLE)  and  thermal  ply 
loading  (TEMP).  The  material  identification  number  (MATJD)  specifies  the  ply  material  type 
(defined  in  the  material  library  section  of  the  database)  for  each  ply.  (Note  that  by  specifying  the 
material  type  for  each  ply  individually,  the  analysis  of  hybrid  laminates  is  easily  accommodated.) 
The  thickness  parameter  (THICKNESS)  is  the  relative  thickness  of  the  ply  in  the  laminate.  The 
ply  angle  (ANGLE)  is  the  orientation  between  the  longitudinal  (1 -direction)  axis  of  the  ply  and 
the  longitudinal  axis  (x-direction)  of  the  laminate  (see  Figure  3).  The  thermal  ply  loading  (TEMP) 
defines  the  difference  between  the  operating  temperature  and  processing  temperature  of  the  ply. 
(Note  that  by  specifying  the  thermal  ply  loading  in  this  manner,  the  laminate  is  not  restricted  to  a 
uniform  thermal  loading  profile.) 

2.8.2. 1 .3  Section  3 :  Material  Library 

The  last  section  of  the  data  base  file  is  structured  in  what  could  be  considered  a 
“material  library.”  This  section  contains  the  three-dimensional  ply  properties,  densities,  and 
strength  parameters  for  all  the  materials  included  in  the  data  base.  The  number  of  materials  listed 
in  this  section  should  be  equal  to  that  defined  in  the  third  card  of  the  first  section.  Each 
individual  material  listed  in  the  data  base  file  occupies  45  cards  (lines),  with  three  columns 
(entries)  for  each  card.  The  45-card  format  for  one  material  is  described  next.  (Additional 
materials  included  in  the  data  base  file  would  be  inserted  in  exactly  the  same  format.) 

The  first  card  contains  the  material  identification  number  (MATJD),  the  material 
type  (MAT_TYPE),  and  the  material  density  (MAT_DENS).  The  first  material  listed  in  the  file 
would  have  an  identification  number  (MATJD)  of  1.0,  the  second  2.0,  and  so  on.  The  material 
type  parameter  (MAT  TYPE)  is  reserved  to  accommodate  future  enhancements  and  should 
always  be  set  to  1.0.  The  material  density  (MAT_DENS)  should  also  be  specified  for  the 
material. 


The  second  card  contains  the  elastic  moduli  for  the  material  defined  in  its  three 
principal  coordinate  directions  (i.e.,  Ei,  E2,  E3).  (Note  that  if  the  material  is  isotropic,  the  same 
value  should  be  entered  for  all  three  moduli.)  The  third  card  contains  the  three  major  Poisson’s 
ratios  for  the  material  (i.e.,  V23,  V13,  V12).  The  fourth  card  contains  the  three  shear  moduli  for  the 
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material  (i.e.,  623,  G13,  G12).  The  fifth  card  contains  the  thermal  expansion  coefficients  in  the 
three  principal  coordinate  directions  (i.e.,  an,  a22,  OC33). 

Cards  6  through  45  contain  all  the  parameters  required  to  define  the  failure  criteria  of 
the  material.  For  each  failure  criteria  (1  through  8),  the  data  base  reserves  five  cards  (lines)  each 
Avith  three  columns.  Since  there  are  eight  failure  criteria,  40  cards  are  reserved  for  each  material  to 
define  all  the  potentially  required  failure  parameters. 

The  Von  Mises  failure  criterion  is  the  first  one  listed.  The  only  datum  the  program 
requires  for  this  failure  criterion  is  the  material  yield  strength,  (Oy),  which  should  be  listed  in  the 
first  column  of  the  first  card  (designated  as  VMl  in  the  sample  data  base  file  in  Appendix  A). 
Since  the  program  does  not  require  any  additional  information  for  this  failure  criterion,  the 
remaining  14  positions  should  all  be  defined  as  0.0.  (Again,  see  the  aluminum  material  in  the  data 
base  file  in  Appendix  A.) 

The  second  failure  criterion  listed  in  the  data  base  file  is  maximum  stress,  which 
requires  nine  constants.  The  nine  constants  (XIT,  XIC,  X2T,  X2C,  X3T,  X3C,  X23,  X13,  and 
X12,  defined  in  Section  2.7,1.2)  are  input  in  the  format  indicated  by  the  IM7-graphite/8551- 
epoxy  material  (MAT_ID=2.0)  entry  listing  in  the  example  data  base  file. 

The  third  failure  criterion  is  maximum  strain  failure,  which  also  requires  nine 
constants.  The  constants  (YIT,  YIC,  Y2T,  Y2C,  Y3T,  Y3C,  Y23,  Y13,  and  Y12,  defined  in 
Section  2.7 . 1 .3)  are  input  in  the  format  indicated  by  the  IM7-graphite/855 1-epoxy  material 
(MAT_ID=2.0)  entry  listing  in  the  example  data  base  file. 

The  fourth  failure  criterion  is  hydrostatic  pressure-adjusted  failure,  which  requires  15 
constants.  The  constants  required  are  defined  in  Section  2.7. 1.4  and  their  format  is  indicated  by 
the  IM7-graphite/8551 -epoxy  material  (MAT_ID=2.0)  entry  listing  in  the  example  data  base  file. 
The  variable  HPD  is  used  to  specify  the  manner  in  which  the  hydrostatic  pressure  state  is 
calculated.  For  HPD=1.0,  the  hydrostatic  pressure  state  is  assumed  to  be  equal  to  the  minimum 
of  |a2|  and  |a3|  (when  they  are  both  compressive)  and  zero  if  either  02  or  03  are  tensile.  For 
HPD=2.0,  the  hydrostatic  pressure  state  is  assumed  to  be  equal  to  the  average  of  ^2  and  03 
(when  this  average  is  compressive)  and  zero  if  this  average  is  tensile. 

The  fifth  failure  criterion  is  the  Tsai-Wu  quadratic  interaction,  which  requires  seven 
parameters  to  be  input  into  the  data  base.  The  constants  required  are  defined  in  Section  2.7. 1 .5 
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and  their  format  is  indicated  by  the  IM7-graphite/8551 -epoxy  material  (MAT_ID-2.0)  entry 
listing  in  the  example  data  base  file. 

Christensen’s  failure  criterion  is  the  sixth  listed,  and  it  requires  13  parameter  entries  in  the 
data  base  file.  The  parameters  El,  E2,  E3,  V23,  V13,  V12,  G23,  G13,  and  G12  are  the  usual  ply 
properties.  The  parameters  K,  a  (ALPHA),  YIT  and  YIC  are  defined  in  Section  2. 7.1. 6  and 
their  format  is  indicated  by  the  IM7-graphite/8551 -epoxy  material  (MAT_ID=2.0)  entry  listing 
in  the  example  data  base  file. 

The  seventh  failure  criterion  is  Feng’s,  which  requires  the  following  six  parameters:  Al, 
A1 1,  A2,  A4,  A5,  and  A55  (see  Section  2.7.1.7).  Again,  the  format  for  these  parameters  is 
indicated  by  the  IM7-graphite/8551 -epoxy  material  (MAT_ID=2.0)  entry  listing  in  the  example 
data  base  file. 


The  eighth  failure  criterion,  modified  Hashin’s  failure,  requires  the  same  nine  parameters  as 
the  maximum  stress  failure  criterion  (XIT,  X2T,  X3T,  XIC,  X2C,  X3C,  X23,  X13,  and  X12). 
The  same  format  for  entry  into  the  database  file  as  the  maximum  stress  criterion  is  also  used. 
(See  the  lM7-graphite/8551 -epoxy  material  (MAT_ID=2.0)  entry  listing  in  the  example  data 
base.) 

2.8.2.2  LAM3D  (version  3.0)  Program  Execution 

Execution  of  LAM3D  is  done  in  an  interactive  mode.  The  following  section  details  the 
interactive  questions  that  are  encountered  in  a  typical  program  run. 


After  appropriate  compilation  and  execution  commands  are  performed,  the  first  required 
user  input  for  LAM3D  is  entered  to  define  the  desired  output  to  be  written  (stored)  to  an 
output  file.  The  following  text  will  be  printed. 


ENTER  (1)  YES  OR  (0)  NO  FOR  THE  FOLLOWING  OUTPUT 
IN  SEQUENCE  ON  ONE  LINE  OF  INPUT  . . . 

(1)  OUTPUT  OF  INPUT 

(2)  CIJ  MATRICES 

(3)  EFFECTIVE  PROPERTIES 

(4)  LOADING 

(5)  MECHANICAL  PLY  STRESS  AND  STRAIN 

(6)  THERMAL  PLY  STRESS  AND  STRAIN 

(7)  FAILURE  ANALYSIS  DETAILS 
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For  each  of  the  seven  items  on  the  list,  the  user  should  enter  a  one  (1)  if  he  or  she  wants  those 
data  to  be  included  in  the  program  output  or  a  zero  (0)  if  those  data  are  not  to  be  included  in  the 
output.  For  example,  to  display  all  the  information  offered  in  the  output,  type 

1, 1, 1, 1, 1, 1, 1  [CR] 

in  which  [CR]  is  the  carriage  return.  LAM3D  will  then  ask  the  user  what  type  of  lamina 
constitutive  model  should  be  assumed  (see  Section  2.2  for  model  details): 

ENTER  TYPE  OF  LAMINA  CONSTITUTIVE  MODEL  > 

(1)  ORTHOTROPIC  (9  I.E.C.) 

(2)  TRANSVERSELY  ISOTROPIC  (5  I.E.C.) 

(3)  IN-PLANE  MODEL  (3  I.E.C.) 

In  this  example  problem,  the  orthotropic  lamina  constitutive  model  will  be  used,  so  the  user 
should  respond  by  typing 

1[CR] 

LAM3D  will  then  ask  the  user  what  type  of  laminate  constitutive  model  will  be  used  (see  Enie 
and  Rizzo  [1970]  for  simplified  model  and  Section  2.4  or  Chou,  Carleone,  and  Hsu  [1972]  for 
exact  model  details): 

ENTER  TYPE  OF  LAMINATE  CONSTITUTIVE  MODEL  > 

(1)  SIMPLIFIED  (RIZZO,  1970) 

(2)  EXACT  (CHOU,  1972) 

In  this  example  problem,  the  exact  laminate  constitutive  model  will  be  employed,  so  the  user 
should  respond  by  typing 

2[CR] 

LAM3D  will  then  ask  the  user  to  enter  the  mechanical  loading  (three-dimensional  stresses)  to  be 
applied  on  the  laminate  in  the  x,  y,  and  z  directions  in  units  of  pounds  per  square  inch: 

ENTER  AVERAGE  STRESS  X,  Y,  Z  (lb/in2)  > 

The  following  response  would  indicate  that  a  load  of  -1000  Ib/in^  was  applied  in  the  x-direction: 
-1000,  0, 0  [CR] 
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LAM3D  will  then  ask  the  user  to  enter  the  applied  shear  stresses  on  the  laminate  in  the  yz-,  xz-, 
and  xy-directions  in  units  of  pounds  per  square  inch: 

ENTER  AVERAGE  SHEAR  STRESSES  YZ,  XZ,  XY  (lb/in2)  > 

The  following  response  would  indicate  an  applied  load  of  10  Ib/in^  in  the  yz-direction,  35  Ib/in^ 
in  the  xz-direction,  and  10  Ib/in^  in  the  xy-direction: 

10, 35, 10  [CR] 

(Recall  that  all  thermal  loads  are  defined  through  the  data  base  file  [see  Section  2. 8.2.1],)  With 
the  following  prompt,  the  name  of  the  pre-established  input  data  base  file  is  entered: 

ENTER  NAME  OF  DATA  BASE  FILE  > 
dbase.in  [CR] 

Finally,  LAM3D  will  ask  the  user  to  enter  the  name  of  the  file  where  all  requested  output  will  be 
stored: 

ENTER  NAME  OF  OUTPUT  FILE  > 
out  [CR] 

After  this  input  is  entered,  the  program  will  begin  to  run.  LAM3D  will  write  various  comments 
on  the  screen  as  it  is  reading  fi-om  the  material  data  base  and  while  it  is  computing.  LAM3D  will 
also  write  the  user-requested  information  during  execution  to  the  indicated  output  file.  During 
the  failure  analysis,  if  the  degradation  of  the  laminate  because  of  progressive  ply  failure  is 
detected,  the  program  will  report  the  following: 

LOCAL  LAMINATE  STRAIN  EXCEEDS  TOLERANCE  . . . 

In  addition,  if  the  initially  applied  thermal  stresses  in  the  laminate  cause  failure  (in  any  ply),  the 
program  will  report  the  following: 

THERMAL  LOADS  EXCEED  ALLOWABLE  FOR  STRESS  . . . 

PROGRAM  TERMINATED  ... 

PLEASE  RESTART  WITH  REDUCED  THERMAL  LOADING  . . . 

The  output  file  format  generated  by  LAM3D  is  discussed  in  the  next  section. 
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2.8.3  Output  File  Format  Summary 

A  sample  output  file  produced  by  LAM3D  (generated  by  the  input  described  in  the 
previous  sections  of  this  report)  is  shown  in  Appendix  B.  The  output  file  can  be  used  to  store  a 
rather  extensive  amount  of  information  related  to  the  laminate  analysis.  For  illustrative  purposes, 
all  the  output  options  have  been  selected  for  this  example.  However,  since  all  this  information  is 
generally  not  required,  the  output  requests  will  generally  be  significantly  reduced.  First,  output 
written  to  the  file  is  a  summary  of  the  input  that  was  used  to  define  the  laminate  analysis.  This 
includes  information  that  is  both  user  input  and  data  base  input  (e.g.,  type  of  analysis  options, 
ply  properties,  failure  criteria,  loading,  etc.). 

As  the  failure  analysis  is  being  executed,  the  program  will  report  the  requested  detailed 
information  about  the  laminate  that  is  “changing”  because  of  the  progressive  ply  failures.  For 
example,  after  each  iteration  in  the  failure  analysis,  the  following  (if  all  are  requested)  are  written 
as  output:  the  Cjj  ply  matrices  in  the  principal  ply  coordinate  system,  the  Cy  ply  matrices  in  the 
global  laminate  coordinate  system,  the  effective  Cy  stif&iess  matrix  for  the  laminate,  the  effective 
engineering  constants  for  the  laminate,  the  applied  mechanical  laminate  stresses,  the  applied 
thermal  ply  loading,  the  effective  mechanical  laminate  strains,  the  effective  thermal  laminate 
strains,  the  mechanical  and  thermal  ply  stresses  in  both  the  principal  ply  and  global  laminate 
coordinate  systems,  and  a  detailed  description  of  the  failure  assessment  for  that  iteration  in  the 
failure  analysis.  The  failure  assessment  at  each  iteration  gives  a  cumulative  summary  of  the  ply 
failure  modes  that  have  occurred  as  well  as  the  associated  safety  factor  for  the  iteration. 

3.  RESULTS  AND  DISCUSSION 

3.1.  Lamina  Properties  and  Strength  Allowables 

In  this  section,  LAM3D  predictions  for  the  effective  thermo-mechanical  laminate 
properties  are  validated  against  results  presented  in  the  published  literature.  Laminate  properties 
and  strength  predictions  are  also  correlated  against  experimental  data.  In  these  correlations,  three 
different  types  of  unidirectional  composite  material  systems  are  considered:  AS4/3501, 
IM7/PEEK,  and  IM9/8551.  The  three-dimensional  unidirectional  lamina  properties  and  strength 
allowables  used  in  the  LAM3D  predictions  are  sximmarized  in  Tables  1  and  2,  respectively.  The 
AS4/3501  lamina  properties  were  taken  from  the  published  literature  (U.S.  Army  Materials 
Technology  Laboratory  1992).  The  unidirectional  lamina  properties  and  strength  allowables  used 
for  the  IM7/PEEK  and  IM9/8551  material  systems  were  generated  under  a  separate 
characterization  effort  (Gillespie  1994). 
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Table  1 


Unidirectional  Composite  Lamina  Properties  Used  in  LAM3D  Predictions 

(E  and  G  in  units  of  MSI) 


AS4/3501 

(MIL-HDBK-17,  1992) 

IM7/PEEK 
(Gillespie  1994) 

IM9/8551 
(Gillespie  1994) 

El 

16.5 

22.19  ±2.83 

22.80 1 2.65 

E2 

1.40 

1.1910.14 

1.2310.08 

Es 

1.40 

1.1910.14 

1.2310.08 

V23 

0.54 

0.45 

0.45 

V13 

0.33 

0.33 

0.33 

Vi2 

0.33 

0.33 

0.33 

G23 

0.45 

0.41 1 0.05 

0.47 1 0.04 

Gi3 

0.87 

1.4510.09 

1.3710.11 

Gi2 

0.87 

1.4510.09 

1.3710.11 

Table  2 

Composite  Lamina  Strength  Allowables  Used  in  LAM3D  Predictions 
(allowables  in  units  of  KSI) 

IM7/PEEK 

IM9/8551 

(Gillespie  1994) 

(Gillespie  1994) 

XIT 

2611.1717.87 

415.05119.85 

X2T 

10.4311.16 

9.7610.19 

X3T 

10.4311.16 

9.76  ±0.19 

XIC 

138.1516.88 

147.26  1 12.89 

X2C 

30.5812.68 

28.9811.26 

X3C 

30.5812.68 

28.981  1.26 

S23 

10.9010.42 

24.97 1 0.23 

S13 

10.90 1 0.42 

24.97 1 0.23 

S12 

10.9010.42 

24.97 1 0.23 
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3.2  Effective  Thermo-Mechanical  Laminate  Property  Predictions 


Effective  thermo-mechanical  laminate  properties  generated  from  the  LAM3D  code  have 
been  validated  against  other  theoretical  predictions  in  the  published  literature  for  some  selected 
laminate  architectures.  Correlation  of  the  effective  laminate  mechanical  properties  are  presented . 
next.  It  is  noted  that,  although  not  presented  here,  the  effective  laminate  thermal  expansion 
coefficients  have  been  validated  against  predictions  reported  by  Pagano  (1974c). 

Theoretical  predictions  of  the  effective  mechanical  engineering  constants  for  several 
AS4/3501  composite  laminate  architectures  by  Pagano  (1974a),  Sim  and  Li  (1988),  Roy  and  Tsai 
(1992)  and  classical  laminated  plate  theory  (CLPT)  (Whitney  1987)  are  compared  with  LAM3D 
predictions  (Chou,  Carleone,  &  Hsu  1972).  In  Table  3,  the  comparison  of  predictions  for  the 
following  laminate  architectures  is  made:  [0/0/90/90/0/0],  [0/90/+45/-45/-45/+45/90/0]  and  [+30/- 
30/-30/+30].  Excellent  agreement  between  the  predictions  gives  indication  that  the  LAM3D 
predictions  are  accurate  and  consistent  with  the  proposed  theory  of  (Chou,  Carleone,  &  Hsu 
1972). 

Correlation  of  LAM3D  effective  mechanical  property  predictions  are  also  made  against  the 
experimental  results  of  Gillespie  (1994)  for  a  [0/0/-45/+45/+45/-45/0/0]  laminate  architecture  in 
both  the  IM7/PEEK  and  IM9/8551  materials.  Reasonable  agreement  between  the  measured 
effective  laminate  moduli  and  the  LAM3D  predictions  is  presented  in  Table  4  for  both  the  x-  and 
y-directions. 

3.3  Ultimate  Laminate  Strength  Predictions 

Uni-axial  strengths  of  the  [0/0/-45/+45/+45/-45/0/0]  laminate  architecture  in  both  the 
IM7/PEEK  and  IM9/8551  material  systems  were  measured  in  tension  and  compression  in  both 
the  X-  and  y-directions.  LAM3D  predictions  for  each  material  and  load  case  were  generated  for 
comparison.  The  maximum  stress  lamina  failure  criterion  was  assumed  in  all  the  predictions 
presented  in  this  correlation.  The  progressive  ply  failure  theory  discussed  in  Section  2.7.2  was 
also  implemented. 
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Table  3 


Predicted  Effective  3D  Properties  for  Various  AS4/3501-7  Laminates 
(E  and  G  in  units  of  MSI) 


CLPT 

Pagano 

Sun 

Roy 

Chou 

(Whitney  1987) 

(Pagano  1974) 

(Sun  and  Li  1988) 

(Roy  and  Tsai  1992) 

(Chou  et  al.  1972) 

[0/0/90/90/0/0]  Laminate  Properties 

Ex 

11.5 

11.5 

11.5 

11.5 

11.5 

Ey 

6.48 

6.47 

6.47 

6.47 

6.47 

Ez 

— 

1.80 

1.80 

1.65 

1.80 

Vyz 

— 

.519 

.520 

.465 

.519 

Vxz 

— 

.488 

.489 

.402 

.488 

■Vxy 

.073 

.073 

.074 

.072 

.073 

Gyz 

— 

.536 

.536 

.503 

.536 

Gxz 

— 

.664 

.664 

.780 

.664 

Gxy 

.870 

.870 

.870 

.870 

.870 

[0/90/+45/-45/-45/+45/90/0]  Laminate  Properties 

Ex 

6.68 

6.68 

6.68 

6.67 

6.68 

Ey 

6.68 

6.68 

6.68 

6.68 

6.68 

Ez 

— 

1.82 

1.82 

1.82 

1.82 

Vyz 

— 

.375 

.376 

.318 

.375 

Vxz 

— 

.375 

.376 

.317 

.375 

Vxy 

.297 

.297 

.298 

.296 

.297 

Gyz 

— 

.593 

.593 

.519 

.593 

Gxz 

— 

.593 

.593 

.627 

.593 

Gxy 

2.58 

2.57 

2.57 

2.57 

2.57 

[30/-30/-30/30]  Laminate  Properties 


Ex 

6.84 

6.84 

6.84 

6.84 

6.84 

Ey 

1.77 

1.77 

1.77 

1.77 

1.77 

Ez 

— 

1.66 

1.66 

1.50 

1.66 

Vyz 

— 

.390 

.390 

.434 

.390 

Vxz 

-.095 

-.094 

-.197 

-.095 

Vxy 

1.14 

1.14 

1.41 

1.13 

1.14 

Gyz 

— 

.512 

.512 

.515 

.512 

Gxz 

.705 

.705 

.708 

.705 

Gxy _ 

3.43 

3.43 

3.42 

3.42 

3.43 

40 


Table  4 


Correlation  of  Measured  Composite  Laminate  Properties 
with  LAM3D  Predictions 
(E  in  units  of  MSI) 


IM7/PEEK 

IM9/855] 

Experimental 
(Gillespie  1994) 

LAM3D 

(Chou  et  al.  1972) 

Experimental 
(Gillespie  1994) 

LAM3D 

(Chou  etal.  1972) 

Ex  11.9811.42 

13.47 

11.5111.21 

13.70 

Ey  2.53  1 0.39 

3.94 

2.90 1 0.22 

3.97 

Correlation  of  LAM3D  failure  analysis  predictions  with  experimental  measurements  for  the 
IM7/PEEK  laminate  is  presented  in  Table  5.  Excellent  agreement  is  noted  for  the  x-direction 
tension  and  y-direction  compression  load  cases.  Agreement  in  each  case  corresponds  to  ultimate 
longitudinal  fiber  direction  failure  (second  ply  failure  in  each  load  situation).  It  is  noted  that 
while  the  x-direction  compression  LAM3D  prediction  is  overly  conservative  (e.g.,  71.18  ksi 
predicted  versus  99.20  ksi  measured),  the  y-direction  tension  prediction  is  exceedingly  high  (e.g., 
65.29  ksi  predicted  versus  34.53  ksi  measured).  In  all  cases,  it  is  seen  that  the  first  ply  failure 
prediction  offers  a  very  conservative  laminate  failure  strength  estimate. 

Correlation  of  LAM3D  failure  analysis  predictions  with  experimental  measurements  for  the 
IM9/8551  laminate  is  presented  in  Table  6.  Reasonable  agreement  is  noted  for  the  x-direction 
tension  and  x-direction  compression  load  cases.  As  with  the  IM7/PEEK  material,  agreement  in 
each  case  corresponds  to  ultimate  longitudinal  fiber  direction  failure.  LAM3D  laminfite  strength 
predictions  are  exceedingly  high  for  the  y-direction  tension  load  case  (e.g.,  73.63  ksi  predicted 
versus  42.39  ksi  measured)  as  well  as  the  y-direction  compression  load  case  (e.g.,  66.57  ksi 
predicted  versus  39.92  ksi  measured).  It  is  noted  that  for  the  IM9/8551  material,  a  first  ply 
failure  prediction  offers  overly  conservative  laminate  strength  predictions  in  three  of  the  four  load 
cases. 

The  discrepancies  between  the  predicted  and  measured  results  are  not  surprising.  It  is  a 
well-accepted  fact  that  the  prediction  of  ultimate  laminate  strength  is  not  always  done  with  great 
accuracy.  This  is  largely  because  the  direct  translation  of  lamina  properties  into  the  laminate 
configuration  is  difficult  to  attain.  Laminate  strength  predictions,  much  more  so  than  the 
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mechanical  property  predictions,  are  subject  to  a  wide  degree  of  scatter.  Laminate  strengths  can 
be  influenced  by  many  factors  such  as  test  specimen  fabrication/preparation  and  specimen 
alignment  in  the  load  fixture  as  well  as  the  type  of  test  method  selected.  Furthermore,  the 
particular  lamina  failure  criterion  employed  as  well  as  the  assumptions  made  in  the  progressive 
ply  failure  analysis  could  have  a  significant  effect  on  the  ultimate  laminate  strength  predictions. 

4.  CONCLUSIONS 

An  analytical  model  for  predicting  the  effective  three-dimensional  thermo-mechanical 
properties  and  ultimate  strength  of  thick  laminated  composite  media  under  applied  mechanical 
and/or  thermal  loadings  was  presented.  The  theoretical  basis  for  the  model  is  based  on  the  work 
of  Chou,  Carleone,  and  Hsu  (1972).  Various  types  of  lamina  constitutive  models  were 
incorporated  in  the  model.  Numerous  lamina  failure  criteria  and  a  laminate  failure  analysis 
methodology  that  employs  a  progressive  ply  failure  strategy  were  also  developed.  Model 
predictions  were  validated  with  comparisons  between  results  in  the  published  literature  and 
available  experimental  data. 

A  major  result  of  this  work  is  the  codification  of  a  three-dimensional  laminated  media  model 
into  a  user-fnendly  computer  program  environment,  entitled  LAM3D,  for  implementation  as  an 
engineering  design  tool  for  thick  laminated  composite  structures.  Several  illustrative  examples  are 
presented  to  demonstrate  the  utility  of  LAM3D  for  predicting  the  effective  properties  and 
ultimate  strength  of  laminated  composite  media. 

This  work  represents  an  important  step  toward  the  development  of  a  computationally 
efficient  engineering  tool,  which  will  significantly  enhance  the  design  and  analysis  capability  of  a 
wide  range  thick  composite  structures  for  future  military  and  commercial  applications. 
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Table  5 


Correlation  of  Measured  Composite  Laminate  Strengths  with  LAM3D  Predictions  for  an 
IM7/PEEK  [O2/+45/-45/-45/+45/O2]  Laminate  (units  of  ksi) 


X-  Direction  Tension 

Experimentally  Measured  Strength: 

LAM3D  Failure  Analysis  Predictions: 
Failure  Iteration  Safety  Factor 

1  64.30 

2  134.57 

135.80±3.30  (Gillespie  1994) 

Ply  Mode 

±45  S12 

0  XIT 

X  -  Direction  Compression 

Experimentally  Measured  Strength: 

99.20±2.95  (Gillespie  1994) 

LAM3D  Failure  Analysis  Predictions: 

Failure  Iteration 

Safety  Factor 

Ply 

Mods 

1 

64.30 

±45 

S12 

2 

71.18 

0 

XIC 

Y  -  Direction  Tension 

Experimentally  Measured  Strength: 

34.53±8.91  (Gillespie  1994) 

LAM3D  Failure  Analysis  Predictions: 

Failure  Iteration 

Safety  Factor 

Ply 

Mode 

1 

25.35 

±45 

S12 

2 

27.32 

0 

X2T 

3 

39.81 

±45 

X2T 

4 

65.29 

±45 

XIT 

Y  -  Direction  Compression 

Experimentally  Measured  Strength: 

49.13±2.43  (Gillespie  1994) 

LAM3D  Failure  Analysis  Predictions: 

Failure  Iteration 

Safety  Factor 

Ply 

Mode 

1 

25.35 

±45 

S12 

2 

45.68 

±45 

XIC 
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Table  6 


Correlation  of  Measured  Composite  Laminate  Strengths  with  LAM3D  Predictions  for  an 
IM9/8551  [O2/+45/-45/-45/+45/O2]  Laminate  (units  of  ksi) 


X-  Direction  Tension 

Experimentally  Measured  Strength;  186.25±7.97  (Gillespie  1994) 

LAM3D  Failure  Analysis  Predictions: 

Failure  Iteration  Safety  Factor  Plv  Mode 

1  156.34  ±45  S12 

2  213.89  0  XIT 


X  -  Direction  Compression 

Experimentally  Measured  Strength:  105.24±5.85  (Gillespie  1994) 

LAM3D  Failure  Analysis  Predictions: 

Failure  Iteration  Safety  Factor  Ply  Mode 

1  88.90  0  XIC 


Y  -  Direction  Tension 


Experimentally  Measured  Strength:  42.39+0.92  (Gillespie  1994) 


LAM3D  Failure  Analysis  Predictions; 
Failure  Iteration  Safety  Factor 

Ply 

Mode 

1 

33.21 

0 

X2T 

2 

24.79 

±45 

S12 

3 

36.89 

±45 

X2T 

4 

73.63 

0 

XIC 

Y  -  Direction  Compression 

Experimentally  Measured  Strength:  39.92±0.72  (Gillespie  1994) 

LAM3D  Failure  Analysis  Predictions: 

Failure  Iteration  Safety  Factor  Ply  Mode 

1  66.57  ±45  XIC 
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MATERIAL  DATA  BASE 


1 

1 

3 


REGIONS  TO  ANALYZE 

NUMBER  OF  REGIONS  IN  THE  DATABASE 

NUMBER  OF  MATERIALS  IN  DATABASE 


1.0, 

2.0, 

5.0, 

0.0 

! 

REG  ID,  FAIL  CRT, 

ITERS, 

ITEXT 

0.0, 

0.0, 

0.0, 

0.0 

! 

BETA,  PHI, 

SI, 

— 

8.0, 

1.0, 

0.0, 

0.0 

! 

NPLY,  REG  TYP, 

— 

2.0, 

0.005, 

0.0, 

-200.0 

! 

MAT  ID,  THICKNESS, 

ANGLE, 

TEMP 

2.0, 

0.005, 

0.0, 

-200.0 

! 

MAT  ID,  THICKNESS, 

ANGLE, 

TEMP 

2.0, 

0.005, 

45.0, 

-200.0 

f 

MAT  ID,  THICKNESS, 

ANGLE, 

TEMP 

2.0, 

0.005, 

-45.0, 

-200.0 

! 

MAT  ID,  THICKNESS, 

ANGLE, 

TEMP 

2.0, 

0.005, 

-45.0, 

-200.0 

! 

MAT  ID,  THICKNESS, 

ANGLE, 

TEMP 

2.0, 

.  0.005, 

45.0, 

-200.0 

! 

MAT  ID,  THICKNESS, 

ANGLE, 

TEMP 

2.0, 

0.005, 

0.0, 

-200.0 

! 

MAT  ID,  THICKNESS, 

ANGLE, 

TEMP 

2.0, 

0.005, 

0.0, 

-200.0 

! 

MAT  ID,  THICKNESS, 

ANGLE, 

TEMP 

1.0, 

1.0, 

O.lOOE+00 

lO.OOOE+06, 

lO.OOOE+06, 

lO.OOOE+06 

0.280E+00, 

0.280E+00, 

0.280E+00 

3.906E+06, 

3.906E+06, 

3.906E+06 

0.0, 

0.0, 

0.0 

80.000E+03, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

0.0, 

0.0, 

0.0 

MATJD,  MAT_TYPE,  MAT_DENS(A1) 


El, 

E2, 

E3 

NU23, 

NU13, 

NU12 

G23, 

G13, 

G12 

Al, 

A2, 

A3 

(l)VMl, 

“"“9 

—  (VonMises) 

“"“5 

”’“'“9 

— 

*""”9 

— 

““*"9 

'‘""'9 

— 

"""9 

— 

(2)X1T, 

X2T, 

X3T  (MaxStress) 

XIC, 

X2C, 

X3C 

X23, 

X13, 

X12 

*""■9 

"""9 

— 

“*”9 

""“9 

— 

(3)  YIT, 

Y2T, 

Y3T  (MaxStrain) 

YIC, 

Y2C, 

Y3C 

Y23, 

Y13, 

Y12 

—9 

“"“9 

— 

"""9 

“"“9 

— 

(4)  XIT, 

X2T, 

HPD  (Hydro-Pressure) 

XIC(O), 

X2C(0), 

S 

MLl, 

ML2, 

LTP 

MTl, 

MT2, 

TTP 

MSI, 

MS2, 

STP 
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0.0, 

0.0, 

0.0 

!  (5)X1T, 

XIC, 

—  (Tsai-Wu) 

0.0, 

0.0, 

0.0 

!  X2T, 

X2C, 

— 

0.0, 

0.0, 

0.0 

!  S12, 

F12, 

F23 

0.0, 

0.0, 

0.0 

•  “""9 

""“9 

— 

0.0, 

0.0, 

0.0 

«  ”"“5 

“”"9 

— 

0.0, 

0.0, 

0.0 

!  (6)E1, 

E2, 

E3  (Christensen's) 

0.0, 

0.0, 

0.0 

!  NU23, 

NU13, 

NU12 

0.0, 

0.0, 

0.0 

!  G23, 

G13, 

G12 

0.0, 

0.0, 

0.0 

!  K, 

ALPHA, 

— 

0.0, 

0.0, 

0.0 

!  YIT, 

YIC, 

— 

0.0, 

0.0, 

0.0 

!  (7)  Al, 

All, 

™  (Feng's) 

0.0, 

0.0, 

0.0 

!  A2, 

A4, 

— 

0.0, 

0.0, 

0.0 

!  A5, 

A55, 

— 

0.0, 

0.0, 

0.0 

““”9 

— 

0.0, 

0.0, 

0.0 

”””9 

— 

0.0, 

0.0, 

0.0 

!  (8)  XIT, 

X2T, 

X3T  (ModifiedHashin) 

0.0, 

0.0, 

0.0 

!  XIC, 

X2C, 

X3C 

0.0, 

0.0, 

0.0 

!  X23, 

X13, 

X12 

0.0, 

0.0, 

0.0 

““9 

0.0, 

0.0, 

0.0 

•  “"“9 

"“"9 

— 

2.0,  1.0, 

1.501E-04 

!  MATJD, 

MAT_TYPE,  MAT_DENS  (Graphitc/8551-  Epoxy) 

2.230E+06, 

1.220E+06, 

1.220E+06 

J 

El,  E2, 

E3 

0.450E+00, 

0.330E+00, 

0.330E+00 

f 

NU23,  NU13, 

NU12 

0.700E+06, 

0.700E+06, 

0.700E+06 

! 

G23,  G13, 

G12 

-5.000E-07, 

1.550E-05, 

1.550E-05 

! 

Al,  A2, 

A3 

0.0, 

0.0, 

0.0 

j 

(l)VMl, 

—  (VonMises) 

0.0, 

0.0, 

0.0 

t 

“““9  — — ^ 

— 

0.0, 

0.0, 

0.0 

J 

"■”"9  """9 

— 

0.0, 

0.0, 

0.0 

! 

-—9  — 

— 

0.0, 

0.0, 

0.0 

! 

—9  -  ^ 

— 

420.0E+03, 

8.7E+03, 

8.7E+03 

! 

(2)  XIT,  X2T, 

X3T  (MaxStress) 

175.0E+03, 

30.0E+03, 

30.0E+03 

! 

XIC,  X2C, 

X3C 

20.0E+03, 

20.0E+03, 

20.0E+03 

? 

X23,  X13, 

X12 

0.0, 

0.0, 

0.0 

f 

'""“9  """9 

— 

0.0, 

0.0, 

0.0 

! 

"”9  “““9 

— 

0.0188, 

0.0071, 

0.0071 

! 

(3)  YIT,  Y2T, 

Y3T  (MaxStrain) 

0.0078, 

0.0246, 

0.0246 

! 

YIC,  Y2C, 

Y3C 

0.0290, 

0.0290, 

0.0290 

1 

Y23,  Y13, 

Y12 

0.0, 

0.0, 

0.0 

! 

— 

0.0, 

0.0, 

0.0 

! 

-.—9  **"”9 

... 

420.0E+03, 

8.7E+03, 

1.0 

f 

(4)  XIT,  X2T, 

HPD  (Hydro-Pressure) 

175.0E+03, 

lOO.OE+03, 

20.30E+03 

1 

XIC(O),  X2C(0), 

S 

2.0, 

2.0, 

0.0 

j 

MLl,  ML2, 

LTP 

1.0, 

1.0, 

0.0 

1 

MTl,  MT2, 

TTP 
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1.0, 

1.0, 

0.0 

j 

MSI, 

MS2, 

STP 

420.0E+03, 

175.0E+03, 

0.0 

! 

(5)  XIT, 

XIC, 

—  (Tsai-Wu) 

8.7E+03, 

30.0E+03, 

0.0 

! 

X2T, 

X2C, 

— 

20.0E+03, 

1.1916E-09, 

-1.916E-09 

! 

S12, 

F12, 

F23 

0.0, 

0.0, 

0.0 

j 

— 

0.0, 

0.0, 

0.0 

? 

””"’9 

2.230E+06, 

1.220E+06, 

1.220E+06 

f 

(6)  El, 

E2, 

E3  (Christensen's) 

0.450E+00, 

0.330E+00, 

0.330E+00 

1 

NU23, 

NU13, 

NU12 

0.700E+06, 

0.700E+06, 

0.700E+06 

! 

G23, 

G13, 

G12 

1.507E-02, 

0.1028, 

0.0 

1 

K, 

ALPHA, 

— 

0.0188, 

0.0078, 

0.0 

f 

YIT, 

YIC, 

— 

452.55, 

101910.0, 

0.0 

j 

(7)  Al, 

All, 

—  (Feng's) 

1189.1, 

1189.1, 

0.0 

j 

A2, 

A4, 

— 

-37.42, 

4819.74, 

0.0 

! 

A5, 

A55, 

— 

0.0, 

0.0, 

0.0 

f 

“""9 

0.0, 

0.0, 

0.0 

! 

""“9 

— 

420.0E+03, 

8.7E+03, 

8.7E+03 

! 

(8)  XIT, 

X2T, 

X3T  (ModifiedHashin) 

175.0E+03, 

30.0E+03, 

30.0E+03 

f 

XIC, 

X2C, 

X3C 

20.0E+03, 

20.0E+03, 

20.0E+03 

! 

X23, 

X13, 

X12 

0.0, 

0.0, 

0.0 

1 

"""’"'9 

... 

0.0, 

0.0, 

0.0 

j 

**'"5 

”““9 

... 

3.0,  1.0, 

1.889E-04  ! 

MAT  ID, 

MAT 

'_TYPE,  MAT_DENS(S2-Glass/3501  -Epoxy) 

7.150E+06, 

2.130E+06, 

2.130E+06 

! 

El,  E2, 

E3 

0.499E+00, 

0.306E+00, 

0.296E+00 

! 

NU23,  NU13, 

NU12 

0.710E+06, 

0.980E+06, 

0.980E+06 

! 

G23,  G13, 

G12 

2.300E-06, 

1.850E-05, 

1.850E-05 

f 

Al,  A2, 

A3 

0.0, 

0.0, 

0.0 

(l)VMl,  -, 

—  (VonMises) 

0.0, 

0.0, 

0.0 

! 

—9  ""”9 

0.0, 

0.0, 

0.0 

f 

“““9  — ^ 

... 

0.0, 

0.0, 

0.0 

! 

"*””9  ”**"9 

0.0, 

0.0, 

0.0 

! 

—-9  “““9 

243.0E+03, 

7.0E+03, 

7.0E+03 

! 

(2)  XIT,  X2T, 

X3T  (MaxStress) 

177.0E+03, 

30.6E+03, 

35.0E+03 

! 

XIC,  X2C, 

X3C 

15.7E+03, 

17.0E+03, 

15.7E+03 

! 

X23,  X13, 

X12 

0.0, 

0.0, 

0.0 

! 

”*"9  ^ 

0.0, 

0.0, 

0.0 

f 

"’**9  """"9 

«... 

0.0340, 

0.0033, 

0.0040 

! 

(3)  YIT,  Y2T, 

Y3T  (MaxStrain) 

0.0248, 

0.0144; 

0.0164 

f 

YIC,  Y2C, 

Y3C 

0.0221, 

0.0174, 

0.0160 

! 

Y23,  Y13, 

Y12 

0.0, 

0.0, 

0.0 

! 

**"*9  *'""9 

0.0, 

0.0, 

0.0 

! 

-—9 

— 

0.0, 

0.0, 

0.0 

! 

(4)  XIT,  X2T, 

HPD  (Hydro-Pressure) 

0.0, 

0.0, 

0.0 

! 

XIC(O),  X2C(0), 

S 

0.0, 

0.0, 

0.0 

! 

MLl,  ML2, 

LTP 
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0.0, 

0.0, 

0.0 

!  MTl, 

MT2, 

TIP 

0.0, 

0.0, 

0.0 

!  MSI, 

MS2, 

STP 

0.0, 

0.0, 

0.0 

!  (5)  XIT, 

XIC, 

—  (Tsai-Wu) 

0.0, 

0.0, 

0.0 

!  X2T, 

X2C, 

— 

0.0, 

0.0, 

0.0 

!  S12, 

F12, 

F23 

0.0, 

0.0, 

0.0 

— 

0.0, 

0.0, 

0.0 

““"9 

0.0, 

0.0, 

0.0 

!  (6)  El, 

E2, 

E3  (Christensen's) 

0.0, 

0.0, 

0.0 

!  NU23, 

NU13, 

NU12 

0.0, 

0.0, 

0.0 

!  G23, 

G13, 

G12 

0.0, 

0.0, 

0.0 

!  K, 

ALPHA, 

— 

0.0, 

0.0, 

0.0 

!  YIT, 

YIC, 

— 

0.0, 

0.0, 

0.0 

!  (7)  Al, 

All, 

—  (Feng's) 

0.0, 

0.0, 

0.0 

!  A2, 

A4, 

— 

0.0, 

0.0, 

0.0 

!  A5, 

A55, 

— 

0.0, 

0.0, 

0.0 

•  “""5 

"“”5 

— 

0.0, 

0.0, 

0.0 

.  — 

*"'"9 

— 

0.0, 

0.0, 

0.0 

!  (8)  XIT, 

X2T, 

X3T  (ModifiedHashin) 

0.0, 

0.0, 

0.0 

!  XIC, 

X2C, 

X3C 

0.0, 

0.0, 

0.0 

!  X23, 

X13, 

X12 

0.0, 

0.0, 

0.0 

“**”9 

— 

0.0, 

0.0, 

0.0 

•  “”5 

"""9 

— 
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SAMPLE  OUTPUT  FILE 


FROM  OUTOFIN.F . 

,tiitiiHt**itiiHH^i^^t»iilritiHtiltiHi4‘*****************************’l' 

*  OUTPUT  FROM  LAM3D  (VERSION  3.0)  * 

***********  **iti4i*«i|i*4i*4i****4c*4ii|i**it<>K****>l>>l>***i|i«i|c****i|i 


TOTAL  NUMBER  OF  PLIES  >  8 

TYPE  OF  LAMINA  CONSTITUTIVE  MODEL  SELECTED  > 

(1)  ORTHOTROPIC  (9  I.E.C.) 

TYPE  OF  LAMINATE  CONSTITUTIVE  MODEL  SELECTED  > 

(2)  EXACT  (REF.  CHOU) 

TYPE  OF  FAILURE  CRITERIA  > 

(2)  MAXIMUM  STRESS 

NUMBER  OF  ITERATIONS  FOR  FAILURE  > 

10 

*********************** 

3-D  LAMINA  PROPERTIES 
*********************** 


PLY=  1 


El 

E2 

E3 

0.223E+08 

0.122E+07 

0.122E+07 

G23 

G13 

G12 

0.700E+06 

0.700E+06 

0.700E+06 

NU23 

NU13 

NU23 

0.450 

0.330 

0.330 

PLY=  2 

El 

E2 

E3 

0.223E+08 

0.122E+07 

0.122E+07 

G23 

G13 

G12 

0.700E+06 

0.700E+06 

0.700E+06 

NU23 

NU13 

NU23 

0.450 

0.330 

0.330 

PLY=  3 


El 

E2 

E3 

0.223E+08 

0.122E+07 

0.122E+07 

G23 

G13 

G12 

0.700E+06 

0.700E+06 

0.700E+06 

NU23 

NU13 

NU23 

0.450 

0.330 

0.330 

PLY=  4 

El 

E2 

E3 

0.223E+08 

0.122E+07 

0.122E+07 

G23 

G13 

G12 

0.700E+06 

0.700E+06 

0.700E+06 

NU23 

NUB 

NU23 

0.450 

0.330 

0.330 
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PLY=  5 


El 

0.223E+08 

G23 

0.700E+06 

NU23 

0.450 

E2 

0.122E+07 

G13 

0.700E-K)6 

NU13 

0.330 

E3 

0.122E+07 

G12 

0.700E+06 

NU23 

0.330 

PLY=  6 

El 

E2 

E3 

0.223E+08 

0.122E+07 

0.122E+07 

G23 

G13 

G12 

0.700E+06 

0.700E+06 

0.700E+06 

NU23 

NU13 

NU23 

0.450 

0.330 

0.330 

PLY=  7 

El 

E2 

E3 

0.223E+08 

0.122E+07 

0.122E+07 

G23 

G13 

G12 

0.700E+06 

0.700E+06 

0.700E+06 

NU23 

NU13 

NU23 

0.450 

0.330 

0.330 

PLY=  8 

El 

E2 

E3 

0.223E+08 

0.122E+07 

0.122E+07 

G23 

G13 

G12 

0.700E+06 

0.700E+06 

0.700E+06 

NU23 

NU13 

NU23 

0.450 

0.330 

0.330 

« l»  *******  i|i  ***  Ik  Hi  * 

FAILURE  PARAMETERS 

**************  Ik  **  mm  Ik 


PLY=  1 

FR-1,1  FR-1,2  FR-1,3 

0.420E+06  0.870E+04  0.870E+04 
FR-2,1  FR-2,2  FR-2,3 

0.175E+06  0.300E+05  0.300E+05 
FR-3,1  FR-3,2  FRO, 3 

0.200E+05  0.200E+05  0.200E+05 
FR-4,1  FR.4,2  FR-4,3 

O.OOOE+00  O.OOOE+00  O.OOOE+00 
FR-5,1  FR-5,2  FR-5,3 

O.OOOE+00  O.OOOE+00  O.OOOE+00 


PLY=  2 
FR-1,1  FR-1,2  FR-1,3 

0.420E+06  0.870E+04  0.870E+04 
FR-2,1  FR-2,2  FR-2,3 

0.175E+06  0.300E+05  0.300E+05 
FR-3,1  FR-3,2  FR-3,3 

0.200E+05  0.200E+05  0.200E+05 
FR-4,1  FR-4,2  FR-4,3 
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O.OOOE+00  O.OOOE+00  O.OOOE+00 
FR-5,1  FR-5,2  FR-5,3 

O.OOOE+00  O.OOOE+00  O.OOOE+00 


PLY=  3 
FR-1,1  FR-1,2  FR-1,3 

0.420E+06  0.870E+04  0.870E+04 
FR-2,1  FR-2,2  FR-2,3 

0.175E+06  0.300E+05  0.300E+05 
FR-3,1  FR-3,2  FR-3,3 

0.200E+05  0.200E+05  0.200E+05 
FR-4,1  FR-4,2  FR-4,3 

O.OOOE+00  O.OOOE+00  O.OOOE+00 
FR-5,1  FR-5,2  FR-5,3 

O.OOOE+00  O.OOOE+00  O.OOOE+00 


PLY=  4 

FR-1,1  FR-1,2  FR-1,3 

0.420E+06  0.870E+04  0.870E+04 
FR-2,1  FR-2,2  FR-2,3 

0.175E+06  0.300E+05  0.300E+05 
FR-3,1  FR-3,2  FR-3,3 

0.200E+05  0.200E+05  0.200E+05 
FR-4,1  FR-4,2  FR-4,3 

O.OOOE+00  O.OOOE+00  O.OOOE+00 
FR-5,1  FR-5,2  FR-5,3 

O.OOOE+00  O.OOOE+00  O.OOOE+00 


PLY=  5 
FR-1,1  FR-1,2  FR-1,3 

0.420E+06  0.870E+04  0.870E+04 
FR-2,1  FR-2,2  FR-2,3 

0.175E+06  0.300E+05  0.300E-K)5 
FR-3,1  FR-3,2  FR-3,3 

0.200E+05  0.200E+05  0.200E+05 
FR-4,1  FR-4,2  FR-4,3 

O.OOOE+00  O.OOOE+00  O.OOOE+00 
FR-5,1  FR-5,2  FR-5,3 

O.OOOE+00  O.OOOE+00  O.OOOE+OO 


PLY=  6 
FR-1,1  FR-1,2  FR-1,3 

0.420E+06  0.870E+04  0.870E+04 
FR-2,1  FR-2,2  FR-2,3 

0.175E+06  0.300E+05  0.300E+05 
FR-3,1  FR-3,2  FR-3,3 

0.200E+05  0.200E+05  0.200E+05 
FR-4,1  FR-4,2  FR-4,3 

O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 
FR-5,1  FR-5,2  FR-5,3 

O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 


PLY=  7 
FR-1,1  FR-1,2  FR-1,3 

0.420E+06  0.870E+04  0.870E+04 
FR-2,1  FR-2,2  FR-2,3 

0.175E+06  0.300E+05  0.300E+05 
FR-3,1  FR-3,2  FR-3,3 


0.200E+05  0.200E+05  0.200E+05 
FR-4,1  FR-4,2  FR-4,3 

O.OOOE+00  O.OOOE+00  O.OOOE+00 
FR-5,1  FR-5,2  FR-5,3 

O.OOOE+00  O.OOOE+00  O.OOOE+OO 


PLY=  8 
FR-1,1  FR-1,2  FR-1,3 

0.420E+06  0.870E+04  0.870E+04 
FR-2,1  FR-2,2  FR-2,3 

0.175E+06  0.300E+05  0.300E+05 
FR-3,1  FR-3,2  FR-3,3 

0.200E+05  0.200E+05  0.200E+05 
FR-4,1  FR-4,2  FR-4,3 

O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 
FR-5,1  FR-5,2  FR-5,3 

O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 


*  *  ]|l  *  *  4l  41  *  *  *  3)1  *  *  *  «  *  *  *  ||I 

LAMINATE  GEOMETRY 

«#3|C3)C*3|I3)|3|B«3)||»3»3)33)13)B3|33|33)33)3 

PLY  ANGLE  (DEGREES)  THICKNESS  (IN)  TEMP  (F) 


1 

O.OOOOE+00 

0.5000E-02 

-200.0 

2 

O.OOOOE+00 

0.5000E-02 

-200.0 

3 

45.00 

0.5000E-02 

-200.0 

4 

-45.00 

0.5000E-02 

-200.0 

5 

-45.00 

0.5000E-02 

-200.0 

6 

45.00 

0.5000E-02 

-200.0 

7 

O.OOOOE+00 

0.5000E-02 

-200.0 

8 

O.OOOOE+00 

0.5000E-02 

-200.0 

3)>  3*  3)3  3)3  3)3  3*  3)3  3»  ])[  1)3  1)3  1)1  *  3)t  3)1 3)31)3  3»  3)3  3)1 3)1 «  jH  3)3  *  3)t  3)1 3)1 3)3  3)1  ]|t  3)3  IK  3)1 3|t «  «  3)1 3)1 3)1 3)C  *  3)1 

LAMINATE  ORIENTATION  WITHIN  GLOBAL  SYSTEM 

3)3  3)1 3)3  1)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3»  3)3  3)3  »  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3)1 1)3  3)3  ««  3)3  3)3  3)3  3)3  1)3  3)3 

BETA=  O.OOOOOOOE+00 
PHI  =  O.OOOOOOOE+00 

SI  =  O.OOOOOOOE+00 

It  1»  1)3  3)3  3)3  1)3  3)3  3)3  3)3  3)3  3)3  3)3  1)3  3)3 1)3  3»  3)3  1)3  3)3  43  3)3  3)3  1)3 1»  1)3  3)3  3)3 1)3  3)3 1)3  3)3  3)3  3)3  3)3  3)3  1)3 

APPLIED  GLOBAL  MECHANICAL  STRESSES 

3)3  3)3  1)3  3)3  3)3  1)3  3)3  3)3  3)3  3)3  3)3  3)3  3)3  3t  3*  3)3  3)3  3)3  1)3  3)3  3)3  3)3  3)3  3)3  3t  *  3)3  3)3  3*  3)3  3)3  3)3  3)3  1)3  3)3  3)3 

X  =  -1000.000 
Y  =  O.OOOOOOOE+00 
Z  =  O.OOOOOOOE+00 
YZ  =  10.00000 

XZ  =  35.00000 

XY  =  10.00000 

43  3)3  1)3  3)3 1)3 1)3  3)3  1)3  *  3)3  1)3  1)3  3)3  1)3  3)3  3)3  3)3  1)3  3)3  1)3  3)3  1)3  1)3  3)3  3)3  3)3  3)3  3)3  3t  3t  3)3  3)3  3)3  1)3  3)3 1)3 

APPLIED  LOCAL  MECHANICAL  STRESSES 

433)33)33)33)3  4343  43433)34b4>4<4s4>433)34B3)3433)343433)33)3433)33)34s43  3)33)3  3)33)3433)3 


X  = 

-1000.000 

Y  = 

O.OOOOOOOE+00 

Z  = 

O.OOOOOOOE+00 

YZ  = 

10.00000 

XZ  = 

35.00000 

XY  = 

10.00000 

60 


FROM  OUTCIJS.F . 

CIJ  PRINCIPAL  MATRIX  > 

PLY:  1 

0.2279E+08  0.7482E+06  0.7482E+06 
0.7482E+06  0.1554E+07  0.7130E+06 
0.7482E+06  0.7130E+06  0.1554E+07 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
PLY:  2 

0.2279E+08  0.7482E+06  0.7482E+06 
0.7482E+06  0.1554E+07  0.7130E+06 
0.7482E+06  0.7130E+06  0.1554E+07 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
PLY:  3 

0.2279E+08  0.7482E+06  0.7482E+06 
0.7482E+06  0.1554E+07  0.7130E+06 
0.7482E+06  0.7130E+06  0.1554E+07 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
PLY:  4 

0.2279E+08  0.7482E+06  0.7482E+06 
0.7482E+06  0.1554E+07  0.7130E+06 
0.7482E+06  0.7130E+06  0.I554E+07 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
PLY:  5 

0.2279E+08  0.7482E+06  0.7482E+06 
0.7482E+06  0.I554E+07  0.7130E+06 
0.7482E+06  0.7130E+06  0.1554E+07 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
PLY:  6 

0.2279E+08  0.7482E+06  0.7482E+06 
0.7482E+06  0.1554E+07  0.7130E+06 
0.7482E+06  0.7130E+06  0.1554E+07 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
PLY:  7 

0.2279E+08  0.7482E+06  0.7482E+06 
0.7482E+06  0.I554E+07  0.7130E+06 
0.7482E+06  0.7130E+06  0.1554E+07 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
PLY:  8 

0.2279E+08  0.7482E+06  0.7482E+06 
0.7482E+06  0.1554E+07  0.7I30E+06 
0.7482E+06  0.7130E+06  0.1554E+07 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 


O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
0.7000E+06  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  0.7000E+06  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  0.7000E+06 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
0.7000E-K)6  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  0.7000E+06  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  0.7000E+06 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
0.7000E+06  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  0.7000E+06  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  0.7000E+06 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
0.7000E+06  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  0.7000E+06  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  0.7000E+06 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
0.7000E+06  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  0.7000E+06  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  0.7000E+06 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
0.7000E+06  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  0.7000E+06  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  0.7000E+06 

O.OOOOE+OO  O.OOOOE+00  O.OOOOE+OO 
O.OOOOE+00  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  0.7000E+06 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 


O.OOOOE+00 

CIJ  GLOBAL 
PLY:  1 

0.2279E+08 
0.7482E+06 
0.7482E+06 
O.OOOOE+00 
O.OOOOE+00 
O.OOOOE+00 
PLY:  2 

0.2279E+08 
0.7482E+06 
0.7482E+06 
O.OOOOE+00 
O.OOOOE+00 
O.OOOOE+00 
PLY:  3 

0.7161E+07 
0.5761E+07 
0.7306E+06 
O.OOOOE+00 
O.OOOOE+00 
0.5310E+07 
PLY:  4 

0.7161E+07 
0.5761E+07 
0.7306E+06 
O.OOOOE+00 
O.OOOOE+00 
-0.5310E+07 
PLY:  5 

0.7161E+07 
0.5761E+07 
0.7306E+06 
O.OOOOE+00 
O.OOOOE+00 
-0.5310E+07 
PLY:  6 

0.7161E+07 
0.5761E+07 
0.7306E+06 
O.OOOOE+00 
O.OOOOE+00 
0.5310E+07 
PLY:  7 

0.2279E+08 
0.7482E+06 
0.7482E+06 
O.OOOOE+00 
O.OOOOE+00 
O.OOOOE+00 
PLY:  8 

0.2279E+08 
0.7482E+06 
0.7482E+06 
O.OOOOE+00 
O.OOOOE+00 


O.OOOOE+00 

MATRIX  > 

0.7482E+06 

0.1554E+07 

0.7130E+06 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7482E+06 

0.1554E+07 

0.7130E+06 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.576  lE+07 
0.7161E+07 
0.7306E+06 
O.OOOOE+00 
O.OOOOE+00 
0.5310E+07 

0.576  lE+07 
0.7161E+07 
0.7306E+06 
O.OOOOE+00 
O.OOOOE+00 
-0.5310E+07 

0.576  lE+07 
0.716  lE+07 
0.7306E+06 
O.OOOOE+00 
O.OOOOE+00 
-0.5310E+07 

0.576  lE+07 
0.7161  E+07 
0.7306E+06 
O.OOOOE+00 
O.OOOOE+00 
0.5310E+07 

0.7482E+06 

0.1554E+07 

0.7130E+06 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7482E+06 

0.1554E+07 

0.7130E+06 

O.OOOOE+00 

O.OOOOE+00 


O.OOOOE+00 


0.7482E+06 

0.7130E+06 

0.1554E+07 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7482E+06 

0.7130E+06 

0.1554E+07 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7306E+06 

0.7306E+06 

0.1554E+07 

O.OOOOE+00 

O.OOOOE+00 

0.1762E+05 

0.7306E+06 

0.7306E+06 

0.1554E+07 

O.OOOOE+00 

O.OOOOE+00 

-0.1762E+05 

0.7306E+06 

0.7306E+06 

0.1554E+07 

O.OOOOE+00 

O.OOOOE+00 

-0.1762E+05 

0.7306E+06 

0.7306E+06 

0.1554E+07 

O.OOOOE+00 

O.OOOOE+00 

0.1762E+05 

0.7482E+06 

0.7130E+06 

0.1554E+07 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7482E+06 

0.7130E+06 

0.1554E+07 

O.OOOOE+00 

O.OOOOE+00 


O.OOOOE+00 


O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7000E+06 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7000E+06 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7000E+06 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7000E+06 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7000E+06 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7000E+06 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7000E+06 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 


O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 


0.7000E+06 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

0.5310E+07 

0.5310E+07 

0.1762E+05 

O.OOOOE+OO 

O.OOOOE+OO 

0.5713E+07 

-0.5310E+07 

-0.5310E+07 

-0.1762E+05 

O.OOOOE+OO 

O.OOOOE+OO 

0.5713E+07 

-0.5310E+07 

-0.5310E+07 

-0.1762E+05 

O.OOOOE+OO 

O.OOOOE+OO 

0.5713E+07 

0.5310E+07 

0.5310E+07 

0.1762E+05 

O.OOOOE+OO 

O.OOOOE+OO 

0.5713E+07 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 


62 


O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06 


ALPHA  GLOBAL  MATRIX  > 

PLY:  1 

-0.5000E-06  0.1550E-04  0.1550E-04  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

PLY:  2 

-0.5000E-06  0.1550E-04  0.1550E-04  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

PLY:  3 

0.7500E-05  0.7500E-05  0.1550E-04  O.OOOOE+00  O.OOOOE+00  -0.1600E-04 

PLY:  4 

0.7500E-05  0.7500E-05  0.1550E-04  O.OOOOE+00  O.OOOOE+00  0.1600E-04 

PLY:  5 

0.7500E-05  0.7500E-05  0.1550E-04  O.OOOOE+00  O.OOOOE+00  0.1600E-04 

PLY:  6 

0.7500E-05  0.7500E-05  0.1550E-04  O.OOOOE+00  O.OOOOE+00  -0.1600E-04 

PLY:  7 

-0.5000E-06  0.1550E-04  0.1550E-04  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

PLY:  8 

-0.5000E-06  0.1550E-04  0.1550E-04  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

CM  GLOBAL  STIFFNESS  MATRIX  > 

0.1498E+08  0.3255E+07  0.7394E-K)6  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.3255E+07  0.4358E+07  0.7218E+06  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.7394E+06  0.7218E+06  0.1554E+07  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+OO  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.3206E+07 

FROM  OUTEFPR.F . 

EFFECTIVE  LAMINATE  PROPERTIES  > 

EX  =  1.2518537E+07 
EY  =  3442994. 

EZ  =  1431594. 

NUYZ  =  0.3696677 
NUXZ  =  0.1396095 
NUXY  =  0.7237686 
NUZY  =  0.1537075 
NUZX  =  1.5965447E-02 
NUYX  =  0.1990592 
GYZ  =  699999.9 
GXZ  =  699999.9 
GXY  =  3206368. 

ALPHAX  = -6.64765 14E-07 
ALPHAY  =  3.6833085E-06 
ALPHAZ  =  2.0974869E-05 
ALPHA  YZ=  O.OOOOOOOE+00 
ALPHAXZ=  O.OOOOOOOE+00 
ALPHAXY  =  -2.9743138E-13 

FROM  OUTLOAD.F . 

AVERAGE  MECHANICAL  STRESSES  (PSI.)  > 

X  =  -1000.000 
Y  =  O.OOOOOOOE+00 
Z  =  O.OOOOOOOE+00 
YZ  =  10.00000 

XZ  =  35.00000 

XY  =  10.00000 
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AVERAGE  THERMAL  STRESSES  (PSI.)  > 
X-TH  =  -3508.091 
Y-TH  =  -2.0268 120E-04 
Z-TH  =  2.0282280E-04 
YZ-TH  =  3.9123173E-05 
XZ-TH  =  5.5890254E-11 
XY-TH  =  7.9843231E-17 

MECHANICAL  LAMINATE  STRAINS  > 
EPSX  =  .7.9881538E-05 
EPSY  =  5.7815741E-05 
EPSZ  =  1.1152222E-05 
EPSYZ=  1. 42857 16E-05 
EPSXZ=  5.0000006E-05 
EPSXY=  3.1187940E-06 

THERMAL  LAMINATE  STRAINS  > 
EPSX-TH  =  -2.0268 120E-04 
EPSY-TH  =  2.0282280E-04 
EPSZ-TH  =  3.9123 173E-05 
EPSYZ-TH  =  5.5890254E-1 1 
EPSXZ-TH=  7.9843231E-I7 
EPSXY-TH  =  2.4901459E-23 

FROM  OUTMESTR.F . 


MECHANICAL  GLOBAL  PLY  STRAIN 


PLY 

EPS-X 

EPS-Y 

EPS-Z 

EPS-YZ 

EPS-XZ 

EPS-XY 

I 

-.799E-04 

0.578E-04 

0.119E-04 

0.143E-04 

0.500E-04 

0.312E-05 

2 

-.799E-04 

0.578E-04 

0.119E-04 

0.143E-04 

0.500E-04 

0.312E-05 

3 

-.799E-04 

0.578E-04 

0.103E-04 

0.143E-04 

0.500E-04 

0.312E-05 

4 

-.799E-04 

0.578E-04 

0.104E-04 

0.143E-04 

0.500E-04 

0.312E-05 

5 

-.799E-04 

0.578E-04 

0.104E-04 

0.143E-04 

0.500E-04 

0.312E-05 

6 

-.799E-04 

0.578E.04 

0.103E-04 

0.143E-04 

0.500E-04 

0.312E-05 

7 

-.799E-04 

0.578E-04 

0.I19E-04 

0.143E-04 

0.500E-04 

0.312E-05 

8 

-.799E-04 

0.578E-04 

0.I19E-04 

0.143E-04 

0.500E-04 

0.312E-05 

MECHANICAL  GLOBAL  PLY  STRESS 


PLY 

SIG-X  SIG-Y 

SIG-Z 

SIG-YZ  SIG-XZ  SIG-XY 

I 

-.I77E+04  38.6 

O.OOOE+00 

10.0 

35.0 

2.18 

2 

-.177E+04  38.6 

O.OOOE+00 

10.0 

35.0 

2.18 

3 

-215.  -22.1 

O.OOOE+00 

10.0 

35.0 

-99.2 

4 

-248.  -55.1 

O.OOOE+00 

10.0 

35.0 

135. 

5 

-248.  -55.1 

O.OOOE+00 

10.0 

35.0 

135. 

6 

-215.  -22.1 

O.OOOE+00 

10.0 

35.0 

-99.2 

7 

-.177E+04  38.6 

O.OOOE+00 

10.0 

35.0 

2.18 

8 

-.177E+04  38.6 

O.OOOE+00 

lO.O 

35.0 

2.18 

MECHANICAL  PRINCIPAL  PLY  STRAIN 


PLY 

EPS-1 

EPS-2 

EPS-3 

EPS-23 

EPS-13 

EPS- 12 

1 

-.799E-04 

0.578E-04 

0.119E-04 

0.143E-04 

0.500E-04 

0.312E-05 

2 

-.799E-04 

0.578E-04 

0.119E-04 

0.143E-04 

0.500E-04 

0.312E-05 

3 

-.947E-05 

-.126E-04 

0.103E-04 

-.253E-04 

0.455E-04 

0.138E-03 

4 

-.126E-04 

-.947E-05 

0.104E-04 

0.455E-04 

0.253E-04 

-.138E-03 

5 

-.126E-04 

-.947E-05 

0.104E-04 

0.455E-04 

0.253E-04 

-.138E-03 

6 

-.947E-05 

-.126E-04 

0.103E-04 

-.253E-04 

0.455E-04 

0.138E-03 

7 

-.799E-04 

0.578E-04 

0.119E-04 

0.143E-04 

0.500E-04 

0.312E-05 

8 

-.799E-04 

0.578E-04 

0.119E-04 

0.143E-04 

0.500E-04 

0.312E-05 
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MECHANICAL  PRINCIPAL  PLY  STRESS 


PLY 

SIG-1  SIG-2 

SIG-3 

SIG-23 

SIG-13 

i  SIG-: 

1 

-.177E+04  38.6 

O.OOOE+00 

10.0 

35.0 

2.18 

2 

-.177E+04  38.6 

O.OOOE+00 

10.0 

35.0 

2.18 

3 

-218.  -19.3 

O.OOOE+00  - 

■17.7 

31.8 

96.4 

4 

-286.  -16.7 

O.OOOE+00 

31.8 

17.7 

■96.4 

5 

-286.  -16.7 

O.OOOE+00 

31.8 

17.7 

■96.4 

6 

-218.  -19.3 

O.OOOE+00  - 

•17.7 

31.8 

96.4 

7 

-.177E+04  38.6 

O.OOOE+00 

10.0 

35.0 

2.18 

8 

-.177E+04  38.6 

O.OOOE+00 

1 0.0 

35.0 

2.18 

FROM  OUTTHSTR.F  . 

THERMAL  GLOBAL  PLY  STRAIN 

PLY  EPS-X  EPS-Y  EPS-Z  EPS-YZ  EPS-XZ  EPS-XY 

1  -.203E-03  0.203E-03  0.453E-05  0.559E-10  0.798E-I6  0.249E-22 

2  -.203E-03  0.203E-03  0.453E-05  0.559E-10  0.798E-16  0.249E-22 

3  -.203E-03  0.203E-03  -.664E-07  0.559E-10  0.798E-16  0.249E-22 

4  -.203E-03  0.203E-03  -.664E-07  0.559E-I0  0.798E-16  0.249E-22 

5  -.203E-03  0.203E-03  -.664E-07  0.559E-10  0.798E-16  0.249E-22 

6  -.203E-03  0.203E-03  -.664E-07  0.559E-10  0.798E-16  0.249E-22 

7  -.203E-03  0.203E-03  0.453E-05  0.559E-I0  0.798E-I6  0.249E-22 

8  -.203E-03  0.203E-03  0.453E-05  0.559E-10  0.798E-I6  0.249E-22 

THERMAL  GLOBAL  PLY  STRESS 

PLY  SIG-X  SIG-Y  SIG-Z  SIG-YZ  SIG-XZ  SIG-XY 

1  -.211E+04  0.712E+04  0.203E-03  0.391E-04  0.559E-10  0.174E-16 

2  -.211E+04  0.712E+04  0.203E-03  0.391E-04  0.559E-I0  0.174E-16 

3  0.437E+04  0.494E+04  0.203E-03  0.391E-04  0.559E-10  -.230E+04 

4  0.437E+04  0.494E+04  0.203E-03  0.391E-04  0.559E-10  0.230E+04 

5  0.437E+04  0.494E+04  0.203E-03  0.391E-04  0.559E-10  0.230E+04 

6  0.437E+04  0.494E+04  0.203E-03  0.391E-04  0.559E-10  -.230E+04 

7  -.211E+04  0.712E+04  0.203E-03  0.391E-04  0.559E-10  0.I74E-16 

8  -.211E+04  0.712E+04  0.203E-03  0.391E-04  0.559E-10  0.174E-16 

THERMAL  PRINCIPAL  PLY  STRAIN 

PLY  EPS-1  EPS-2  EPS-3  EPS-23  EPS-13  EPS-12 

1  -.203E-03  0.203E-03  0.453E-05  0.559E-10  0.798E-16  0.249E-22 

2  -.203E-03  0.203E-03  0.453E-05  0.559E-10  0.798E-16  0.249E-22 

3  0.708E-07  0.708E-07  -.664E-07  0.395E-10  0.395E-10  0.406E-03 

4  0.708E-07  0.708E-07  -.664E-07  0.395E-10  -.395E-10  -.406E-03 

5  0.708E-07  0.708E-07  -.664E-07  0.395E-10  -.395E-10  -.406E-03 

6  0.708E-07  0.708E-07  -.664E-07  0.395E-10  0.395E-10  0.406E-03 

7  -.203E-03  0.203E-03  0.453E-05  0.559E-10  0.798E-16  0.249E-22 

8  -.203E-03  0.203E-03  0.453E-05  0.559E-10  0.798E-16  0.249E-22 

THERMAL  PRINCIPAL  PLY  STRESS 

PLY  SIG-1  SIG-2  SIG-3  SIG-23  SIG-13  SIG-12 

1  -.211E+04  0.712E+04  0.203E-03  0.391E-04  0.559E-10  0.174E-16 

2  -.211E+04  0.712E+04  0.203E-03  0.391E-04  0.559E-10  0.174E-16 

3  0.236E+04  0.695E+04  0.203E-03  0.277E-04  0.277E-04  284. 

4  0.236E+04  0.695E+04  0.203E-03  0.277E-04  -.277E-04  -284. 

5  0.236E+04  0.695E+04  0.203E-03  0.277E-04  -.277E-04  -284. 

6  0.236E+04  0.695E-I-04  0.203E-03  0.277E-04  0.277E-04  284. 

7  -.211E+04  0.712E+04  0.203E-03  0.391E-04  0.559E-10  0.174E-16 

8  -.211E+04  0.712E+04  0.203E-03  0.391E-04  0.559E-10  0.174E-16 

FROM  OUTFAIL.F . 

ITERATION  NUMBER  >  1 
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FAILURE  ASSESSMENT- 
PLY  >  1  MODE  >  X2T 

PLY  >  2  MODE  >  X2T 

PLY  >  7  MODE  >  X2T 

PLY  >  8  MODE  >  X2T 

SAFTY  FACTOR  >  40.91054 

FROM  OUTCIJS.F . 

CIJ  PRINCIPAL  MATRIX  > 

PLY:  1 

0.2230E+08  -0.7333  -2.363  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

-0.7333  -0.4024E-05  -2.222  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

-2.363  -2.222  -4.938  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06 

PLY:  2 

0.2230E+08  -0.7333  -2.363  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

-0.7333  -0.4024E-05  -2.222  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

-2.363  -2.222  -4.938  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06 

PLY:  3 

0.2279E+08  0.7482E+06  0.7482E+06  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.7482E+06  0.1554E+07  0.7130E+06  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.7482E+06  0.7130E+06  0.1554E+07  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06 

PLY:  4 

0.2279E+08  0.7482E+06  0.7482E+06  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.7482E+06  0.1554E+07  0.7130E+06  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.7482E+06  0.7130E+06  0.1554E+07  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06 

PLY:  5 

0.2279E+08  0.7482E+06  0.7482E+06  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.7482E+06  0.I554E+07  0.7130E+06  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.7482E+06  0.7130E+06  0.1554E+07  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06 

PLY:  6 

0.2279E+08  0.7482E+06  0.7482E+06  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.7482E+06  0.1554E+07  0.7130E+06  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.7482E+06  0.7130E+06  0.1554E+07  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06 

PLY:  7 

0.2230E+08  -0.7333  -2.363  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

-0.7333  -0.4024E-05  -2.222  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

-2.363  -2.222  -4.938  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06 
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PLY:  8 

0.2230E+08  -0.7333  -2.363  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

-0.7333  -0.4024E-05  -2.222  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

-2.363  -2.222  -4.938  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00 

O.OOOOE+OO  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06 

CIJ  GLOBAL  MATRIX  > 

PLY:  1 

0.2230E+08  -0.7333  -2.363  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

-0.7333  -0.4024E-05  -2.222  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

-2.363  -2.222  -4.938  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06 

PLY:  2 

0.2230E+08  -0.7333  -2.363  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

-0.7333  -0.4024E-05  -2.222  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

-2.363  -2.222  -4.938  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06 

PLY:  3 

0.7161E+07  0.5761E+07  0.7306E+06  O.OOOOE+OO  O.OOOOE+OO  0.5310E+07 

0.5761E+07  0.7161E+07  0.7306E+06  O.OOOOE+OO  O.OOOOE+OO  0.5310E+07 

0.7306E+06  0.7306E+06  0.1554E+07  O.OOOOE+OO  O.OOOOE+OO  0.1762E+05 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 

0.5310E+07  0.5310E+07  0.1762E+05  O.OOOOE+OO  O.OOOOE+OO  0.5713E+07 

PLY:  4 

0.7161E+07  0.5761E+07  0.7306E+06  O.OOOOE+OO  O.OOOOE+OO  -0.5310E+07 

0.5761E+07  0.7161E+07  0.7306E+06  O.OOOOE+OO  O.OOOOE+OO  -0.5310E+07 

0.7306E+06  0.7306E+06  0.1554E+07  O.OOOOE+OO  O.OOOOE+OO  -0.1762E+05 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 

-0.5310E+07  -0.5310E+07  -0.1762E+05  O.OOOOE+OO  O.OOOOE+OO  0.5713E+07 

PLY:  5 

0.7161E+07  0.5761E+07  0.7306E+06  O.OOOOE+OO  O.OOOOE+OO  -0.5310E+07 

0.5761E+07  0.7161E+07  0.7306E+06  O.OOOOE+OO  O.OOOOE+OO  -0.5310E+07 

0.7306E+06  0.7306E+06  0.1554E+07  O.OOOOE+OO  O.OOOOE+OO  -0.1762E+05 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 

-0.5310E+07  -0.5310E+07  -0.1762E+05  O.OOOOE+OO  O.OOOOE+OO  0.5713E+07 

PLY:  6 

0.7161E+07  0.5761E+07  0.7306E+06  O.OOOOE+OO  O.OOOOE+OO  0.5310E+07 

0.5761E+07  0.7161E+07  0.7306E+06  O.OOOOE+OO  O.OOOOE+OO  0.5310E+07 

0.7306E+06  0.7306E+06  0.1554E+07  O.OOOOE+OO  O.OOOOE+OO  0.1762E+05 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 

0.5310E+07  0.5310E+07  0.1762E+05  O.OOOOE+OO  O.OOOOE+OO  0.5713E+07 

PLY:  7 

0.2230E+08  -0.7333  -2.363  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

-0.7333  -0.4024E-05  -2.222  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

-2.363  -2.222  -4.938  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06 
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PLY:  8 

0.2230E+08  -0.7333  -2.363  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

-0.7333  -0.4024E-05  -2.222  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

-2.363  -2.222  -4.938  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06 

ALPHA  GLOBAL  MATRIX  > 

PLY:  1 


-0.5000E-06  0.1550E-04 

0.1550E-04 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

PLY:  2 

-0.5000E-06  0.1550E-04 

0.1550E-04 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

PLY:  3 

0.7500E-05  0.7500E-05 

0.1550E-04 

O.OOOOE+00 

O.OOOOE+00 

-0.1600E-04 

PLY:  4 

0.7500E-05  0.7500E-05 

0.1550E-04 

O.OOOOE+00 

O.OOOOE+00 

0.1600E-04 

PLY:  5 

0.7500E-05  0.7500E-05 

0.1550E-04 

O.OOOOE+00 

O.OOOOE+00 

0.1600E-04 

PLY:  6 

0.7500E-05  0.7500E-05 

0.1550E-04 

O.OOOOE+00 

O.OOOOE+00 

-0.1600E-04 

PLY:  7 

-0.5000E-06  0.1550E-04 

0.1550E-04 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

PLY:  8 

-0.5000E-06  0.1550E-04 

0.1550E-04 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

CM  GLOBAL  STIFFNESS  MATRIX  > 
0.1456E+08  0.2709E+07  -4.684 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.2709E+07  0.3409E+07 

-4.543 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

-4.684  -4.543  -9.814  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E-K)6  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.7000E+06  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00  0.3206E+07 

FROM  OUTEFPR.F . 

EFFECTIVE  LAMINATE  PROPERTIES  > 

EX  =  1.2406259E-K)7 

EY  =  2904853. 

EZ  =  -9.814123 
NUYZ  =  0.3741388 
NUXZ  =  0.1094047 
NUXY  =  0.7946534 
NUZY  =-1.2640380E-06 
NUZX  =  -8.6545938E-08 
NUYX  =  0.1860634 
GYZ  =  699999.9 
GXZ  =  699999.9 
GXY  =  3206367. 

ALPHAX  = -4.0998452E-07 
ALPHAY  =  1.3799854E-06 
ALPHAZ  =  -1.771376 
ALPHAYZ=  O.OOOOOOOE+00 
ALPHAXZ=  O.OOOOOOOE+00 
ALPHAXY  =  -2.9743 140E-13 

FROM  OUTLOAD.F . 

AVERAGE  MECHANICAL  STRESSES  (PSI.)  > 

X  =  -1000.000 
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Y  =  O.OOOOOOOE+00 
Z  =  O.OOOOOOOE+00 
YZ  =  10.00000 

XZ  =  35.00000 

XY  =  10.00000 

AVERAGE  THERMAL  STRESSES  (PSl.)  > 

X-TH  =  -1213.327 
Y-TH  =  -6.7138564E-05 
Z-TH  =  7.7716759E-05 
YZ-TH  =  2.7808765E-06 
XZ-TH  =  3.972681  lE-12 
XY-TH  =  5.6752597E-18 

MECHANICAL  LAMINATE  STRAINS  > 

EPSX  =  -8.0604470E-05 
EPSY  =  6.40526 12E-05 
EPSZ  =  8.8185079E-06 
EPSYZ=  1.4285716E-05 
EPSXZ=  5.0000006E-05 
EPSXY=  3.1187942E-06 

THERMAL  LAMINATE  STRAINS  > 

EPSX-TH  =  -6.7138564E-05 
EPS  Y-TH  =  7.7716759E-05 
EPSZ-TH  =  2.7808765E-06 
EPSYZ-TH=  3.972681  lE-12 
EPSXZ-TH=  5.6752597E-18 
EPSXY-TH  =  1.7699968E-24 

FROM  OUTMESTR.F . 

MECHANICAL  GLOBAL  PLY  STRAIN 

PLY  EPS-X  EPS-Y  EPS-Z  EPS-YZ  EPS-XZ  EPS-XY 

1  -.806E-04  0.641E-04  0.975E-05  0.143E-04  0.500E-04  0.312E-05 

2  -.806E-04  0.641E-04  0.975E-05  0.143E-04  0.500E-04  0.312E-05 

3  -.806E-04  0.641E-04  0.774E-05  0.143E-04  0.500E-04  0.312E-05 

4  -.806E-04  0.641E-04  0.782E-05  0.143E-04  0.500E-04  0.312E-05 

5  -.806E-04  0.641E-04  0.782E-05  0.143E-04  0.500E-04  0.312E-05 

6  -.806E-04  0.641E-04  0.774E-05  0.143E-04  0.500E-04  0.312E-05 

7  -.806E-04  0.641E-04  0.975E-05  0.143E-04  0.500E-04  0.312E-05 

8  -.806E-04  0.641E-04  0.975E-05  0.143E-04  0.500E-04  0.312E-05 

MECHANICAL  GLOBAL  PLY  STRESS 

PLY  SIG-X  SIG-Y  SIG-Z  SIG-YZ  SIG-XZ  SIG-XY 

1  -.180E+04  0.375E-04  O.OOOE+00  10.0  35.0  2.18 

2  -.180E+04  0.375E-04  O.OOOE+00  10.0  35.0  2.18 


3 

-186. 

16.5 

O.OOOE+00 

10.0 

35.0 

-69.9 

4 

-219. 

-16.5 

O.OOOE+00 

10.0 

35.0 

106. 

5 

-219. 

-16.5 

O.OOOE+00 

10.0 

35.0 

106. 

6 

-186. 

16.5 

O.OOOE+00 

10.0 

35.0 

-69.9 

7  -,180E+04  0.375E-04  O.OOOE+00  10.0  35.0  2.18 

8  -.180E+04  0.375E-04  O.OOOE+00  10.0  35.0  2.18 

MECHANICAL  PRINCIPAL  PLY  STRAIN 

PLY  EPS-1  EPS-2  EPS-3  EPS-23  EPS-13  EPS-12 

1  -.806E-04  0.641E-04  0.975E-05  0.143E-04  0.500E-04  0.312E-05 

2  -.806E-04  0.641E-04  0.975E-05  0.143E-04  0.500E-04  0.312E-05 

3  -.672E-05  -.984E-05  0.774E-05  -.253E-04  0.455E-04  0.145E-03 
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4  -.984E-05  -.672E-05  0.782E-05 

5  -.984E-05  -.672E-05  0.782E-05 

6  -.672E-05  -.984E-05  0.774E-05 

7  -.806E-04  0.641E-04  0.975E-05 

8  -.806E-04  0.641E-04  0.975E-05 


0.455E-04  0.253E-04  -.145E-03 
0.455E-04  0.253E-04  -.145E-03 
-.253E-04  0.455E-04  0.145E-03 
0.143E-04  0.500E-04  0.312E-05 
0.143E-04  0.500E-04  0.312E-05 


MECHANICAL  PRINCIPAL  PLY  STRESS 
PLY  SIG-1  SIG-2  SIG-3  SIG-23 

1  -.180E+04  0.375E-04  O.OOOE+00  10.0 

2  -.I80E+04  0.375E-04  O.OOOE+00  10.0 

3  -155.  -14.8  O.OOOE+00  -17.7 

4  -223.  -12.2  O.OOOE+00  31.8 

5  -223.  -12.2  O.OOOE+00  31.8 

6  -155.  -14.8  O.OOOE+00  -17.7 

7  -.180E+04  0.375E-04  O.OOOE+00  10.0 

8  -.180E+04  0.375E-04  O.OOOE+00  10.0 


SIG-13  SIG-12 

35.0  2.18 

35.0  2.18 

31.8  101. 

17.7  -101. 

17.7  -101. 

31.8  101. 

35.0  2.18 

35.0  2.18 


FROM  OUTTHSTR.F . 

THERMAL  GLOBAL  PLY  STRAIN 


PLY 

EPS-X 

EPS-Y 

EPS-Z 

1 

-.671E-04 

0.777E-04 

-.186E-04 

2 

-.671E-04 

0.777E-04 

-.186E-04 

3 

-.671E-04 

0.777E-04 

-.497E-05 

4 

-.671E-04 

0.777E-04 

-.497E-05 

5 

-.671E-04 

0.777E-04 

-.497E-05 

6 

-.671E-04 

0.777E-04 

-.497E-05 

7 

-.671E-04 

0.777E-04 

-.186E-04 

8 

-.671E-04 

0.777E-04 

-.186E-04 

EPS-YZ  EPS-XZ  EPS-XY 
0.397E-11  0.568E-17  0.177E-23 
0.397E-11  0.568E-17  0.177E-23 
0.397E-11  0.568E-17  0.177E-23 
0.397E-11  0.568E-17  0.177E-23 
0.397E-11  0.568E-17  0.177E-23 
0.397E-11  0.568E-17  0.177E-23 
0.397E-11  0.568E-17  0.177E-23 
0.397E-11  0.568E-17  0.177E-23 


THERMAL  GLOBAL  PLY  STRESS 

PLY  SIG-X  SIG-Y  SIG-Z  SIG-YZ  SIG-XZ  SIG-XY 

1  -.373E+04  -.673E-02  0.777E-04  0.278E-05  0.397E-11  0.124E-17 

2  -.373E+04  -.673E-02  0.777E-04  0.278E-05  0.397E-11  0.124E-17 

3  0.462E+04  0.482E+04  0.777E-04  0.278E-05  0.397E-11  -.224E+04 

4  0.462E+04  0.482E+04  0.777E-04  0.278E-05  0.397E-11  0.224E+04 

5  0.462E+04  0.482E+04  0.777E-04  0.278E-05  0.397E-11  0.224E+04 

6  0.462E+04  0.482E+04  0.777E-04  0.278E-05  0.397E-11  -.224E+04 

7  -.373E+04  -.673E-02  0.777E-04  0.278E-05  0.397E-11  0.124E-17 

8  -.373E+04  -.673E-02  0.777E-04  0.278E-05  0.397E-11  0.124E-17 


THERMAL 
PLY  EPS- 
1  -.671E 

-.671E 
0.529E 
0.529E 
0.529E 
0.529E 
-.671E-I 
-.671E-I 


PRINCIPAL  PLY  STRAIN 


2 

3 

4 

5 

6 

7 

8 


1 

04 

04 

-05 

-05 

-05 

-05 

04 

04 


EPS-2 

0.777E-04 

0.777E-04 

0.529E-05 

0.529E-05 

0.529E-05 

0.529E-05 

0.777E-04 

0.777E-04 


EPS-3 

-.186E-04 

-.186E-04 

-.497E-05 

-.497E-05 

-.497E-05 

-.497E-05 

-.186E-04 

-.186E-04 


EPS-23 
0.397E-11 
0.397E-11 
0.281E-11 
0.281E-11 
0.28  IE- 11 
0.281E-11 
0.397E-11 
0.397E-11 


EPS-13 
0.568E-17 
0.568E-17 
0.28  IE- 11 
-.28  IE- 11 
-.281E-11 
0.281E-11 
0.568E-17 
0.568E-17 


EPS-12 

0.177E-23 

0.177E-23 

0.145E-03 

-.145E-03 

-.145E-03 

0.145E-03 

0.177E-23 

0.177E-23 


THERMAL  PRINCIPAL  PLY  STRESS 

PLY  SIG-1  SIG-2  SIG-3  SIG-23  SIG-13  SIG-12 

1  -.373E+04  -.673E-02  0.777E-04  0.278E-05  0.397E-11  0.124E-17 

2  -.373E+04  -.673E-02  0.777E-04  0.278E-05  0.397E-11  0.124E-17 

3  0.248E+04  0.696E+04  0.777E-04  0.197E-05  0.197E-05  101. 

4  0.248E+04  0.696E+04  0.777E-04  0.197E-05  -.197E-05  -101. 

5  0.248E+04  0.696E+04  0.777E-04  0.197E-05  -.197E-05  -101. 
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6 

7 

8 


0.248E+04 

-.373E+04 

-.373E+04 


0.696E+04  0.777E-04  0.197E-05  0.197E-05  101. 

-.673E-02  0.777E-04  0.278E-05  0.397E-11  0.124E-17 
-.673E-02  0.777E-04  0.278E-05  0.397E-11  0.124E-17 


FROM  OUTFAIL.F . 

ITERATION  NUMBER  >  2 

FAILURE  ASSESSMENT- 
PLY  >  1  MODE  >  XIC 

PLY  >  1  MODE  >  X2T 

PLY  >  2  MODE  >  XIC 

PLY  >  2  MODE  >  X2T 

PLY  >  7  MODE  >  XIC 

PLY  >  7  MODE  >  X2T 

PLY  >  8  MODE  >  XIC 

PLY  >  8  MODE  >  X2T 

SAFTY  FACTOR  >  95.28498 


FROM  OUTCIJS.F 


CIJ  PRINCIPAL  MATRIX 

PLY:  1 

0.4946 

-0.3627 

-0.3627 

0.2660 

-1.169 

-1.365 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

PLY:  2 

0.4946 

-0.3627 

-0.3627 

0.2660 

-1.169 

-1.365 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

PLY:  3 

0.2279E+08 

0.7482E+06 

0.7482E+06 

0.1554E+07 

0.7482E+06 

0.7130E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

PLY:  4 

0.2279E+08 

0.7482E+06 

0.7482E+06 

0.1554E+07 

0.7482E+06 

0.7130E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

PLY:  5 

0.2279E+08 

0.7482E+06 

0.7482E+06 

0.1554E+07 

0.7482E+06 

0.7130E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

PLY:  6 

0.2279E+08 

0.7482E+06 

0.7482E+06 

0.1554E+07 

0.7482E+06 

0.7130E+06 

> 


-1.169 

-1.365 

-2.177 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 


-1.169 

-1.365 

-2.177 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7482E+06 

0.7130E+06 

0.1554E+07 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7482E+06 

0.7130E+06 

0.1554E+07 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0J482E+06 

0.7130E+06 

0.1554E+07 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

0.7482E+06 

0.7130E+06 

0.1554E+07 


O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.7000E+06  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  0.7000E+06  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+OO  0.7000E+06 


O.OOOOE+00  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  0.7000E+06 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 
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O.OOOOE+00 
O.OOOOE+00 
O.OOOOE+00 
PLY:  7 

0.4946 
-0.3627 
-1.169 
O.OOOOE+00 
O.OOOOE+00 
O.OOOOE+00 
PLY:  8 

0.4946 
-0.3627 
-1.169 
O.OOOOE+00 
O.OOOOE+00 
O.OOOOE+00 


O.OOOOE+00 

O.OOOOE+OO 

O.OOOOE+00 

-0.3627 

0.2660 

-1.365 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

-0.3627 

0.2660 

-1.365 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

-1.169 

-1.365 

-2.177 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

-1.169 

-1.365 

-2.177 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 


-1.169 

-1.365 

-2.177 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

-1.169 

-1.365 

-2.177 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7306E+06 

0.7306E+06 

0.1554E+07 

O.OOOOE+OO 

O.OOOOE+OO 

0.1762E+05 

0.7306E+06 

0.7306E+06 

0.1554E+07 

O.OOOOE+OO 

O.OOOOE+OO 

-0.1762E+05 

0.7306E+06 

0.7306E+06 

0.1554E+07 

O.OOOOE+OO 

O.OOOOE+OO 

-0.1762E+05 

0.7306E+06 

0.7306E+06 

0.1554E+07 


0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 


O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06 


0.5310E+07 

0.5310E+07 

0.1762E+05 

O.OOOOE+OO 

O.OOOOE+OO 

0.5713E+07 

-0.5310E+07 

-0.5310E+07 

-0.1762E+05 

O.OOOOE+OO 

O.OOOOE+OO 

0.5713E+07 

-0.5310E+07 

-0.5310E+07 

-0.1762E+05 

O.OOOOE+OO 

O.OOOOE+OO 

0.5713E+07 

0.5310E+07 

0.5310E+07 

0.1762E+05 


CIJ  GLOBAL 
PLY:  1 

0.4946 
-0.3627 
-1.169 
O.OOOOE+OO 
O.OOOOE+OO 
O.OOOOE+OO 
PLY:  2 

0.4946 
-0.3627 
-1.169 
O.OOOOE+OO 
O.OOOOE+OO 
O.OOOOE+OO 
PLY:  3 

0.7161E+07 
0.5761E+07 
0.7306E+06 
O.OOOOE+OO 
O.OOOOE+OO 
0.5310E+07 
PLY:  4 

0.7161E+07 
0.5761E+07 
0.7306E+06 
O.OOOOE+OO 
O.OOOOE+OO 
-0.5310E+07 
PLY:  5 

0.7161E+07 
0.5761E+07 
0.7306E+06 
O.OOOOE+OO 
O.OOOOE+OO 
-0.5310E+07 
PLY:  6 

0.7161E+07 
0.576  lE+07 
0.7306E+06 


MATRIX  > 

-0.3627 

0.2660 

-1.365 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

-0.3627 

0.2660 

-1.365 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.5761E+07 

0.7161E+07 

0.7306E+06 

O.OOOOE+OO 

O.OOOOE+OO 

0.5310E+07 

0.576  lE+07 
0.7161E+07 
0.7306E+06 
O.OOOOE+OO 
O.OOOOE+OO 
-0.5310E+07 

0.576 1E+07 
0.7161E+07 
0.7306E+06 
O.OOOOE+OO 
O.OOOOE+OO 
-0.5310E+07 

0.576  lE+07 
0.7161E+07 
0.7306E+06 


O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  0.7000E+06 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  0.7000E+06 


O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  0.7000E+06 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 
0.7000E+06  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOE+OO  0.7000E+06 


72 


O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.5310E+07 

0.5310E+07 

0.1762E+05 

PLY:  7 

0.4946 

-0.3627 

-1.169 

-0.3627 

0.2660 

-1.365 

-1.169 

-1.365 

-2.177 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

PLY:  8 

0.4946 

-0.3627 

-1.169 

-0.3627 

0.2660 

-1.365 

-1.169 

-1.365 

-2.177 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7000E+06  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  0.7000E+06  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  0.5713E+07 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
0.7000E+06  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  0.7000E+06  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  0.7000E+06 

O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 
0.7000E+06  O.OOOOE+00  O.OOOOE+00 
O.OOOOE+00  0.7000E+06  O.OOOOE+00 
O.OOOOE+00  O.OOOOE+00  0.7000E+06 


ALPHA  GLOBAL  MATRIX  > 
PLY:  1 


-0.5000E-06  0.1550E-04 
PLY:  2 

0.1550E-04 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

-0.5000E-06  0.1550E-04 
PLY:  3 

0.1550E-04 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.7500E-05  0.7500E-05 
PLY:  4 

0.1550E-04 

O.OOOOE+OO 

O.OOOOE+OO 

-0.I600E-04 

0.7500E-05  0.7500E-05 
PLY:  5 

0.1550E-04 

O.OOOOE+OO 

O.OOOOE+OO 

0.1600E-04 

0.7500E-05  0.7500E-05 
PLY:  6 

0.1550E-04 

O.OOOOE+OO 

O.OOOOE+OO 

0.1600E-04 

0.7500E-05  0.7500E-05 
PLY:  7 

0.1550E-04 

O.OOOOE+OO 

O.OOOOE+OO 

-0.1600E-04 

-0.5000E-06  0.1550E-04 
PLY:  8 

0.1550E-04 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

-0.5000E-06  0.1550E-04  0.1550E-04 

CM  GLOBAL  STIFFNESS  MATRIX  > 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0.3409E+07  0.2709E+07  -2.192  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

0.2709E+07  0.3409E+07  -2.388  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

-2.192  -2.388  -4.291  O.OOOOE+00  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+OO  O.OOOOE+00  0.7000E+06  O.OOOOE+00  O.OOOOE+00 

O.OOOOE+00  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.7000E+06  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.3206E+07 

FROM  OUTEFPR.F  . 

EFFECTIVE  LAMINATE  PROPERTIES  > 

EX  =  1256259. 

EY  =  1256259. 

EZ  =  -4.290609 

NUYZ  =  0.1506686 
NUXZ  =  6.8527453E-02 
NUXY  =  0.7946534 
NUZY  =  -5.1459142E-07 
NUZX  =  -2.3404780E-07 
NUYX  =  0.7946537 
GYZ  =  699999.8 
GXZ  =  699999.8 
GXY  =  3206367. 
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ALPHAX  =  9.5462826E-07 
ALPHAY  =-1.8207174E-07 
ALPHAZ  =  -4.051761 
ALPHAYZ  =  O.OOOOOOOE+00 
ALPHAXZ=  O.OOOOOOOE+00 
ALPHAXY  =  -2.9743 140E-13 

FROM  OUTLOAD.F . 

AVERAGE  MECHANICAL  STRESSES  (PSI.)  > 
X  =  -1000.000 
Y  =  O.OOOOOOOE+00 
Z  =  O.OOOOOOOE+00 
YZ  =  10.00000 

XZ  =  35.00000 

XY  =  10.00000 

AVERAGE  THERMAL  STRESSES  (PSI.)  > 
X-TH  =  -2328.330 
Y-TH  =  -1.6637220E-03 
Z-TH  =  1.4727962E-03 
YZ-TH  =  -2.1625223E-04 
XZ-TH  =  -3.0893 182E-10 
XY-TH  =  -4.4133 127E-16 

MECHANICAL  LAMINATE  STRAINS  > 

EPSX  =  -7.9601456E-04 
EPSY  =  6.3255552E-04 
EPSZ  =  5.4548957E-05 
EPSYZ=  1. 42857 17E-05 
EPSXZ  =  5.00000  lOE-05 
EPSXY=  3.1187942E-06 

THERMAL  LAMINATE  STRAINS  > 

EPSX-TH  =  -1.6637220E-03 
EPSY-TH  =  1.4727962E-03 
EPSZ-TH  =  -2.1625223E-04 
EPSYZ-TH  =  -3.0893 1 82E- 1 0 
EPSXZ-TH=  -4.4133127E-16 
EPSXY-TH  =  -1.3764214E-22 


FROM  OUTMESTR.F  . 

MECHANICAL  GLOBAL  PLY  STRAIN 


PLY 

EPS-X 

EPS-Y 

EPS-Z 

EPS-YZ 

EPS-XZ 

EPS-XY 

1 

-.796E-03 

0.633E-03 

0.307E-04 

0.143E-04 

0.500E-04 

0.312E-05 

2 

-.796E-03 

0.633E-03 

0.307E-04 

0.143E-04 

0.500E-04 

0.312E-05 

3 

-.796E-03 

0.633E-03 

0.768E-04 

0.143E-04 

0.500E-04 

0.312E-05 

4 

-.796E-03 

0.633E-03 

0.769E-04 

0.143E-04 

0.500E-04 

0.312E-05 

5 

-.796E-03 

0.633E-03 

0.769E-04 

0.143E-04 

0.500E-04 

0.312E-05 

6 

-.796E-03 

0.633E-03 

0.768E-04 

0.143E-04 

0.500E-04 

0.312E-05 

7 

..796E-03 

0.633E-03 

0.307E-04 

0.143E-04 

0.500E-04 

0.312E-05 

8 

-.796E-03 

0.633E-03 

0.307E-04 

0.143E-04 

0.500E-04 

0.312E-05 

MECHANICAL 

GLOBAL  PLY  STRESS 

PLY 

SIG-X 

SIG-Y 

SIG-Z 

SIG-YZ 

SIG-XZ 

SIG-XY 

1 

-.659E-03 

0.415E-03 

O.OOOE+00 

10.0 

35.0  2.18 

2 

-.659E-03 

0.415E-03 

O.OOOE+00 

10.0 

35.0  2.18 

3 

-.198E+04 

16.5 

O.OOOE+00 

10.0  35.0  -849. 

4 

-.202E+04 

-16.5 

O.OOOE+00 

10.0  35.0  884. 
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5 

-.202E+04 

-16.5 

O.OOOE+00 

10.0 

35.0 

884. 

6 

-.I98E+04 

16.5 

O.OOOE+00 

10.0 

35.0 

-849. 

7 

-.659E-03 

0.415E-03 

O.OOOE+00 

10.0 

35.0 

2.18 

8 

-.659E-03 

0.415E-03 

O.OOOE+00 

10.0 

35.0 

2.18 

MECHANICAL 

PRINCIPAL  PLY  STRAIN 

PLY 

EPS-1 

EPS-2 

EPS-3 

EPS-23 

EPS-13  EPS-12 

I 

-.796E-03 

0.633E-03 

0.307E-04 

0.143E-04 

0.500E-04  0.312E-05 

2 

-.796E-03 

0.633E-03 

0.307E-04 

0.143E-04 

0.500E-04  0.312E-05 

3 

-.802E-04 

-.833E-04 

0.768E-04 

-.253E-04 

0.455E-04  0.143E.02 

4 

-.833E-04 

-.802E-04 

0.769E-04 

0.455E-04 

0.253E-04  -.143E-02 

5 

-.833E-04 

-.802E-04 

0.769E-04 

0.455E-04 

0.253E-04  -.143E-02 

6 

-.802E-04 

-.833E-04 

0.768E-04 

-.253E-04 

0.455E-04  0.143E-02 

7 

-.796E-03 

0.633E-03 

0.307E-04 

0.143E-04 

0.500E-04  0.312E-05 

8 

-.796E-03 

0.633E-03 

0.307E-04 

0.143E-04 

0.500E-04  0.312E-05 

MECHANICAL 

PRINCIPAL  PLY  STRESS 

PLY 

SIG-1 

SIG-2 

SIG-3 

SIG-23 

SIG-13 

SIG-12 

I 

-.659E-03 

0.415E-03 

O.OOOE+00 

10.0 

35.0 

2.18 

2 

-.659E-03 

0.415E-03 

O.OOOE+00 

10.0 

35.0 

2.18 

3 

-.183E+04 

-135. 

O.OOOE+00 

-17.7 

31.8 

1000. 

4 

-.190E+04 

-132. 

O.OOOE+00 

31.8 

17.7 

***** 

5 

-.190E+04 

-132. 

O.OOOE+00 

31.8 

17.7 

***** 

6 

-.183E-(-04 

-135. 

O.OOOE+00 

-17.7 

31.8 

1000. 

7 

-.659E-03 

0.415E-03 

O.OOOE+00 

10.0 

35.0 

2.18 

8 

-.659E-03 

0.415E-03 

O.OOOE+00 

10.0 

35.0 

2.18 

FROM  OUTTHSTR.F 


THERMAL  GLOBAL  PLY  STRAIN 


PLY 

EPS-X 

EPS-Y 

EPS-Z  EPS-YZ 

EPS-XZ  EPS-XY 

1 

-.166E-02 

0.147E-02 

-.707E-03  -.309E-09 

-.441E-15 

-.138E-21 

2 

-.166E-02 

0.147E-02 

-.707E-03  -.309E-09 

-.441E-15 

-.138E-21 

3 

-.166E-02 

0.147E-02 

0.897E-04  -.309E-09 

-.441E-15 

-.138E-21 

4 

-.166E-02 

0.147E-02 

0.897E-04  -.309E-09 

-.441E-15 

-.138E-21 

5 

-.166E-02 

0.147E-02 

0.897E-04  -.309E-09 

-.441E-15 

-.138E-21 

6 

-.166E-02 

0.147E-02 

0.897E-04  -.309E-09 

-.441E-15 

-.138E-21 

7 

-.166E-02 

0.147E-02 

-.707E-03  -.309E-09 

-.441E-15 

-.138E-21 

8 

-.166E-02 

0.147E-02 

-.707E-03  -.309E-09 

-.441E-15 

-.138E-21 

THERMAL  GLOBAL  PLY 

STRESS 

PLY 

SIG-X 

SIG-Y 

SIG-Z  SIG-YZ 

SIG-XZ 

SIG-XY 

1 

-.533E-02 

-.141E-02 

0.147E-02  -.216E-03 

-.309E-09 

-.963E-16 

2 

-.533E-02 

-.141E-02 

0.147E-02  -.216E-03 

-.309E-09 

-.963E-16 

3 

0.129E+04 

0.568E+04  0.147E-02  -.216E-03  -.309E-09  -.331E+04 

4 

0.129E+04 

0.568E+04  0.147E-02  -.216E-03  -.309E-09  0.331E+04 

5 

0.129E-K)4 

0.568E+04  0.147E-02  -.216E-03  -.309E-09  0.331E+04 

6 

0.129E+04 

0.568E+04  0.147E-02  -.216E-03  -.309E-09  -.331E-K)4 

7 

-.533E-02 

-.141E-02 

0.147E-02  -.216E-03 

-.309E-09 

-.963E-16 

8 

-.533E-02 

..141E-02 

0.147E-02  -.216E-03 

..309E-09 

-.963E-16 

THERMAL  PRINCIPAL  PLY  STRAIN 

PLY 

EPS-1 

EPS-2 

EPS-3  EPS-23 

EPS- 13 

EPS- 12 

1 

-.166E-02 

0.147E-02 

-.707E-03  ■.309E-09 

-.441E-15 

-.138E-21 

2 

-.166E-02 

0.147E-02 

-.707E-03  -.309E-09 

-.441E-15 

-.138E-21 

3 

-.955E-04 

-.955E-04 

0.897E-04  -.218E-09 

-.218E-09 

0.314E-02 

4 

-.955E-04 

-.955E-04 

0.897E-04  -.218E-09 

0.218E-09 

-.314E-02 

5 

-.955E-04 

-.955E-04 

0.897E-04  -.218E-09 

0.218E-09 

-.314E-02 

6 

-.955E.04 

-.955E-04 

0.897E-04  -.218E-09 

-.218E-09 

0.314E-02 

75 


7  -.166E-02 

0.147E-02  -.707E-03 

-.309E-09 

-.441E-15 

-.138E-21 

8  -.166E-02 

0.147E-02  -.707E-03 

-.309E-09 

-.44  IE- 15 

-.138E-21 

THERMAL  PRINCIPAL  PLY  STRESS 

PLY  SIG-1 

SIG-2  SIG-3 

SIG-23 

SIG-13 

SIG-12 

1  -.533E-02 

-.141E-02  0.147E-02 

-.216E-03 

-.309E-09 

-.963E-16 

2  -.533E-02 

-.141E-02  0.147E-02 

-.216E-03 

-.309E-09 

-.963E-16 

3  179. 

0.680E+04  0.147E-02 

-.153E-03 

-.153E-03 

0.220E+04 

4  179. 

0.680E+04  0.147E-02 

-.153E-03 

0.153E-03 

-.220E+04 

5  179. 

0.680E+04  0.147E-02 

-.153E-03 

0.153E-03 

-.220E+04 

6  179. 

0.680E+04  0.147E-02 

-.153E-03 

-.153E-03 

0.220E+04 

7  -.533E-02 

-.141E-02  0.147E-02 

-.216E-03 

-.309E-09 

-.963E-16 

8  -.533E-02 

-.141E-02  0.147E-02 

-.216E-03 

-.309E-09 

-.963E-16 

FROM  OUTFAIL.F . 

ITERATION  NUMBER  >  3 

FAILURE  ASSESSMENT... 

PLY> 

1  MODE  >  XIC 

PLY> 

I  MODE  >  X2T 

PLY> 

2  MODE  >  XIC 

PLY> 

2  MODE  >  X2T 

PLY> 

3  MODE  >  S12 

PLY> 

4  MODE  >  S12 

PLY> 

5  MODE  >  S12 

PLY> 

6  MODE  >  S12 

PLY> 

7  MODE  >  XIC 

PLY> 

7  MODE  >  X2T 

PLY> 

8  MODE  >  XIC 

PLY> 

8  MODE  >  X2T 

SAFTY  FACTOR  >  17.80445 


LOCAL  LAMINATE  STRAIN  EXCEEDS  TOLERANCE... 


FAE.URE  ANALYSIS  SUMMARY: 


ITERATIONS  =  4 

LOAD  DROPS  =  1 

SF  =  95.28498 

MODE  =  2 

PLY  =  8.000000 
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COPIES  ORGANIZATION 

1  HQDA  SARD  TT  DR  F  MILTON 
WASHINGTON  DC  20310-0103 

1  HQDA  SARD  TT  MR  J  APPEL 
WASHINGTON  DC  20310-0103 


1  HQDA  SARD  TT  MS  C  NASH 
WASHINGTON  DC  20310-0103 


1  HQDA  SARD  TR  DR  R  CHAIT 

WASHINGTON  DC  20310-0103 

1  HQDA  SARD  TR  MS  K  KQMINOS 

WASHINGTON  DC  201310-0103 

1  DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRL  CP  CA  D  SNIDER 
2800  POWDER  MILL  RD 
ADELPHI MD  20783 

6  DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRL  MAP 
L  JOHNSON 
BHALPIN 
TCHOU 
AMSRL  MA  PA 

D  GRANVILLE 
•  W  HASKELL 

AMSRL  MA  MA  G  HAGNAUER 
WATERTOWN  M A  02172-0001 

1  COMMANDER 

US  ARMY  ARMAMENT  RESEARCH, 
DEVELOPMENT  &  ENGINEERING  CTR 
ATTN  SMCARFSE  T  GORA 
PICATINNY  ARSENAL  NJ  07806-5000 

3  COMMANDER 

US  ARMY  ARMAMENT  RESEARCH, 
DEVELOPMENT  &  ENGINEERING  CTR 
ATTN  SMCARTD 
R  PRICE 
V  LINDNER 
C  SPINELLI 

PICATINNY  ARSENAL  NJ  07806-5000 

1  US  ARMY  TACOM 

ATTN  AMSTA  JSK  SAM  GOODMAN 
WARREN  MI  48397-5000 


NO.  OF 

COPIES  ORGANIZATION 

1  COMMANDER 

US  ARMY  ARMAMENT  RESEARCH 
&  DEVELOPMENT  CTR 
ATTN  F  MCLAUGHLIN 
PICATINNY  ARSENAL  NJ  07806 

5  COMMANDER 

US  ARMY  ARMAMENT  RESERACH 
DEVELOPMENT  &  ENGINEERING  CTR 
ATTN  SMCAR  CCH  T 
S  MUSALLI 
P  CHRISTIAN 
R  CARR 
N  KRASNOW 

PICATINNY  ARSENAL  NJ  07806-5000 

1  COMMANDER 

US  ARMY  ARMAMENT  RESEARCH 
DEVELOPMENT  &  ENGINEERING  CTR 
ATTN  SMCAR  CCH  V  E  FENNELL 
PICATINNY  ARSENAL  NJ  07806-5000 

1  COMMANDER 

US  ARMY  ARMAMENT  RESEARCH 
DEVELOPMENT  &  ENGINEERING  CTR 
ATTN  SMCAR  CCH  J  DELORENZO 
PICATINNY  ARSENAL  NJ  07806-5000 

2  COMMANDER 

US  ARMY  ARMAMENT  RESEARCH 
DEVELOPMENT  &  ENGINEERING  CTR 
ATTN  SMCAR  CC 
JHEDDERICH 
COL  SINCLAIR 

PICATINNY  ARSENAL  NJ  07806-5000 

1  COMMANDER 

US  ARMY  ARMAMENT  RESEARCH 
DEVELOPMENT  &  ENGINEERING  CTR 
ATTN  SMCAR  CCH  P  J  LUTZ 
PICATINNY  ARSENAL  NJ  07806-5000 

2  COMMANDER 

US  ARMY  ARMAMENT  RESEARCH 
DEVELOPMENT  &  ENGINEERING  CTR 
ATTN  SMCAR  FSAM 
D  DEMELLA 
F  DIORIO 

PICATINNY  ARSENAL  NJ  07806-5000 
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1  COMMANDER 

US  ARMY  ARMAMENT  RESEARCH 
DEVELOPMENT  &  ENGINEERING  CTR 
ATTN  SMCARFSA 
A  WARNASH 
BMACHAK 

PICATINNY  ARSENAL  NJ  07806-5000 

10  DIRECTOR 

BENET  LABORATORIES 
ATTN  SMCARCCB 
C  KITCHENS 
J KEANE 
J  BATTAGLIA 
JVASILAKIS 
G  FFIAR 
T  SIMKINS 
V  MONTVORI 
J  WRZOCHALSKI 
RHASENBEIN 
SMCAR  CCB  R  S  SOPOK 
WATERVLIETNY  12189 

1  COMMANDER 

ATTN  SMCWV  QAE  Q  C  HOWD 
BLDG  44  WATERVLIET  ARSENAL 
WATERVLIETNY  12189-4050 

1  COMMANDER 

ATTN  SMCWV  SPM  T  MCCLOSKEY 
BLDG  25  3  WATERVLIET  ARSENAL 
WATERVLIETNY  12189-4050 

1  COMMANDER 

WATERVLIET  ARSENAL 
ATTN  SMCWV  QA  QS  K  INSCO 
WATERVLIETNY  12189-4050 

1  COMMANDER 

PRODUCTION  BASE  MODERN  ACTY 
US  ARMY  ARMAMENT  RESEARCH 
DEVELOPMENT  &  ENGINEERING  CTR 
ATTN  AMSMCPBMK 
PICATINNY  ARSENAL  NJ  07806-5000 

1  COMMANDER 

US  ARMY  BEL  VOIR  RD&E  CTR 
ATTN  STRBEJBC 
FTBELVOIRVA  22060-5606 

1  US  ARMY  COLD  REGIONS  RESEARCH 

&  ENGINEERING  LABORATORY 
ATTN  P  DUTTA 
72  LYME  RD 
HANOVER  NH  03755 


1  DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRLWTL  D  WOODBURY 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 145 

3  COMMANDER 

US  ARMY  MISSILE  COMMAND 
ATTN  AMSMI  RD  W  MCCORKLE 
AMSMI  RD  ST  P  DOYLE 
AMSMI  RD  ST  CN  T  VANDIVER 
REDSTONE  ARSENAL  AL  35898 

2  US  ARMY  RESEARCH  OFFICE 
ATTN  A  CROWSON 

J CHANDRA 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

3  US  ARMY  RESEARCH  OFFICE 

ENGINEERING  SCIENCES  DIV 
ATTN  R  SINGLETON 

G  ANDERSON 
KIYER 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

2  PROJECT  MANAGER 

SADARM 

PICATINNY  ARSENAL  NJ  07806-5000 

2  PROJECT  MANAGER 

TANK  MAIN  ARMAMENT  SYSTEMS 
ATTN  SFAEARTMA  COL  BREGARD 
KKIMKER 

PICATINNY  ARSENAL  NJ  07806-5000 

1  PROJECT  MANAGER 

TANK  MAIN  ARMAMENT  SYSTEMS 
ATTN  SFAEARTMAMD 
R  KOWALSKI 

PICATINNY  ARSENAL  NJ  07806-5000 

2  PEO  FIELD  ARTILLERY  SYSTEMS 
ATTN  SFAEFASPM 

D  ADAMS 
T  MCWILLIAMS 

PICATINNY  ARSENAL  NJ  07806-5000 

1  PEO  FIELD  ARTILLERY  SYSTEMS 

ATTN  SFAE  FAS  PM  H.  GOLDMAN 
PICATINNY  ARSENAL  NJ  07806 
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PROJECT  MANAGER  AFAS 
ATTN  GDELCOCO 
J  SHIELDS 

PICATINNY  ARSENAL  NJ  07806-5000 

NASA  LANGLEY  RESEARCH  CTR 
MS  266 

ATTN  AMSRLVS 
W  ELDER 
F  BARTLETT  JR 
HAMPTON  VA  23681-0001 

COMMANDER 

DARPA 

ATTN  J  KELLY 
B  WILCOX 
3701  N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

COMMANDER 

WRIGHT  PATTERSON  AIR  FORCE  BASE 
ATTN  WLFIV  A  MAYER 
DAYTON  OH  45433 

NAVAL  SURFACE  WARFARE  CTR 

DAHLGRENDIV 

CODE  G33 

DAHLGREN  VA  224488 

NAVAL  RESEARCH  LABORATORY 
CODE  6383 
ATTN  I WOLOCK 
WASHINGTON  DC  20375-5000 

OFFICE  OF  NAVAL  RESEARCH 
MECH  DIV  CODE  1 132SM 
ATTN  YAPA  RAJAPAKSE 
ARLINGTON  VA  22217 

NAVAL  ORDNANCE  STATION 
ADVANCED  SYSTEMS  TECHNOLOGY  BR 
ATTN  D  HOLMES 
CODE  20 11 

LOUISVILLE  KY  40214-5245 

DAVID  TAYLOR  RESEARCH  CTR 
SHIP  STRUCTURES  &  PROTECT  DEPT 
ATTN  J  CORRADO  CODE  1702 
BETHESDA  MD  20084 

DAVID  TAYLOR  RESEARCH  CTR 
ATTN  R  ROCKWELL 
W  PHYILLAIER 
BETHESDA  MD  20054-5000 


1  DEFENSE  NUCLEAR  AGENCY 

INNOVATIVE  CONCEPTS  DIV 
ATTN  DR  R  ROHR 
6801  TELEGRAPH  RD 
ALEXANDRIA  VA  22310-3398 

1  DR  FRANK  SHOUP 

EXPEDITIONARY  WARFARE  DIV  N85 
2000  NAVY  PENTAGON 
WASHINGTON  DC  20350-2000 

1  OFFICE  OF  NAVAL  RESEARCH 

ATTN  MR  DAVID  SIEGEL  3  5 1 
800  N  QUINCY  ST 
ARLINGTON  VA  22217-5660 

1  JOSEPH  H  FRANCIS 

NAVAL  SURFACE  WARFARE  CTR 
ATTN  CODE  G30 
DAHLGREN  VA  22448 

1  JOHN  FRAYSSE 

NAVAL  SURFACE  WARFARE  CTR 
ATTN  CODE  G33 
DAHLGREN  VA  22448 

1  LTC  JYUJI  D  HEWITT 

DEFENSE  NUCLEAR  AGENCY 
INNOVATIVE  CONCEPTS  DIV 
6801  TELEGRAPH  RD 
ALEXANDRIA  VA  22448 

1  CDR,  NAVAL  SEA  SYSTEMS  CMD 
ATTN  D  LIESE 

2531  JEFFERSON  DAVIS  HIGHWAY 
ARLINGTON  VA  22242-5160 

1  NAVAL  SURFACE  WARFARE 

MARY  E  LACY  CODE  D4 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448 

5  DIRECTOR 

LAWRENCE  LIVERMORE  NATIONAL 
LABORATORY 
ATTN  R  CHRISTENSEN 
S  DETERESA 
F  MAGNESS 
M  FINGER 
S.  GROVE 
PO  BOX  808 
LIVERMORE  CA  94550 
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DIRECTOR 

LOS  ALAMOS  NATIONAL  LAB 

ATTN  D  RABERN 

MEE  13  MS  J  576 

PO  BOX  1633 

LOS  ALAMOS  NM  87545 

LOS  ALAMOS  NATIONAL  LAB 
ATTN  J  REPPA  MS  F668 
PO  BOX  1663 
LOS  ALAMOS  NM  87545 

OAK  RIDGE  NATIONAL  LABORATORY 
ATTN  R  M  DAVIS 
PO  BOX  2008 

OAK  RIDGE  TN  37831-6195 
DIRECTOR 

SANDIA  NATIONAL  LABORATORIES 
APPLIED  MECHANICS  DEPARTMENT 
DIVISION  8241 
ATTN  C  ROBINSON 
G  BENEDETTI 
W  KAWAHARA 
K  PERANO 
D  DAWSON 
P  NIELAN 
PO  BOX  969 

LIVERMORE  CA  94550-0096 
BATTELLE 

ATTN  C  R  HARGREAVES 
505  KING  AVE 
COLUMBUS  OH  43201-2681 

PACIFIC  NORTHWEST  LABORATORY 
ATTN  M  SMITH 
PO  BOX  999 
RICHLAND  WA  99352 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN  M  MURPHY 
PO  BOX  808  L  282 
LIVERMORE  CA  94550 

DREXEL  UNIVERSITY 
ATTN  ALBERT  SD  WANG 
32ND  AND  CHESTNUT  STREETS 
PHILADELPHIA  PA  19104 


1  PENNSYLVANIA  STATE  UNIVERSITY 

ATTN  DAVID  W.  JENSEN 
223  N  HAMMOND 
UNIVERSITY  PARK  PA  16802 

1  PENNSYLVANIA  STATE  UNIVERSITY 

ATTN  RICHARD  MCNITT 
227  HAMMOND  BLDG 
UNIVERISTY  PARK  PA  16802 

1  PENNSYLVANIA  STATE  UNIVERSITY 

ATTN  RENATAS  ENGEL 
245  HAMMOND  BLDG 
UNIVERSITY  PARK  PA  16801 

1  PURDUE  UNIVERSITY 

SCHOOL  OF  AERO  &  ASTRO 
ATTN  CT  SUN 
W  LAFAYETTE  IN  47907-1282 

1  STANFORD  UNIVERSITY 

DEPARTMENT  OF  AERONAUTICS  AND 
AEROBALLISTICS  DURANT  BUILDING 
ATTN  S  TSAI 
STANFORD  CA  94305 

1  UCLA 

MANE  DEPT  ENGR IV 
ATTN  H  THOMAS  HAHN 
LOS  ANGELES  CA  90024-1597 

2  U  OF  DAYTON  RESEARCH  INSTITUTE 
ATTN  RAN  Y  KIM 

AJIT  K  ROY 

300  COLLEGE  PARK  AVE 
DAYTON  OH  45469-0 168 

1  UNIVERISTY  OF  DAYTON 
ATTN  JAMES  M  WHITNEY 
300  COLLEGE  PARK  AVE 
DAYTON  OH  45469-0240 

2  UNIVERSITY  OF  DELAWARE 

CTR  FOR  COMPOSITE  MATERIALS 
ATTN  J  GILLESPIE 
M  SANTARE 

201  SPENCER  LABORATORY 
NEWARK  DE  19716 


NORTH  CAROLINA  STATE  UNIV 
CIVIL  ENGINEERING  DEPARTMENT 
ATTN  W  RASDORF 
L  SPAINHOUR 
PO  BOX  7908 
RALEIGH  NC  27696-7908 
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UNIV  OF  ILLINOIS  AT  URBANA-  2 

CHAMPAIGN 

NATIONAL  CENTER  FOR  COMPOSITE 
MATERIALS  RESEARCH 
216  TALBOT  LABORATORY 
ATTN  J  ECONOMY 
104  S  WRIGHT  STREET 
URBANA  IL  61801  1 

UNIVERSITY  OF  KENTUCKY 
ATTN  LYNN  PENN 
763  ANDERSON  HALL 

LEXINGTON  KY  40506-0046  1 

THE  UNIV  OF  TEXAS  AT  AUSTIN 
CENTER  FOR  ELECTROMECHANICS 
ATTN  J  PRICE 

10100  BURNET  RD  1 

AUSTIN  TX  78758-4497 

UNIVERSITY  OF  UTAH 
DEPARTMENT  OF  MECHANICAL  AND 
INDUSTRIAL  ENGINEERING  1 

ATTN  S  SWANSON 
SALT  LAKE  CITY  UT  84112 

VIRGINIA  POLYTECHNICAL  INSTITUTE 
&  STATE  UNIVERSITY 

DEPT  OF  ESM  4 

ATTN  MICHAEL  W  HYER 

KENNETH  L  REIFSNIDER 
BLACKSBURG  VA  24061-0219 

DR  ANTHONY  J  VIZZINI 

DEPT  OF  AEROSPACE  ENGINEERING 

UNIVERSITY  OF  MARYLAND 

COLLEGE  PARK  MD  20742  1 

AAI  CORPORATION 
PO  BOX  126 

HUNT  VALLEY  MD  21030-0126 

1 

JOHN  HEBERT 
PO  BOX  1072 

HUNT  VALLEY  MD  21030-0126 

ARMTEC  DEFENSE  PRODUCTS  1 

ATTN  STEVE  DYER 
85  901  AVE  53 
PO  BOX  848 
COACHELLA  CA  92236 


ADVANCED  COMPOSITE  MATERIALS 
CORPORATION 
ATTN  P  HOOD 
J  RHODES 

1525  S  BUNCOMBE  RD 
GREER  SC  29651-9208 

DAN  DAKIN 
SAIC 

2200  POWELL  ST  STE  1090 
EMERYVILLE  CA  94608 

MILES  PALMER 
SAIC 

2109  AIR  PARK  RD  SE 
ALBUQUERQUE  NM  87106 

ROBERT  ACEBAL 
SAIC 

1225  JOHNSON  FERRY  RD  STE  100 
MARIETTA  GA  30068 

DR  GEORGE  CHRYSSOMALLIS 
SAIC 

3800  W  80th  STREET 
STE  1090 

BLOOMINGTON  MN  55431 

ALLIANT  TECHSYSTEMS  INC. 

ATTN  C  CANDLAND 
JBODE 
R  BECKER 
KWARD 
600  2ND  ST  NE 
HOPKINS  MN  55343-8367 

AMOCO  PERFORMANCE  PRODUCTS  INC 
ATTN  M  MICHNO  JR 
4500  MCGINNIS  FERRY  RD 
ALPHARETTA  GA  30202-3944 

APPLIED  COMPOSITES 
ATTN  W  GRISCH 
333  NORTH  SIXTH  ST 
ST  CHARLES  IL  60174 

BRUNSWICK  DEFENSE 
ATTN  T  HARRIS 
STE  410 

1745  JEFFERSON  DAVIS  HWY 
ARLINGTON  VA  22202 
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1  CHAMBERLAIN  MANUF  CORP 

RESEARCH  AND  DEVELOPMENT  DIV 

ATTN  M  TOWNSEND 

PO  BOX  2545 

550  ESTHER  ST 

WATERLOO  lA  50704 

I  CUSTOM  ANALYTICAL  ENGINEERING 

SYSTEMS  INC 
ATTN  A  ALEXANDER 
STAR  ROUTE  BOX  4A 
FLINTSTONEMD  21530 

1  PROJECTILE  TECHNOLOGY  INC 

515  GILES  ST 

HAVRE  DE  GRACE  MD  21078 

1  ALLEN  BOUTZ 

NOESIS  INC 
1 1  ION  GLEBE RD 
STE  250 

ARLINGTON  VA  22201-4795 

1  ARROW  TECH  ASSO 

1233  SHELBURNE  RD  STE  D  8 
SOUTH  BURLINGTON  VT  05403-7700 

7  CIVIL  ENGRRESERACH  FOUNDATION 

ATTN  H  BERNSTEIN  PRESIDENT 
C  MAGNELL 
K  ALMOND 
R  BELLE 
M  WILLETT 
EDELO 
B  MATTES 

1015  15TH  ST  NW  STE  600 
WASHINGTON  DC  20005 

1  NATIONAL  INSTITUTE  OF  STANDARDS 

AND  TECHNOLOGY 
STRUCTURE  &  MECHANICS  GROUP 
POLYMER  DIV  POLYMERS  RM  A209 
ATTN  GREGORY  MCKENNA 
GAITHERSBURG  MD  20899 

1  DUPONT  COMPANY 

COMPOSITES  ARAMID  FIBERS 
ATTN  S  BORLESKE  DEVELOP  MGR 
CHESNUT  RUN  PLAZA 
PO  BOX  80702 

WILMINGTON  DE  19880-0702 


1  GENERAL  DYNAMICS 

LAND  SYSTEMS  DIVISION 
ATTN  DBARTLE 
PO  BOX  1901 
WARREN  MI  48090 

3  HERCULES  INC 

ATTN  R  BOE 

F  POLICELLI 
J  POESCH 
PO  BOX  98 
MAGNA  UT  84044 

3  HERCULES  INC 

ATTN  G  KUEBELER 

J  VERMEYCHUK 
B  MANDERVILLE  JR 
HERCULES  PLAZA 
WILMINGTON  DE  19894 

1  HEXCEL 

ATTN  M  SHELENDICH 
11555  DUBLIN  BLVD 
PO  BOX  2312 
DUBLIN  CA  94568-0705 

1  lAP  RESEARCH  INC 

ATTN  A  CHALLITA 
2763  CULVER  AVE 
DAYTON  OH  45429 

3  INST  FOR  ADVANCED  TECHNOLOGY 

ATTN  T  KIEHNE 
HFAIR 
P  SULLIVAN 
4030  2  W  BRAKER  LN 
AUSTIN  TX  78759 

1  TECHNICAL  DIRECTOR  AEROPHYSICS 

INST  FOR  ADVANCED  TECHNOLOGY 
THE  UNIV  OF  TEXAS  AT  AUSTIN 
ATTN  DR  W  REINECKE 
4320-2  W  BRAKER  LN 
AUSTIN  TX  78759-5329 

1  INTEGRATED  COMPOSITE  TECH 

ATTN  H  PERKINSON  JR 
PO  BOX  397 

YORK  NEW  SALEM  PA  17371-0397 

1  INTERFEROMETRICS  INC 

ATTN  R  LARRIVA  VICE  PRESIDENT 
8150  LEESBURG  PIKE 
VIENNA  VA  22100 


t 


1 


c 


NO.  OF 

COPIES  ORGANIZATION 

1  KAMAN  SCIENCES  CORPORATION 

PO  BOX  7463 

COLORADO  SPRINGS  CO  80933 

1  PM  ADVANCED  CONCEPTS 
LORAL  VOUGHT  SYSTEMS 
ATTN  J TAYLOR 

PO  BOX  650003 
MS  WT21 

DALLAS  TX  76265-0003 

2  LORAL  VOUGHT  SYSTEMS 
ATTN  G  JACKSON 

KCOOK 

1701  W  MARSHALL  DR 
GRAND  PRAIRIE  TX  75051 

1  MR  JOE  BACKOFEN 

BRIGS  CO 

2668  PETERBOROUGH  ST 
HERDON  VA  22071-2443 

1  SOUTHWEST  RESERACH  INSTITUTE 

MR  JACK  RIEGEL 
ENG  &  MATERIAL  SCIENCES  DIV 
6220  CULEBRA  RD 
PO  DRAWER  28510 
SAN  ANTONIO  TX  78228-0510 

1  DR  LOUIS  ZERNOW 

ZERNOW  TECHNICAL  SERVICES 
425  W  BONITA  AVE  SUITE  208 
SAN  DIMAS  CA  91773 

1  DR  ROBERT  EICHELBERGER  CONS 

409  W  CATHERINE  ST 
BEL  AIR  MD  21014-3613 

1  DR  PEI  CHI  CHOU 

DYNA  EAST  CORPORATION 
3201  ARCH  ST 

PHILADELPHIA  PA  1 9 1 04-27 1 1 

2  MARTIN  MARIETTA  CORP 
ATTN  P  DEWAR 

L  SPONAR 

230  EAST  GODDARD  BLVD 
KING  OF  PRUSSIA  PA  19406 

2  OLIN  CORPORATION 

FLINCHBAUGHDIV 
ATTN  E  STEINER 
B  STEWART 
PO  BOX  127 
RED  LION  PA  17356 


NO.  OF 

COPIES  organization 

1  OLIN  CORPORATION 
ATTN  L  WHITMORE 
10101  9TH  ST  NORTH 
ST  PETERSBURG  FL  33702 

1  RENNSAELER  POLYTECHNIC  INST 

ATTN  R  B  PIPES 

PRESIDENT  OFC  PITTSBURGH  BLDG 
TROY  NY  12180-3590 

1  SPARTA  INC 
ATTN  J  GLATZ 

9455  TOWNE  CTR  DRIVE 
SAN  DIEGO  CA  92121-1964 

2  UNITED  DEFENSE  LP 
ATTN  P  PARA 

G  THOMAS 

1 107  COLEMAN  AVE  BOX  367 
SAN  JOSE  CA  95103 

2  US  ARMY  RESEARCH  LABORATORY 

UNIVERSITY  OF  DELAWARE 
CENTER  FOR  COMPOSITE  MATERIALS 
ATTN  T  BOGETTI 
BFINK 

NEWARK  DE  19716 

2  ADMINISTRATOR 

DEFENSE  TECHNICAL  INFO  CENTER 
ATTN  DTICDDA 
8725  JOHN  J  KINGMAN  RD  STE  0944 
FT  BEL  VOIR  VA  22060-6218 

1  DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRL  OP  SD  TL  (TECH  LIB) 
ADELPHI  MD  20783-1145 

1  DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRL  OP  SD  TP  (TECH  PUB  OFC) 
ADELPHI  MD  20783-1145 

1  DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRL  OP  SD  TA  (REC  MGMT) 
ADELPHI  MD  20783-1145 
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ABERDEEN  PROVING  GROUND 
DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRL  OP  AP  L  (TECH  LIB) 

BLDG  305  APG-AA 

DIR  USARL 
ATTN  AMSRL  Cl 

C  MERMEGAN  394 
AMSRL  Cl  C 

WSTUREK1121 
AMSRL  Cl  CB 

R  KASTE  394 
AMSRL  Cl  S 

A  MARK  309 
AMSRL  SL  B 

P  DIETZ  328 
AMSRL  SL  BA 

J  WALBERT  1065 
AMSRL  SL  BL 
D  BELY  328 
AMSRL  SL  I 

DHASKILL  1065 
AMSRL  WT  P 

A  HORST  390A 
E  SCHMIDT  390A 
AMSRL  WT  PA 

T  MINOR  390 
C  LEVERITT  390 
D  KOOKER  390A 
AMSRL  WT  PB 

P  PLOSTINS  390 
AMSRL  WT  PC 

R  FIFER  390A 
AMSRL  WT  PD 

B  BURNS  390 
W  DRYSDALE  390 
J  BENDER  390 
R  MURRAY  390 
RKIRKENDALL390 
T  ERLINE  390 
D  HOPKINS  390 
S  WILKERSON  390 
D  HENRY  390 
R  KASTE  390 
L  BURTON  390 
J  TZENG  390 
R  LIEB  390 
G  GAZONAS  390 
M  LEADORE  390 
J  NEWILL  390 
C  HOPPEL  390 


ORGANIZATION 

AMSRL  WT  PD  ALC 
A  ABRAHAMIAN 
K BARNES 
M  BERMAN 
A  FRYDMAN 
TLI 

W  MCINTOSH 

E  SZYMANSKI  » 

AMSRL  WT  T 

W  MORRISON  309 

AMSRL  WT  TA  4 

W  GILLICH  390 
W  BRUCHEY  390 
AMSRL  WT  TC 

K  KIMSEY  309 
R  COATES  309 
W  DE  ROSSET  309 
AMSRL  WT  TD 

D  DIETRICH  309 
G  RANDERS  PEHRSON  309 
J  HUFFINGTON  309 
A  DAS  GUPTA  309 
J  SANTIAGO  309 
AMSRL  WT  W 

C  MURPHY  120 
AMSRL  WT  WA 

H  ROGERS  394 
B  MOORE  394 
AMSRL  WT  WB 

F  BRANDON  120 
W  D  AMICO  120 
AMSRL  WT  WC 

J  BORNSTEIN  120 
AMSRL  WT  WD 

A  ZIELINSKI  120 
J  POWELL  120 
AMSRL  WT  WE 

J  LACETERA  120 
J  THOMAS  394 
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1  MR  PETER  N  JONES 
DRA  FORT  HALSTEAD 
W7  DIVISION  BLDG  A20 
SEVENOAKS  KENT  TN  147BP 
UNITED  KINGDOM 

1  MR  FRANCOIS  LESAGE 
DEFENSE  RES  ESTAB  VALCARTIER 
PO  BOX  8800 

COURCELETTE  QUEBEC  COA 
IRO  CANADA 

2  ROYAL  MILIT  COLLEGE  OF  SCIENCE 
SHRTVENHAM 

ATTN  DR  DAVID  BULMAN 
DR  BRIAN  LAWTON 
SWINDON  WILTS  SN6  8LA 
UNITED  KINGDOM 

1  SWISS  FEDERAL  ARMAMENTS  WORKS 

ATTN  WALTER  LANZ 
ALLMENDSTRASSE  86 
3602  THUN 
SWITZERLAND 

1  PROFESSOR  SOL  BODNER 

ISRAEL  INSTITUTE  OF  TECHNOLOGY 
FACULTY  OF  MECHANICAL  ENG 
HAIFA  3200 
ISRAEL 

1  DSTO 

MATERIALS  RESEARCH  LABORATORY 
DR  NORBERT  BURMAN 
NAVAL  PLATFORM  VULNERABILITY 
SCHIP  STRUCTURES  &  MATERIALS  DIV 
PO  BOX  50 

ASCOT  VALE  VICTORIA 
AUSTRALIA  3032 

1  PRPFESSPR  EDWARD  CELENS 

ECOLE  ROYAL  MILITAIRE 
AVE  DE  LA  RENAISSANCE  30 
1040BRUXELLE 
BELGIQUE 

1  DEF  RES  ESTABLISHMENT  VALCARTIER 

ALAIN  DUPUIS 

2459  BOULEVARD  PIE  XI  NORTH 
VALCARTIER  QUEBEC 
CANADA 

PO  BOX  8800  COURCELETTE 
GOA  IRO  QUEBEC 
CANADA 


NO.  OF 

COPIES  ORGANIZATION 

1  INSTITUT  FRANCO  ALLEMAND  DE 

RECHERCHES  DE  SANIT  LOUIS 
DR  CLAUDE  FAUQUIGNON 
5  RUE  DU  GENERAL  CASSAGNOU 
BOITE  POSTALE  34 
F  68301  SAINT-LOUIS  CEDEX 
FRANCE 

1  INSTITUT  FRANCO  ALLEMAND  DE 

RECHERCHES  DE  SANIT  LOUIS 
DE  MARC  GIRAUD 
5  RUE  DU  GENERAL  CASSAGNOU 
BOITE  POSTALE  34 
F  68301  SAINT  LOUIS  CEDEX 
FRANCE 

1  TNO  PRINS  MAURITS  LABORATORY 

DR  ROB  IJSSELSTEIN 
LANGE  KLEIWEG  137 
PO  BOX  45 
2280  AA  RIJSWIJK 
THE  NETHERLANDS 

1  FOA  NATL  DEFENSE  RESEARCH  ESTAB 
DRBOJANZON 

DIR  DEPT  OF  WEAPONS  &  PROTECTION 

S  172  90  STOCKHOLM 

SWEDEN 

2  DEFENSE  TECH  &  PROC  AGENCY  GRD 
MR  GERHARD  LAUBE 

GENERAL  HERZOG  HAUS 
3602  THUN 
SWITZERLAND 

1  ROYAL  MIL  COLLEGE  OF  SCIENCE 

PROF  J  D  MACKWORTH 
SHRTVENHAM 
SWINDON  WILTS  SN6  8LA 
UNITED  KINGDOM 

1  MINISTRY  OF  DEFENCE 

RAFAEL 

DR  MEIR  MAYSELESS 
ARMAMENT  DEVELOPMENT  AUTHORITY 
PO  BOX  2250 
HAIFA  31021  ISRAEL 

1  DR  AKE  PERSSON 

DYNAMEC  RESEARCH  AB 
PARADISGRND  7 
S  151  36  SODERTALJE 
SWEDEN 
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1  ERNST  MACH  INSTITUT  EMI 

DR  GUSTAV  ADOLF  SCHRODER 
HAUPTSTRASSE  18 
79576  WEIL  AM  RHEIN 
GERMANY 

1  DRA 

DR  DAVE  SCOTT  TECH  MGR 

LAUNCH  SYSTEMS 

FT  HALSTEAD 

SEVENOAKS 

KENT  TN14  7BP 

ENGLAND 

1  ERNST  MACH  INSTITUT  EMI 

DR  ALOIS  STILP 
ECKERSTRASSE  4 
7800  FREIBURG 
GERMANY 

1  DR  IR  HANS  PASMAN 

TNO  DEFENSE  RESEARCH 
POSTBUS  6006 
2600  JA  DELFT 
THE  NETHERLANDS 

1  DR  BITAN  HIRSCH 

TACHKEMONY  ST  6 
NETAMUA  42611 
ISRAEL 

1  PROF  DR  MANFRED  HELD 

DEUTSCHE  AEROSPACE  AG 
DYNAMICS  SYSTEMS 
PO  BOX  1340 

D  86523  SCHROBENHAUSEN 
GERMANY 


